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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions, ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 

using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. V1, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, P 221, 1932). Ingeneral 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I1.A.U. practice. 


4- Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should a . Blocks are reta:ned by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alteration exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 
8. When submitting pres authors are requested to indicate whether they will be 


willing and able tu read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


9. Postcards giving the programme of each meeting are issued some days hefore the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 


“4 
| 
fee 
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OF THE 


ROYAL ASTRONOMICAL SOCIETY 
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MEETING OF 1956 JANUARY 13 


Professor Sir Harold Jetfreys, President, in the Chair 


The President announced the death of Professor W. M. H. Greaves, Past 
President of the Society and Astronomer-Royal for .Scotland, and read the 
following Resolution, passed by the Council that afternoon, the Fellows 
standing :— 

‘With the untimely death of William Michael Herbert Greaves, British 
astronomy has suffered a severe loss. Whilst he will be remembered chiefly 
for his pioneer work on stellar spectrophotometry, his contributions to dynamical 
astronomy and to geomagnetism were substantial. 

‘* He was always glad to serve the Society: he was a member of Council 
for an aggregate of 16 years, being Secretary for 7 years and holding the office 
of President from 1947 to 1949. Always an active member of the national 
and international committees on which he served, his considered advice will 
be greatly missed.” 


The President announced that the Council had awarded the Gold Medal 
of the Society to Professor Thomas George Cowling for his distinguished 
contributions to theoretical astrophysics and particularly for his work on the 
stability of stars. 


The election by the Council of the following Fellows was duly confirmed : 
Bruce Alan Bolt, M.Sc., Department of Applied Mathematics, University 
of Sydney, Australia (proposed by K. E. Bullen); 
*Michael William Walter Farrow, ‘Trinity College, Cambridge (proposed by 
H. Jeffreys); 
Victor Walter Raymond Charles Lionheart-Richardson, 213 Milkwood 
Road, Herne Hill, London, S.E.24 (proposed by J. F. H. Carr-Gregg); 
Richard Fred Pashley, M.A., 72 Brocklehurst Avenue, Norton, Sheffield 8 
(proposed by R. R. 5. Cox); 
Geoffrey ‘Turner, 43 Green Street, Sandbach, Cheshire (proposed by 
W. H. Steavenson); and 
*Roger D. Wade, 76 Starford Road, Hoddesdon, Hertfordshire (proposed by 
P. Moore). 


* Transferred from Junior Membership. 
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Meeting of 1956 January 13 Vol. 116 


Ninety-eight presents were announced as having been received since the 
last meeting, including :— 

O. C. Mohler, A table of solar spectrum wave-lengths (presented by the 
McMath-Hulbert Observatory); 

G. van den Bergh, Periodicity and variation of solar (and lunar) eclipses 
(presented by the Editor of Physics Abstracts); 

A. Dauvillier, Cosmologie et Chemie (presented by the author) ; 

W. J. Luyten, A catalogue of 1849 stars with proper motions exceeding 0-5 
annually (presented by the author); and 


V. de Callatay, Atlas du Ciel (presented by M. A. de Visscher). 
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SPECTROSCOPIC AND PHOTOMETRIC OBSERVATIONS 
OF S DORADUS 


A. F. Wesselink 


(Communicated by the Radcliffe Observer) 


(Received 1956 January 12 


Summary 


The P Cygni character of the spectrum of S Dor has been confirmed 
The radial velocity does not vary during the three years covered by the 
observations. The mean emission-line velocity is 295 km/sec and 
therefore the membership of the Large Magellanic Cloud is beyond 
question. The mean displacement of the absorption lines is 82 km/sec 
towards the violet; this is unusually small for a P Cygni star. The presence 
of helium is doubtful. The spectral type is Ao. ‘The Ca 1 lines are mainly 
stellar. Lines due to H and Catt are strong. For 74 lines wave-lengths 
and intensities are presented together with probable identifications. Photo- 
electric observations in an interval of seventeen months show a decline in 
brightness of o™-8 without change in colour. During this same interval the 
emission character of the spectrum became more pronounced, while for- 
bidden lines of Fe 11, which are diffuse, strengthened with respect to the 
permitted lines. In 1955 February the apparent photographic magnitude 
was 10-8; the corresponding absolute magnitude is —9™:2. S Dor is not 
an eclipsing variable but varies irregularly. 


1. The present paper presents wave-lengths and identifications of lines 
measured in the spectrum of S Dor with the 2-prism Cassegrain spectrograph 
attached to the 74-inch Radcliffe reflector. 


TABLE I 


Plate and Exp. Emulsion Region Radial velocity 
Dispersion Date ‘Time Kodak Measured Emission Absorption 
(min) km/sec km/sec 


c 404 1951 Nov. 9 60 103a0) - 3973-48660 297 {191 


c 418 1951 Nov. 13 75 " 3934-4924 296 196 
c1740 1953 Oct. 14 3892-4866 306 214 
b1763 1953 Oct. 29 126 re 4105-4866 296 227 
b1779 1953 Nov. 16 = 63 9 3930-4594 295 226 
c1796 1953 Nov. 29 120 1o3ak 3930-0568 315 22 
b1830 1954 Jan. 19 130 103a0) 3930-4929 293 224 
b1848 1954 Feb. 3 135 ws 3892-5023 293 228 
C2274 1955 Jan. 5 70 103aB 3913-0566 291 215 
C2279 1955 Jan. 8 go 103a0) 4104-4866 289 205 
C2281 1955 Jan. 11 35 - 3973-4866 285 196 
C2285 1955 Jan. 14 120 ms 3840-5023 283 214 
+ 295 213 


‘Table I presents the journal of the spectrographic observations. ‘lhe prefix 
to the plate number shows which camera was used ; the dispersions at Hy of the 
corresponding spectra are as follows: (b) 29, (c) 49 A/mm. ‘The various columns 
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are self-explanatory. 
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from emission and absorption lines respectively. 


Wave- 
length 
3537 °64 
3589°13 
3933 °66 
3905 
3970°21 
4101°7 
4122° 
4125" 
4173" 
4173" 
4233° 
4243" 
4257" 
4270" 
4287° 
4289" 
4296° 
4300°0: 
4303" 
4314" 
4340": 
4351" 
4353" 
4354" 
4358" 
4394° 
4413° 
44160" 
4443" 
4408-2 
4471" 
4472" 
4476" 
4451" 
4459: 
4490" 
4505" 
4515": 
4519" 
4522°0: 
4533" 
4541°23 
4549°45 
4555°55 
4555°89 
4571°95 
457012 
4550°34 
455385 
4555°52 
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The last two columns give the radial velocity as derived 


Mean wave-lengths and intensities of emission and absorption lines 


Type 


a 

ea,a 
ea 

¢a 

e 

ea,a 
ea,eé,a 
ea,e 
ea,e 
ea,e@,a 
ea,a 

e 

e 

a 
ea,eé,a 
ea,a 
ea 


ea,eé,a 
ea,a 
ea,a 
a, 
a 
e,a 
a 
e 
ea,a 
ea,eé,a 
ea,e 
ea,e 
ca,a 
ea,e,a 
ea,a 
ea,eé,a 
ea,e,a 
ca 
ea,a 
€,d 
ea,eé 
ea,a 


Intensity 


51 53/54 '55 


3 
1,1 


2,2 
6,4 


I, 
10,3 11,6 
as 
4,23 


2, 
I, 


11,6 


aN 


= 
~ 


w 


Main Contributor 
Element R.M.T. 
H 5°39; 
H 9°05; 
Car(s) 3°66; 
Cati(1) 8°47; 
H (1) 0°07; 

(1) 1°74; 
11 (28) 2°64; 
11 (27) 8°74; 
11 (27) 3°45; 
‘e 11 (28) 8°86; 

11 (27) 3°17; 
11 (21F) 3°98; 
Fe 11 (28) 8:16; 
fe 11 (21F) 6°83; 
Fe 11 (7F) 7°40; 
Ti 11 (41) 0°22; 
Fe 11 (28) 6°57; 
Ti 11 (41) 0°05; 
Fe 11 (27) 3°17; 
Fe 11 (32) 4°29; 
H 0°47; 
Fe 11 (27) 1°76; 
Fe 11(7F) 9°34; 
Fe 11 (27) 5°38; 
?He 1 (51) 7°93; 
Tit (19) 5°03; 
Fe 11 (7F) 3.78; 
Fe 11 (6F) 6°27; 
Ti (19) 3°80; 
Ti 11 (31) 8°49; 
?He 1 (14) 1°48; 
Fe 11 (37) 2°92; 


Mg 11 (4) 1°26; 
Fe 11 (37) 9°18; 
Fe 11 (37) 1°40; 
11 (38) 8°28; 
11 (37) 5°34; 
‘e 11 (37) 0°22; 
“e 11 (38) 2°63; 
"e 11 (37) 4°17; 
11 (38) 1°52; 
11 (38) 9°47; 
11 (37) 5°89; 
11 (44) 8°66; 
111 (82) 1°97; 
Fe 11 (38) 6°33; 
11 (26) 0°05; 
11 (38) 3°83; 
Cr 11 (44) 8°22; 


Minor Contributor 
Element R.M.T 


Fe 11 (3) 9°40 


Fe 11 (23F) 8-95 
Fe 11 (21F) 1°56 
Fe 11 (21F) 4°81 


Cr 11 (31) 5°57 


Ti 11 (41) 4.98 
Fe 11 (21F) 2.78 


Fe 11 (21F) 8°37 


Ti 11 (51) 4°06 
Fe 11 (32) 3.60 
Fe 11 (27) 6°82 


Fe 11 (6F) 8°75 


Ti 11 (50) 3°97 


Ti 11 (82) 9°62 
Cr 11 (44) 5°02 


° 
n 
3 
1,1 1,2 8 
2,2 6 
2,1 1,2 
6,6 5,3 I 
I, 
1,0 
38 4, 
83 
2,2 2,1 I, I 
2,2 
1,1 
I, 2, 
Goat : 1,1 
1,1 0,1 
e 
ca,eé I 
ea,e 
ea 1,1 
ea,a 
a 
0,0 i 
| 
I, »2 
2, 
0,0 
2,1 1,1 2,2 
2,3 1,3 
1,0 
22 
1,0 
38 
= 2.3 
1,0 I 
1,0 2 
2 
4:3 3:2 2 
2 


and photometric observations of S Doradus 


Intensity Main Contributor Minor Contributor 
length 53/54 Element R.M.T Element R.M.T. 
4592'07 ?Cr 11 (44) 2°09; 
4616°99 Cr 11 (44) 6°64; 
4620°12 ‘ “e 11 (38) ov51; 
4628-96 8 “e 11 (37) 9°34; 
4033 ? 11 (44) 4°11; 
4657°93 ?Fe 111 (3F) 
4666°23 oe 11 (37) 6°75; 
4730°89: 11 (43) 1°44; 
4814°16 2 ‘ “e 11 (20F) 4°55; 
486160 
492369 11 (42) 3°92; 
5017°85 11 (42) 8°43; 
5168-66 ? 11 (42) 9°03; 
5262°77 (19F) 1°6; 
5275°73 11 (49) 5°99; 
§317°14 11 (49) 6°61; Fe 11 (48) 6°78 
5407-06 ‘ell 6°94; unclassified 
5890°55 Na (1) 
§990°02_ 11 (46) 1°38; 
6330°'20 *e 11 (199) 1°97; 
6431°14 *e¢ 11 (40) 2°65; 
6493°79e ‘ell unclassified 
6516'10 e 11 (40) 6°05; 
6560°31 2:8: 


2. Table II presents the results of measurements in the spectrum of 5 Dor. 

First column: The mean wave-length after reduction for radial velocity. 
A colon indicates uncertainty in the measures. 

Second column: ‘The various types of the line that have been found on at 
least one plate: ea means a P Cygni line, thus showing emission and absorption, 
é means emission without absorption, whereas a indicates that the line occurs in 
absorption only. 

Third column: Mean estimated intensity. It was found difficult to represent 
the intensities of the lines for all the plates on the same scale because of the 
variations that took place with time. ‘lhe estimates on the various plates were 
adjusted so that those lines that were measured on all the plates showed no 
variation. ‘The intensities are given for three groups of plates taken respectively 
during the 1951-52, 1953-54 and 1954-55 oppositions. ‘lhe intensity of the 
emission line is always to the left of that of the absorption component, the two 


values being separated by a comma. ‘Though these intensities are only rough, 
they are usually sufficiently accurate to indicate the behaviour in time of a 
particular line relative to the other lines. 

Red-sensitive plates were only used in the second of the three intervals. The 


absence of intensities for lines with wave-lengths larger than 4900 A in the first 
and third intervals does not therefore necessarily mean these lines to be either 
weak or absent. 

Fifth column: Number of plates on which the line was recorded. 

Sixth column: Suggested identifications giving the responsible ion followed 
by the multiplet number in the Revised Multiplet Table and wave-length (also 
from the R.M.7.). Forbidden lines are distinguished by F following the multiplet 
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number. In the case of blends the probable major contributor is given first. 
The identifications were arrived at in the usual manner by systematic searches 
for the strongest lines of each ion in the R.M.T. 

3. Not only is the emission character of the spectrum on later Radcliffe plates 
more pronounced than on the earlier plates but there is evidence that [Fe] 
emission has strengthened relative to permitted Fe in the course of time. 
As an instance we mention the relative intensity of [Fe 1] 4244 and Fei 4233. 
The effect is clearly visible in Plate 1 where an early (c) and a late (d) spectrum 
are shown. It is very remarkable that the forbidden lines are definitely more 
diffuse than the permitted ones; so far we have been unable to find a satisfactory 
explanation for this fact. 

Dr A. D. Thackeray, in a recent study (1) of the spectrum of RR ‘Tel, finds 
the ratio of [Fe11] 4244 and Fe! 4233 to increase steadily with time while the 
emission lines of higher stages of ionization have been emerging. ‘There is no 
evidence that [Fe] lines in RR'Tel are more diffuse than Fett (as they are in 
S Dor) but increasing diffuseness in RR ‘Tel is shown by the lines of higher 
ionization. If the diffuseness of [Feit] lines in S Dor is to be attributed to 
turbulent broadening then one seems driven to conclude that the sharp Fe! 
lines originate in an entirely different atmospheric level relatively free from 
turbulence. Both the lines 4471 and 4388 are very weak and have provisionally 
been attributed to Het. ‘The presence of He is therefore rather doubtful and 
consequently the usual classification of the star as Be (2) does not seem to accord 
with Radcliffe spectra. It is, however, possible that the evidence for He has 
been stronger in the early Harvard spectra; for the present we classify the star 
as Ao. Our plates show the lines H and K of Catt in emission as well as in 
absorption. ‘This observation shows conclusively that the lines are mainly 
stellar. H. J. Smith (Harvard Thesis, 1955) has the same opinion about the 
origin of these lines, though he has found them only in absorption. His inference 
must be based on their great strength and the advanced spectral type of the 
tar. It is further remarkable that Smith finds He only in absorption whereas 
all the Radcliffe spectra show the line in absorption 1s well as in emission. 
\s a possible explanation of this difference in results we suggest the higher 
dispersion of our plates. ‘The line 4481 Mgit appears strong and in absorption 
only. In contrast with the other absorption lines 4481 is very diffuse and much 
less displaced. ‘This behaviour of 4481 seems to be similar in other stars of the 
P Cygni class (3). 

Miss Cannon (4) found 5 Dor to be a P Cygni star on objective prism plates. 
H. Shapley (5) then concluded, from the position of the star in the sky and from 
fact that galactic P Cygni stars are strongly concentrated towards the galactic 
plane, that 5 Dor was a member of the Large Magellanic Cloud. 

According to Beals (2) it is very difficult to determine the true radial velocity 
of a P Cygns star as the position of the lines is not uniquely determined by th 
velocity. ‘lhe intrinsic displacement is however smallest for the emission li 
he mean value of the radial velocity derived from the emission lines is 
+295 km/sec. ‘This result is sufficiently close to the radial velocity of the 
Large Magellanic Cloud at the position of S Dor to prove its membership 
without question. ‘lhe absorption-line velocities given in ‘Table I are systemati- 
cally higher when derived from spectra with the larger dispersion. ‘his is 
undoubtedly a consequence of the greater resolution obtained with the higl 
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dispersion. ‘lhe mean difference between the emission and absorption velocities 
given in ‘Table I is —82 km/sec. This is unusually small for a P Cygni star. 


TaAB.e ITI] 


Mean displacement 
with respect to 
emission lines 

km/sec 


Fe 
Cri 
H 

Call 
Mg 11 


4. In Table III the mean displacements of the absorption lines are given 
for the different ions. It is seen that there is little or no significant variation with 
the ion, with the exception of Mg for which the figure is due to 4481, already 
discussed in Section 3. Beals (2) found clear evidence of stratification in the 
atmosphere of P Cygni manifesting itself as a variation of the displacement 
with the ionization potential. Such a relation could not be found from the 
measures in the spectrum of 5 Doradus, probably because of the small range in 
ionization potential of the ions present. 

5. Photoelectric photometry.—On several occasions the brightness of 5S Dor 
was measured with the Cassegrain photoelectric photometer. A few remarks 
on this photometer seem to be called for. When the telescope is used for 
spectroscopy, the light of the star reaches the spectrograph slit direct from the 
Cassegrain mirror. In the case of Cassegrainian photoelectric photometry a 
flat mirror a few inches above the spectrograph slit deflects the light sideways 
into the photometer. In order to focus the stellar image on the diaphragm in 
the photometer it is necessary to move the Cassegrain mirror about half an inch 


I\ 
Brightness of S Dora 


Date 


30. 9.1953 2434051 

17. 9.1954 5003 9°50 9°75 
13.12.1954 5090 9°94 10°13 
19.12.1954 93 Q Ol 

6. 2.1955 5145 10°04 10°23 o'19g 


closer to the primary mirror. ‘lable IV shows the photographic and photovisual 
magnitudes, so far obtained, reduced to the Cape 5S system (5). ‘The brightness 
is seen to vary markedly during recent years ; the star is of photographic magnitude 
10™-o in 1955 February. ‘lhe colour-index (—o"-1g in the Cape 5 system) is 
seen to be remarkably constant. Its value in the international system is — 0-07 
After reduction for galactic reddening we obtain an intrinsic colour index of 

-o-1g in the international system; this is rather blue for an Ao star. The 
brightness of 10-0 photographic is about o"'-8 fainter than the faintest brightness 


S2 
96 
55 
62 
76 
21 
m 
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recorded by 5. Gaposchkin (7). Though part of this amount may be due to a 
difference in photometric system, it may also be caused by the light of neigh- 
bouring stars contributing to that of S Dor on the small-scale Harvard photo- 
graphs used by him. ‘The recent drop in brightness is inconsistent with the 
ephemeris given by Gaposchkin on the hypothesis that the star is an eclipsing 
variable. ‘The suggestion of Hubble and Sandage that S Dor belongs to the 
same class of bright blue variables as discovered by them in M31 and M33 
seems more likely (8). 


+ 
2:0 

9.0r-— 
all 
Mpg a 
= 


é 

243 4000 5000 


Fic. 1 
Ordinate: in upper part, the number of emission lines divided by the number of absorption lines; 
in lower part, the photographic magnitude. 
Abscissa; in either part, the Julian Date. 


The gradual fading of the light of S Dor in recent times is shown graphically 
in the lower half of Fig. r. The upper half of the diagram exhibits the ratio 
between the number of emission lines to the number of absorption lines measured 
on Radcliffe spectra. ‘The marked increase of the ratio with time could of course 
be attributed to the general fading of the continuum as shown by the photoelectric 
measures, the emission remaining constant; but other explanations seem to be 
possible. 


‘The author wishes to record his indebtedness to Dr A. D. ‘Thackeray and 
Dr M. W. Feast for invaluable advice. | am deeply indebted to H.M. Astronomer 
at the Cape, Dr R. H. Stoy, for the use of an amplifier and to Dr A. W. J. Cousins 
for much support in photoelectric matters. 
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RESULTS OF A SEARCH FOR BRIGHT 8 CEPHEI VARIABLES 
IN THE SOUTHERN SKY 


B. E. Pagel 


(Communicated by the Radcliffe Observer) 


(Received 1956 January 3) 


Summary 


The radial velocities of seven bright southern early B-type stars, known 
to have variable velocity, have been systematically examined at the Radcliffe 
Observatory in an effort to detect the presence of short periods (3 to 7 hours) 
characteristic of / Cephei variables. Positive results were obtained for three 
stars: #$ Cru (o<2K<14 km/sec), 7! Lup (2K <11 km/sec), and «a Lup 
(2K 15 km/sec). Definitive periods have been derived for 7' Lup and 
x Lup, while a provisional period is suggested for # Cru. 

( Cru and 7! Lup have been observed photoelectrically by the staff of the 
Cape Observatory and found to have short-period variations in light with a 
range of about o-o4 In the case of 7! Lup, the light variation shows the 
same period as the velocity variation, with the phase relationship typical of 
fb Cephei variables. In the case of / Cru, all that can be said is that the light 
variation has approximately the period of the velocity variation. «a Lup was 
observed photoelectrically by A. B. Muller at the Leiden Observatory 
Southern Station; no variation in excess of o™-o1 has so far been detected 


T' Lup fits closely the period—luminosity relations derived for / Cephei 
variables by Blaauw and Savedoff and by Petrie. It is clearly a typical raember 
of the group with constant velocity amplitude. a Lup comes tolerably near 


to the period—luminosity relation, but, until some variation in light can be 
established, one cannot quite exclude the possibility that it is a single-lined 
spectroscopic binary. / Cru, in spite of its early spectral type (Bo), would 
seem to be related to the / Cephei group, though it is anomalous in a number 
of ways 


1. Introduction.—\n recent years, there has been a revival of interest in 
variable stars of the 8 Cephei or B Canis Majoris group. ‘These are early 
B-type stars, somewhat above the main sequence in the HR diagram, which 
vary in light and radial velocity with periods between 3 and 7 hours. At present, 
ten stars visible from northern observatories are known to be members of this 
class (see, for example, Struve, 1955 b) and have been studied extensively. ‘They 
appear to fit a period-luminosity curve consistent with the normal period—density 
relation for pulsating stars (Blaauw and Savedoff, 1953; Petrie, 19544) and a 
period-colour relation (McNamara and Williams, 1955), and their variations in 
light and radial velocity can usually be described by one-sine or two-sine curves ; 
in each component oscillation, light maximum 1s closely coincident with the 
descent of the radial-velocity curve through its mean value. 

‘Together with these general features, there is an immense variety in individual 
behaviour. ‘Three stars with small variations in velocity and light—y Peg, 6 Cet 
and £4 CMa—behave in a comparatively simple fashion with a single period and 
apparently constant velocity amplitude (McNamara, 1955). On the other hand, 


| 
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BCep has an irregularly varying amplitude and five other members of the group 
have a well-marked secondary oscillation giving rise to a beat phenomenon 
(Struve, 1952). These latter, and BW Vul (Petrie, 1954b; McNamara, Struve 
and Berthiau, 1955; Odgers, 1955), show a periodic variation in line width 
—identified with high dispersion as doubling of the lines—together with 
departures from a sinusoidal velocity curve (e.g. a stillstand or intermediate 
branch on the descending portion) and a phase lag between the velocities 
determined from the hydrogen and other lines (Struve, 1955); in a few cases, 
a secular variation in period has been detected (Struve, 1955 a). 

The general features of the group suggest a radial pulsation analogous to 
that of cepheids and this hypothesis fits well with observations of 8 Cep in six 
colours (Stebbins and Kron, 1954). However, a similar investigation of BW Vul 
has been inconclusive (Walker, 1954a), and the pulsation hypothesis seems to 
break down completely for 16 Lac (Walker, 1954b). 

This variety of phenomena suggests a need for further efforts to identify 
and study new § Cephei variables. ‘The high time resolution required for the 
elucidation of details of the velocity curves sets a limiting magnitude to variables 
detectable spectroscopically with a given telescope and dispersion; hence a study 
of bright early B-type stars of known variable velocity in the southern sky 
immediately suggests itself. In this paper, we describe a spectrographic study 
of a few promising stars carried out at the Radcliffe Observatory between 1955 
February and 1955 July; in individual cases, the radial-velocity data were 
supplemented with photoelectric observations carried out by the staff of the 
Cape Observatory and by A. B. Muller at the Leiden Observatory Southern 
Station. 

2. Observational details.—'Vhe spectra were obtained with the two-prism 
Cassegrain spectrograph of the Radcliffe 74-inch reflector, using the a camera 
(f/8: dispersion 21 A/mm at Hy) with slit width o-o30mm (projected width 
o'o17mm). Emulsions were selected to give exposure times between 5 and 
10 minutes. Comparison spectra were obtained from an iron are exposed for 
equal periods at the beginning and end of each stellar exposure and the plates 
calibrated by means of a tube sensitometer. Lhe plates were measured using 
the Casella projection instrument of the Radcliffe Observatory, each set of 
measurements being repeated with the plate reversed. Lines were selected for 
measurement in accordance with current practice at the Radcliffe Observatory 
(Feast, ‘Thackeray and Wesselink, 1955)—-based on the I.A.U. standards and 
on the recommendations of Petrie (1953b)—but with a few modifications. 
Details of the lines measured and wave-lengths adopted in this investigation are 
given in ‘lable I. 

[he effective wave-lengths of the lines A4088 Sitv and A4471 Het were 
chosen empirically to bring about the best possible agreement between the 


radial velocities derived from these lines and from the others (for the wave-length 


cf. Petrie, 1953 b) Ihe spectral classes in the fourth column of 
| are the MKK classes of stars fe which the relevant line was found 

eful in this particular investigation and are not intended as general recommen- 
atiol 
(1899) in the form described by Wesselink (1952). ‘The standard diurnal 


he reductions to the Sun were carried out by the method of Schlesinger 


correction was also applied. 
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TABLE I 


Lines used for Radial velocity measurement 


Laboratory Identi- Adopted Range of Remarks 
wave-length fication wave-length spectral 
type 
3918°977 c 3915°977 B 2,3 } I.A.U. standards not in Petrie’s list 
3920°677 C iu 3920°677 B 2,3 
3933 °064 Ca tl 3933°664 B2,3 
3964°727 He 1 3964°727 all 
3970°075 H « 3970°O75 &* Cen only 
3994'°996 N 1 3994996 all 
4026°1%9 He! 4026-189 all Lab. wave-length seems satisfactory 
4069°794 Oi 4069°794 all 
4088-862 a 4088-998 Bo { Not in Petrie’s list. Lab. wave-length 
increased by 0:136A (10 km/sec) 
4101°738 H 0 4101°738 all 
4116°103 Si Iv 4116°103 Bo Not in Petrie’s list. 
4120°837 Het 4120°812 all I.A.U. wave-length used. 
4143°759 He 1 4143°759 all 
4267°167 Ci 4267°160 all 
4340°468 H y 4340°466 all 
43387°928 He! 4387°928 all 
4437°549 He 1 4437°549 all 1.A.U. standard not in Petrie’s list 
ane all { Lab. wave-length increased by 0-075 A 
(5 km/sec). 

4451°228 Mg 1 4481°228 B 2,3 
4552°662 Si ut 4552°622 all 
4507°841 Si 45607°841 all 


As a check on the general procedure of measurement, occasional plates of 
three 1.A.U, standard stars of late spectral type were taken and measured, using 
standard wave-lengths recommended by Petrie (1946), apart from the line 
A4494°575 Fei, which was found to be unreliable and therefore rejected. 


TABLe II 


Measured radial velocities of standard stars 


Star Spectral Date, U.T. Measured R.V. 
‘Type R.V. (Wilson) 
km/sec km) sec 
693 
(6 Cet) F 5 1955 July g‘181 12°9+ 1°0 14°8 


81797 March 


‘74 4°7 ‘4 
(a Hya) March 30°777 4°1+0'4 
April 7°822 
April 13°747 6-9+0'4 
May 18°700 
Vlay 21°702 
Mean O° 


150795 March I 2°7 
Tri A) May 1°174 
May 22°015 2°4+1°0 
Mean 2‘1+0°6 


* Taken with the + dispersion (29 A/mm) 


: % 
i 
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‘The results are given in Table II, where they are compared with those given in 
R. E. Wilson’s catalogue (R. E. Wilson, 1953). ‘lhe errors given in this and 
later tables are standard deviations per plate. 

The data are too few in number for a detailed analysis of the residuals, but 
it is satisfactory to note that they can be of either sign, and that they are of the 
order of magnitude of the uncertainty in the catalogue values (1-2 km/sec). 

3. Results of radial-velocity and photoelectric measurements.—The selection of 
promising stars for investigation was based on a list supplied by D. H. McNamara. 
‘They are early B-type stars noted as variable velocity stars in the Lick Observatory 
Catalogue (Lick Observatory, 1932). ‘The stars considered are given in ‘Table III, 
together with a summary of the results obtained from radial-velocity measure- 
ments. ‘The MKK spectral types given in the table have been taken from 
Cape Mimeogram 1 (Stoy, 1953). 


Tasie III 


Summary of results of radial-velocity measurements 


1950 Spectral No. No. Conclusions 
Type of of 
(MKK) plates nights 
28873 B 3s} Not observed owing to unsuitable 
(0 Cae) R.A. 


106490* 12" 12M 58° B2V 3 Lines too diffuse for velocity 
(0 Cru) measurements of adequate accuracy. 


111123 12% 44m- 5’ Bolll 79 Sharp lines: definite variation 

Cru) with period between o@-15 and 
o"-2; exact period doubtful (may be 
variable). 2K variable, sometimes 
reaching 14 km/sec. 


113791 13! 04™- Sharp lines: single-line spectro- 

(€? Cen) scopic binary (orbit by Neubauer, 
1931) with no velocity variation 
detectable during period of observ- 
ation, 


26341* 14! 22m. B 2 (IV)t Sharp lines: definite variation 
(7! Lup) with period o4-177365 + 15. 
2K 10°6 km/sec. 
129056 14" 38™. Sharp lines: definite variation 
(x Lup} with period 259882 4 20. 
2K = 14°7 km/sec. 


132200 14" 55m. 5 Sharp lines: no velocity variation 
(« Cen) detectable during period of observ- 
ations. 


134687 15" Double-line spectroscopic binary 

(0 Lup) with sharp and diffuse components 
Velocity of sharp component varies 
from night to night, but not on one 
night. 


* Added to McNamara’s list by A. D. Thackeray. 

+ Derived by the writer from a single low-dispersion spectrum. The luminosity class 
is uncertain. 

t H.D. spectral class ; MKK class not available. 
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The results indicate that the three stars 8 Cru, +! Lup and « Lup definitely 
possess a short-term variation in radial velocity, which suggests that they may 
be 6 Cephei variables of small amplitude (comparable with y Peg, ‘CMa and 
6Cet). The spectroscopic binaries €*Cen and Lup definitely are not, while 
the lines in the spectrum of 6 Cru are too diffuse for any conclusion to be drawn. 
«Cen merits further detailed study, since variability is suggested by the early 
Lick observations (Campbell and Moore, 1928), but a short period seems to be 
excluded by the present observational material. 

We now give the detailed results for individual stars. 

(i) €*Cen, «Cen and 6 Lup. ‘The measured radial velocities of these stars, 
together with the velocities in Wilson’s catalogue, are given in Table IV. 


TABLE IV 


Radial velocities of €* Cen, x Cen, and 6 Lup 


Mean Date, U.T. Observed 
Velocity Velocity 
(Wilson) 

14°3 April 30°927 +29'2+06 
‘949 04 

974 30°3t0°5 

‘999 +20'9+0°6 

5°742 +29'2+0°7 

‘768 |} 28:6+0°8 

"882 +26°3+1°0 


“149 3°9t 
"154 + 
‘069 +6°34 
‘109 
135 
"156 rs 
"174 +494 
963 
"930 
‘O31 +5 ‘8+ 
‘075 + 6-94 
B15 


4 Lup 
(sharp component) 


(ii) BCru. ‘The measured radial velocities of 8 Cru are presented graphically 
in Figs. 1, 2. The epochs were converted into heliocentric Julian days using 
Prager’s tables (Prager, 1932); as the exposure times were between 5 and 
10 minutes, the epochs are accurate to about 04-001. Considerable difficulty 
was experienced in determining a period. ‘The amplitude is variable (Fig. 1), 
and the apparent departures from strict periodicity may be either genuine or 
due to errors of observation in cycles when the amplitude was small. ‘The radial- 


§ Star 
&* Cen 
« Cen 1955 April o's 
I 0°7 
3 
3 o's 
3 
3 
3 
May 
I o's 
I 
I o'9 
1°3 
I 
I 
1955 April 4°137 
4°163 +18-st+1°2 
30°097 72+1°8 
30°120 77+1'8 
an 
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velocity observations have been supplemented by photoelectric observations 
carried out at the Cape Observatory on four nights when conditions were not 
considered good enough for the general photometric programmes. On one of 
these nights (1955 March 25=J.D. 2435192), both photoelectric and _ radial- 
velocity measurements were carried out; these data suggest that light maximum 
coincides approximately with radial velocity descending through y. If we 
assume that this always holds, we obtain a period P=o0"-160474 + 50, but this 
period is not established with any certainty. ‘he deviations are large (Fig. 2) 
and it has not proved possible to find any interfering oscillation that explains 
them. It is quite possible that the oscillations of 8 Cru are not strictly periodic ; 
neither the radial-velocity data nor the photometric data are sufficiently accurate 
or numerous to decide this point. When combined by means of the period 
given above, the radial-velocity observations scatter about a symmetrical 
curve; the sine curve obtained by least squares is 

V = + 12°44 + 2:20 sin (27t/0*- 160474), 

+0°25+0°50 (1) 

t=J.D.—(24352119-357 + 07-160474£), 

+ 04-006 (2) 
where E is an integer (Fig. 2). Owing to the considerable discrepancies in the 
phases, the amplitude in equation (1) is somewhat less than the mean of the 
observed amplitudes. ‘The r.m.s. residual is 2-2 km/sec. 


T 
ole 


Fic. 1 Observed radial ve of B Cru 
Abscissae: Heliocentric 7.D.—2435000 


Ordinates: Radial velocity 


locitu 


Vertical lines represent the observations, with the standard error estimated for each plate; the 


dotted curve represents the least-squares solution, assuming the period 0@-160474 


! 
| 735 45 55 204.40 60 
20r- {2 
r + + 4 
- 
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4 
| 20530 40 2u.30 4 22795 45 
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| } | 
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Cycle beginning JD 2435 188.949 

188.409 

192.261 

195,149 

195.310 

195.470 

212.390 

212.480 

hic. 2.—Radial-velocity and photoelectric observations of B Cru, combined into one cycle. 

Abscissae: Phase, in days. 

Ordinates: Upper figure, radial velocity. 

Lower figure, Am. 


The circles and crosses represent the observations; the full line in the upper figure represents the 
least-squares solution. The limitations of the solution are obvious. 


The results of the Cape photoelectric observations, combined into one cycle 
with the aid of equation (2), are illustrated in the lower half of Fig. 2. y Cen, 
5Cen and pCen were used as comparison stars and observations were carried 
out in blue and yellow light (on the natural colour system of the photometer). 
The ordinate in the diagram is Am, the mean of the departures in blue and yellow 
light from the respective mean magnitudes m,=0™-793 and m,=1"™-650; the 
uncertainty in the zero-points ranges from o™-o1 to 0-03. 

The scatter in the individual observations of this series is too large to make 
any definite deductions concerning the shape and periodicity of the light curve, 
but they are sufficient to demonstrate a range of about o™-o4 between maximum 
and minimum and a period between o"-15 and o%-20, in agreement with the 
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radial velocity observations. More definite conclusions must be postponed 
until a greater number of simultaneous photoelectric and spectrographi 
observations can be carried out on first-class photometric nights. 

No changes were observed in the appearance of the spectrum of 8 Cru. 

(11) 7! Lup. ‘The measured radial velocities of r' Lup are given in Figs. 3, 4. 
The observations show a well-defined periodicity with P= 04177365 +15, and, 
when combined into one cycle, lie close to a simple sine curve (Fig. 4). ‘The 
curve obtained by least squares is 


V = —21°51 + 5°31 sin (27t/o"- 177365), 
+ 0°22 +0°40 (3) 


t= J.D. —(2435250"-394 + 071773652), 

+ 04-002 (4) 
where FE is an integer (Fig. 4). ‘The r.m.s. residual from the least-squares solution 
is 1-9 km/sec, slightly larger than might be expected from observational inaccuracy ; 
if the discrepancies are genuine, they could represent either an interfering cycle 
of a variation in amplitude, but they are clearly not large enough to deserve 
detailed analysis. 


| 
20 4 t - fe 


+ 
| 


Fic. 3.--Observed radial velocities of Lup. 
Abscissae: Heliocentric J.D. — 2435000. 
Ordinates: Radial velocity. 
Vertical lines represent the observations, with the standard error estimated for each plate; the 


dotted curve represents the least-squares solution, assuming the period 04 -177365. 
2 
2 


-20 


ctric observations of 7 into one cycl 
Phase, in da 
Ordinates; ipper fw 
Lower 
T/ ivcles represent th ) tions; the full line in the up, figure represents the least-square 
lution Th phase lation hip heteween the Wil ¢ 4 clearly visible 


t' Lup was observed photoelectrically in two colours at the Cape Observatory 
on 1955 May 13 and 14, and the means of the results obtained with blue and 
yellow filters are illustrated in Fig. 4. Am is the mean departure from m, = 4"+127 
and m, = 4"-637. ‘Lhe uncertainty in the zero-point is about o™-or. 

he photoelectric observations indicate a range of o™-03 in Am and a period 
in perfect agreement with that derived from the radial-velocity data; furthermore 
the phase of light maximum coincides with that of the descent of the radial- 
velocity curve through the y value (Fig. 4), in accordance with the typical 
behaviour of 8 Cephei variables. +' Lup has thus been shown to be a 8 Cephei 
star of small amplitude, similar in behaviour to y Peg, Cet and €' CMa. The 
scatter in individual observations is too large for the detection of any change 
in colour during the cycle. A comparison of the results of 1955 May 14 
(= J.D. 2435242) with some Fabry photometric observations carried out during 
the past few years suggests that the range in light may vary by + o0™-o2 (A. W. J. 


Cousins, private communication) 


Che velocity amplitude of the star is 40 per cent below the value corresponding 
to the linear period-amplitude relation found by McNamara (1955). 

A visual examination of the spectrum of r! Lup shows a slight intensification 
of the O11 doublet A 4415 at light maximum (suggesting slightly earlier spectral 
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type) during certain cycles. Such intensification was found on April 6 
(=J.D. 2435204), April 29 (=J.D. 2435227), May 11 (=J.D. 2435239), and 
July 8 (=J.D. 2435297), but not on the other nights when spectra were taken. 
No changes were observed in the shapes of line profiles. 

(iv) « Lup. ‘The measured radial velocities of « Lup are illustrated in Figs. 5, 6. 
The observations show a well-defined periodicity with P =o0*-259882 + 20, and 
when combined into a single cycle, line close to a simple sine curve (Fig. 6). 
The curve obtained by least squares is 


2 E), 
(6) 


where £ is an integer (Fig. 6). ‘The r.m.s. residual from the least-squares solution 


ed radial 1 elon if 
Absa ! Heliocentric F.D 
Ordinates Radial velocits 


Vertical lines represent the observations, with the standard error « 
dotted curve represents the least-squares solution, assuming th 


V = + 3°98 + 7°33 sin (27t/o"-259882), 
+ 0°22 + 0°37 (5) 
t=J.D. —(2435297"-409 + 0425988 
+ 01-002 
| 
243500 
timated for each pl ile; the 
period 
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Phase 


bic. 6.—-Radial-velocity observations of « Lup, combined into one cycle, 


Abscissae: Phase, in days. 


Ordinates: Radial velocity. 


The open circles represent the observations; the full line represents the least-squares solution. 


is 1°83 km/sec, closely similar to that obtained for 7! Lup. This coincidence 
suggests that the discrepancies are probably due to observational inaccuracy in 
both cases. 

z\L.up was observed photoelectrically in yellow light on 1955 July 20 by 
A. B. Muller at the Leiden Observatory Southern Station, using the neighbouring 
fifth-magnitude star HD 129858 for comparison, « Lup being measured through 
a 2-5 neutral filter. ‘The observations were continued for 34 hours, during 
which time «Lup was not found to vary by more than +0™-o1; even this 
variation behaves in such a random fashion that a systematic variation of more 
than +0™-005 can be ruled out. ‘The apparent constancy in light of « Lup is 
confirmed by reference to some earlier observations of this star carried out at 
the Cape Observatory in 1948 and 1950 (A. W. J. Cousins, private communi- 
cation). A. B. Muller’s observations, averaged in groups of five successive 
measurements, are given in ‘lable V. 

It is conceivable that «Lup has variable light amplitude, like €'CMa 
(McNamara, 1955), so that some variation might be detectable on a more 
favourable occasion; however, the present data indicate constancy in light to 
a remarkably high accuracy and would be perfectly consistent with the hypothesis 
that « Lup is a spectroscopic binary. 

4. Absolute magnitudes of BCru, +! Lup and «Lup.—The purpose of the 
present section is to determine the absolute magnitudes of the three stars that 
have been found to have a short-period variation in radial velocity and to compare 
the results with the two period-luminosity relations that have been obtained 
for 8 Cephei variables. Owing to the very small range in intrinsic colour of the 
group, it is not possible at present to determine the intrinsic colours of the three 
stars on the U, B, V system sufficiently accurately for a detailed comparison 
with the results of McNamara and Williams (1955), apart from remarking that 
fh Cru must be intrinsically much bluer than any well-established 8 Cephei star. 
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TABLE V 


Leiden 


Photoelectric measurements of « Lup 


Heliocentric Am Hehocentric \m 
J.D 2435000 (vellow filter) J.D 2435000 (vellow filter 


31395 


1955 July 20 1955 July 20 (continued) 


412 
415 


6 Cru and « Lup belong to the Scorpio—Centaurus stream, and their absolute 
magnitudes have been determined by Blaauw (1946). However, these magnitudes 


may require some correction for interstellar absorption; the correction may be 


estimated in an indirect manner using the equation 


Am = 7(E, + 0"-05) (7) 


where £, is the colour excess on the scale of Stebbins, Huffer and Whitford ( 1940). 
Equation (7) is based on the work of Oort (1938) and Oosterhoft (1951). A com- 
parison of the colour excesses derived by Stebbins and his collaborators with 


those of Johnson and Morgan (1953) suggests the approximate relation 


t+o™-os=1E (8) 


“ 


where E, is the colour excess in (B — V) of an early B-type star. From the colour 


index C,,, given in Cape Mimeogram 1, we can now find (B V) using the relation 


C,,, = 0°964(B — V)—o™-188 (9). 


(Stoy, 1953) and compare it with the normal value for its spectral class given by 


Johnson and Morgan (1953). Equation (g) was tested for seven early B stars 


common to both lists and found to give a maximum error of o™-or. 


t' Lup does not belong to the Scorpio—Centaurus stream, and no parallax is 


available. Its absolute visual magnitude has therefore been estimated from the 


equivalent width of Hy, following Petrie (19534). For this purpose, the Hy line 


in suitable spectra of +! Lup and of a number of other stars of known absolute 


magnitude (Blaauw, 1946; corrected for absorption) was traced with the 


microphotometer at the Radcliffe Observatory (using 24 calibration spots from 


the tube sensitometer on each plate) and the magnitudes calculated from Petrie’s 


formula. ‘lhe results are given in ‘lable VI. 


Unfortunately, the stars studied are so intrinsically luminous that the 


magnitudes determined from Hy are very sensitive to errors in the equivalent 
width; the results do show, however, that the absolute magnitude of 7! Lup is 
not likely to differ from that estimated from the Hy measurements by much more 
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than the standard error determined trom their intrinsic scatter. Furthermore, 
our value for the equivalent width of Hy in the spectrum of 7 Sco is in tolerable 
agreement with Petrie’s value (3-3A). We may therefore use the results of 
Table VI to compare the observational data for 8 Cru, 7! Lup and « Lup with 
the period-luminosity relations of Blaauw and Savedoff (1953) and of Petrie 
(1954a). ‘This is done in Fig. 7; it is seen that r' Lup fits either of the two 
curves very well, « Lup comes tolerably close to that of Petrie, and 8 Cru departs 


from both curves by more than a magnitude. 


7 Periods and luminosities of B Cru, r' 1 nd a Lup, 
Period, in day 
Ordinates lhsolute visual 
The wertu lines represent the observational results for the thy fas 1 mited limits of error, 
Broken line: Period-luminosity cur of Blaa 1 Loff (1953) 
Full line: Period-luminosity curve of Pet I 


‘These results suggest the following conclusions. 
(i) 7’ Lup appears to be a typical small-amplitude 8 Cephei star. 

(ii) In view of the failure of all efforts hitherto to establish a variation in light 
of «Lup, it is not possible to conclude from the velocity variation that it is 
definitely an intrinsic variable rather than a single-lined spectroscopic binary, 
However, if we assume the mass of the primary component to be of the order 
of 10 solar masses, it is easily shown that a binary system completing an orbit 
in 04-26 would have positive total energy if the mean mutual distance of the 
centres of the components were greater than 3-3 x 10!!cm, or about 5 solar radii. 
As the radius of a B2 giant is probably of the order of to solar radii (cf. the 
discussion of BCep by Stebbins and Kron, 1954), it is clear that the binary 
hypothesis is not a very likely one. Clearly further studies of « Lup are needed, 
If it should turn out to be a 8 Cephei star, it would be of great significance in 
fixing the period—luminosity relation because its period is the longest known and 
its parallax relatively well determined. 

(iii) As has been mentioned, it is possible that the period of B Cru may 
require revision in the light of future observations; alternatively, the anomalies 
that have been found here may be connected with its early spectral type, which 
suggests that 8 Cru cannot in any case be a “typical” BCephei star. Further 
study of this star involving simultaneous observations of light and radial velocity, 
would be of the greatest interest (see note overleaf). 
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C‘ambridec 


1956 January 2. 


Note added in proof.—Since the completion of this paper, a measurement 
of the radial velocity of « Lup, carried out by W. Buscombe at the Mount 
Stromlo Observatory, Canberra, has been brought to my notice. Buscombe’s 
result, for 1954 March 24.701 U.T. (heliocentric J.D. 2434826.204), 1s 
(—2+1°2)km/sec, in perfect agreement with equation (5). I am grateful to 
Dr Buscombe for his permission to quote this result. 
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THE EVOLUTION OF UNMIXED STARS 


R. 7. Tayler 


(Received 1955 November 10) 


Summary 


The earlier treatment of this problern by Roy is shown to suffer from 


difficulties with regard to the transfer of energy in the convective core 


These can be eliminated if Roy’s assumed discontinuity of composition is 


replaced by a continuous gradient of molecular weight. Six such models 


have been obtained with central hydrogen content changing from 99 to 2 
9 


per cent. The evolutionary tracks of these models are compared with those 


obtained by Roy’s method and there is seen to be no significant change in the 


observable properties of the stars. ‘The H. R. Diagram for stars of varying 


masses but of the same age gives qualitative agreement with the lower part of 
the observed H 


R. Diagram for globular clusters 


1. Introduction.— This paper may be regarded as a sequel to one on the 
evolution of massive stars. [‘layler (18). ‘This will be referred to as Paper | 
in what follows.] It is also a reconsideration of a problem previously studied 
by Roy (14). Paper I contains a list of references to previous papers on stellar 
evolution; papers that have appeared since include Harm and Schwarzschild 
(7, 8), Hoyle and Schwarzschild (g) and Reiz (13). 


If a star evolves without mixing currents, and either it possesses no convective 


core or the size of the convective core decreases as it evolves, a region with a 


continuous gradient of molecular weight will be set up in the star. Sequences 
of models obtained by Schonberg and Chandrasekhar (16), Roy (14), and 
‘Tayler (17) have assumed a discontinuity of composition at the surface of the 


core, but all these series of models have decreasing convective cores so that 


zones of variable composition must be present. In Paper | models have been 


obtained for massive stars with zones of variable composition, but for electron 
scattering opacity (which predominates in such stars) no models with a dis- 


continuous change in composition at the surface of the core can exist. Hence 
no estimate could be made of the effect of the neglect of a transition zone. In 
this paper models are obtained similar to those of (14) but with a zone of variable 
composition replacing the discontinuity. At the same time a difficulty in the 
transport of energy, which occurs in the cores of some of the models of (14), 


is removed. ‘lhis will be discussed in the next section 


2. Outline of problem.—In (14) Roy considered the evolution of unmixed 
stars with a photoelectric opacity law and with convective cores. He assumed 


that as the stars evolved, with conversion of hydrogen into helium in the energy- 


producing regions, a discontinuity of composition was built up at the surface 


of the convective core. The outer envelope retained its initial composition while 
the helium content of the core gradually increased. As was pointed out in Paper I, 
the assumption of a discontinuity of composition is essentially approximate, 


3 
i 
{ 
i 


26 R. F. Tayler Vol. 116 


when the convective core mass decreases as the star evolves. In this case 
convection dies out in the outer regions of the core and, as this material has 
already converted some of its hydrogen into helium, a radiative zone of variable 
composition is formed. Since the change in size of the convective core in the 
models of (14) is very small, one might expect that the replacement of the zone 
of variable composition by a discontinuity would have little effect on the models. 
gee ‘There is however another difficulty associated with the transfer of energy in the 
core. 


At the surface of the convective core all the energy flowmg out from the 
centre of the star must be carried by radiation. One of the conditions used at 
the surtace of discontinuity of models in (14) is that the flow of radiation should 
be continuous across the interface and equal to the surface luminosity. Because 


of this the energy carried by radiation at any point of the convective core must 
be less than that at the surface of the core; otherwise convection currents would 
be required to transport energy against the temperature gradient. ‘lhe con- 
vective core is the Lane~Emden polytrope (m= 3/2) which is tabulated in (4) and 
for any given opacity law we can plot the radiative energy flux against Emden 
radius (or the Bondis’ variable S used in Paper I (14) and the present paper). 
‘The transter of radiant energy in convective cores has previously been considered 
by Opik (12) and the results shown here are essentially those given by him. 
Fig. 1 shows the results for Electron and Kramers opacity laws. We see from 


3.0 1.4 0.6 0.3 

1-0} 
Lp | 

Lu | 

/ 
0.5 ECTRON 


0.0 f = 
x 
bic. t.-Energy transfer by radiation in convective cores The curves show the energy 
arrived by radiation at points in the convective core of a star for two different opacity laws. For 
Bes : cach curve the flux is expressed in terms of its maximum value and 1s plotted against the Emden 


radius x and the homology variable S. 


this figure that there is a maximum value.of the Emden radius at which the 
boundary of the core can occur; for the two opacity laws we have 


Electron x <1°893, S>o-6, 


Kramers S>1°4. 


his condition in a more refined form has been used by Naur and Osterbrock (11) 


> 
5 
; 
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to obtain a necessary condition for the existence of convective cores in stars. 
They require that the energy flow by radiation inside the core must be less than 
the flow across the surtace of the core by an amount at least equal to that produced 
in the intervening region. We shall assume that we are dealing with sufficiently 
concentrated energy production for this correction to our criterion to be 
unimportant. 

An equivalent way of looking at this problem is that of examining the radiative 
equations near the interface. At the interface the value of N (= 1/(m + 1)) jumps 
to the convective value o-4, but in order that convective conditions shall persist 
N must remain greater than or equal to 0-4 as we go further into the core. ‘The 
radiative equations for Kramers’ opacity give (14) 


V(S+2-—85N) 


dS S(i +2S— 40) 
When V=o-4 on the Emden solution, 1+2S—4O0 -o and 
dN o-4(S— 1: 
4 4) (2.2) 


dS S(1+2S—40)° 
Thus if S<1-4, dN/dS is negative just inside the core, and the temperature 
gradient gives stability against convection once again. ‘This condition that, 
at the boundary of a convective region, the radiative temperature gradient must 
remain unstable has been previously stated by Ledoux (10). In (14) we have 
Kramers’ opacity and after the first three models of the sequence the convective 
core boundary occurs with S 1-4, implying convective transport of energy 
against the temperature gradient. 

In (17) it was shown that this same difhculty occurs with models having a 
modified Kramers’ opacity law « ~ p'’’°/7*° but no attempt was made to remove 
the difficulty. In the present paper, by taking account of the zone of variable 
composition formed when the convective core shrinks in size, we have also 
obtained models in which there is no difficulty about the energy transport in 


the core. Kramers’ rather than modified Kramers’ opacity has been used to 
facilitate direct comparison with (14). 

3. Assumptions and equations. ‘\Vhe models here considered consist of a 
helium enriched core in convective equilibrium surrounded by a zone of variable 
composition in radiative equilibrium and an envelope of constant composition 


also in radiative equilibrium. ‘lhe following assumptions have been made. 


(1) nergy production is by the carbon cycle and is confined to the convective 
core, 
(2) xcept in the convective regions there 1s no mixing. 
(3) Only photoelectric opacity is important. 
(4) The initial composition is 9g per cent hydrogen and 1 per cent heavy 
elements 
(5) lhe central temperature remains constant as the star evolves. 
(6) Radiation pressure is negligible. 
(7) Degeneracy never becomes important. 
Homology variables. ‘Vhe solutions were obtained in homology variables 
introduced by Bondi and Bondi (2) and defined by 
47 Pr PdT 


GMp’ ~ TaP (3-1) 


G M 
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where the third variable replaces the Bondi variable N=1—(P/p)(dp/dP) to 


which it reduces when the composition is uniform. 

Subscripts. s=surface; c=centre; 1=inner interface between the con- 
vective core and the zone of variable composition; 2 = outer interface between 
intermediate zone and envelope. 


Equations. In a region of variable composition, the hydrostatic equation 
and the equation of mass conservation take the form 


Sd0 S+N-O 5S dX 
1+2S—40 2-99+5XdS (3-2) 
where 
In a radiative region with opacity 
K + X)p/T*, (3.4) 
S dA 2018 dX 
dS 1+2S—4O0 (3-5) 
While, in a convective region, 
y=5/3, =0'4. (3-6) 


Limit of convection. ‘Vhe condition for instability against convection will be 


(cf. (18)) 
dT dT 
( dr ( dr (3-7) 


A A 


‘This relation reduces to 


adial (3.5) 
At the boundary of the convective core 


N =0'4. (3-9) 
4. Method of solution.—'Vhe method of construction of the models is 
essentially the same as that used in Paper I. The equations for the envelope 
possess a one-parameter family of solutions. In (14) Roy has solved these in 
terms of the variable 
u In(M/M,) (4.1) 
and near the surface (u =o) 
S = ult — 3°238100 + 2°752406" + 0°82594 +u+...] 
O=6[1 — 1:0119050 — 0°352246? — 0'01172 + 0'84314u+...] (4.2) 
= 0°235294 + ulo'078431 — 0157008 — 0:0565 567] 
where 6= Au’, Roy has obtained envelope solutions for eight values of the 
parameter A and we use these solutions in the construction of our models. 
In the convective regions there is a single solution satisfying the central 
boundary conditions 
O=5/3 + O(1), 


MN =0'4. (4-3) 
This is the Lane-Emden solution for polytropic index n=1°5. ‘This is tabulated 
in (4), and Roy has tabulated it in terms of Bondi variables in (14). We have 
both used Roy’s table and retabulated the function for smaller intervals of the 
Bondi variables, at different stages in the problem. 


he 
; 
= 5 
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For one value of the parameter A, the envelope solution will fit smoothly 
to the core solution, and provide a model for a homogeneous star. Such a model 
has been obtained by several authors including Bondi and Bondi (2) and 
Gardiner (6). Roy extrapolated from his calculated solutions to obtain the 
homogeneous model which was the starting model for his evolutionary sequence 
(he obtained good agreement with the model of Gardiner) and we adopted the 
same procedure. 

Once the homogeneous model has been obtained (A = 0-234), the second 
model of the sequence was constructed as follows. An envelope with slightly 
higher A( =0-2415) was taken and was used to the point at which it contained 
the same mass as the radiative part of the first model; this point was reached 
before condition (3.9) for convective instability was satisfied. Beyond this 
point we are integrating through a region which was originally in the convective 
core and has therefore converted some of its hydrogen into helium, and thus 
we are in a zone of variable hydrogen content. As a first approximation to the 
variation in this zone we assume that the hydrogen content Y is linear with u, 
so that in the transition zone 


dX 

=const. 

dX dX 

dS 1+2S—40 du’ (4-4) 


‘The integration was carried out for a trial value of dX/du. Initially our solution 
lies above the Lane-Emden solution in the S—O plane. For some values of 
dX/du, the variable 4 reaches the value 0-4 while our solution still lies above 
the Emden solution; for other values the curve crosses the Emden solution 
before .¥ =o-4. For one value of dX/du, .4 =0-4 at the point at which the 
solution crosses the Emden solution. At this point all the conditions for a fit 
to a convective core are satisfied. Several trial values of dX/du must be used 
and the required result obtained by interpolation between them. 

The variation of composition in this transition zone defines the variation in 
the first part of the transition zone of the next model, and the remainder of the 
transition zone is constructed as above. In all, including the homogeneous 
starting model, six models have been obtained and the mathematical character- 
istics of the models are listed in ‘Table I. It may be noted that in each case the 
fit to a convective core occurs with S > 1-4 so that there is no maximum of radiative 


energy transfer inside the core. Fig. 2 shows the change in hydrogen content 


I 
Il 


TABLI 


IV V VI 


A 0'234 0O°2415 0'2513 0'2568 0°2635 0°2680 
Se 2°120) 1°7757 1°3655 1°1559 0°9212 0°7774 
O, o'8104 0°7895 0°7740 0°7507 0°7318 
Ns 0°4000 0°3546 0°3636 0°3515 0°3362 0°3259 
1°3269 1°5035 1°4538 1°4970 1°4523 
O; 0'7220 0°6130 0°'5960 o'6106 O°5954 
0'4000 0°4000 0°4000 0°4000 0'4000 
My 1°9548 2°0448 2°2073 2°5690 2°5790 
(dX 6°30 0554) 0503 


X, 0°9900 0°7802 0°5225 0°3738 0°1734 O'O175 


1 
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and the retreat of the convective core as the star evolves. It can be seen that as 
the first fifty per cent of the central hydrogen is burnt the transition zone is fairly 
adequately represented by a discontinuity; after that the size of the transition 
region increases rapidly. 


cent ymposition ¢ 
for all models of the sequence, d thi the | tal line at the 
unbered horizontal line lro é ent and ext 
nodel lhe hori fal lin ue ) hoon 


composition in the intermediate zone. 


5. Physical properties of models.—Once we have obtained our solutions in 
terms of homology variables, we require to pass to the solution in physical variables. 
If the total mass, M,, and the central temperature, 7), are given, we require 
expressions for L,, P., p,, r, and 7,. From the definitions of the homology 


variables 

Pa P, 
+Iin— 
Gp,M, 


l 


>, 4 
cls l 


In S, =In +In 


Vi 


1 
—2In 
where we have taken these relations at the boundary of the convective core. 


In this equation 


In(M,/M,)= 
In(P,/P.), In(p,/p,) are tabulated in the B.A. ‘Tables (4) [B and A in those tables] 


and 


In(r,/r,)= — | (O/S) du 


must be obtained by quadrature. Using 
P.=Ap,T.)p 


we obtain expressions for r,, P., p,- 


0 H | 
0.6 
x H4 | 
0 
| 
3 0 vi 
0. 0.1 0.2 0. 
) M/Ms 5 
| 
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L, is obtained from the ‘Mass Luminosity”’ relation which in this case 


takes the form 
4 485 = 3K00Z,(1 + he 1/2) 
17 A 167ac (4 Ly M, (5.5) 


The effective temperature can then be found from the relation 
L,=nacr,?T 4. (5.6) 
The total mass of hydrogen that has been converted into helium (mass burned) 


for any model of the sequence can be obtained by integrating the hydrogen 
content over the mass 


-M, 
M,= (o-99 —-X)dM. (5.7) 
“0 
The time taken in the evolution is found from 
‘Time = 0-007¢* [aM,/L. (5.8) 


Once values are taken for kyo, 7, and M,, the physical properties of the models 
are completely determined. However the application of this particular com- 
bination of opacity and law of energy generation to any model is uncertain (3), 
and in ‘lable II we have expressed the physical properties non-dimensionally 
showing changes from a zero point which is only fixed when ko, 7,, and M, are 
given. We could have fitted our models to some of the solar parameters but 
that might have given quite misleading values for other solar properties. 


Taare I] 


140 


*Qgoo 


u 


By calculating | (O/S)du for values of u other than u,, and by using 


J 
equations (5.1) and (5.2), we can obtain the variation of r, P and p with M fo: 
different models of the sequence. ‘Table I1] shows the results for Models II-\ 
Model | has been omitted because of the extrapolation involved, and Model V1 
because its very low central hydrogen content implies that some of our basi 
assumptions are invalid. In Fig. 3 we plot the actual movement of mass shell 
against mean composition of the models. We see that most of the increase in 
radius occurs for the outer 66 per cent of the mass as would be expected fron 
the small change in central density. 
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“ABLE III 
Model III 
"904 ‘929 
B43 "BBO 
‘734 ‘ ‘796 
‘595 650 


"509 
‘500 
"410 
‘286 

195 

132 
‘083 
‘047 
'024 

007 
‘002 
‘000 


‘000 
‘S70 
"799 
‘O45 


097 47° 
100 ‘612 ‘441 
116 0°352 
14! "25! ‘250 
180 ‘169 


215 ‘110 
‘249 ‘10% ‘069 
"042 
"023 
‘O10 
401 ze ‘003 
00 000 ‘000 


020000900 


We can also consider the relation between luminosity and mean molecular 
weight for the models of the sequence. It is found that all but Model VI quite 
closely follow a power law 


(5-9) 
‘This result can be compared with the mass—luminosity—-molecular weight relation 
for unmixed stars of the same opacity, and with results for both mixed and unmixed 
stars with electron opacity. ‘These results previously considered in (19) are 
shown in ‘Table IV. It will be seen that for both opacity laws the increase in 
brightness is more rapid for the unmixed stars, the effect being larger in the case 
of Kramers’ opacity. 

We are now in a position to consider the validity of the assumptions made 
in Section 3. As has been mentioned previously, it is uncertain whether the 
combination of Carbon Cycle energy production and photoelectric opacity 
considered in this paper does occur in any stars. Bondi (3) considers that Carbon 


P'P, 
0°525 
O'142 9°150 0'667 
0°20 0'175 0506 O'153  0°353 0°334 
0°30 O'212 0°397 0°266 0°213 
as 0'40 0°246 0°307 ( 0°226 0°200 O°142 
0°60 0°316 0°163 0°297 0°056 
0°95 0°555 o’oll 0°540 0'006 o'0014 
0°00 0°000 1000 0°000 1°900 1*000 
0°O42 0'g20 0'034 0*goo 0°833 
0°05 0'074 0°769 ( 0°5607 
0'071 0-619 0°449 
o°146 O'101 0'209 0251 
201 o85 0'068 
: 236 060 0'042 
74 049 0°025 
316 025 0°'013 
369 O14 0'006 
150 006 00016 
522 002 00005 
100 000 #000 
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Cycle production predominates in stars somewhat less massive than the Sun and 
in these stars photoelectric opacity would be most important. On the other hand 
Epstein (5) believes that the P—P reaction predominates in the Sun. Bound 
up with this uncertainty is that of the existence and size of convective cores. 
For these reasons we shall not attempt to justify our choice of laws of energy 
generation and opacity ; we have been interested in the possible effects of replacing 
a discontinuity by a transition zone, and we have chosen conditions agreeing 


ooL 


0.5 1.5 2.0 


Fic. 3.—Motion of mass shells during the evolution. The curves show the motion during the 
evolution of surfaces containing given fractions of the star's mass. Radial distance, in terms of the 
radius of Model II, is plotted against fractional mass of hydrogen converted into helium, The 
numbers at the top show the percentage of the mass contained within the surfaces considered, 


with those used in (14), and avoiding the mathematical difficulties of mixed 
opacities and energy generation. In Paper I, where the opacity and energy 
generation were fairly clear, we were unable to compare our models with those 
having a discontinuity at the surface of the core, as such models do not exist 
for electron scattering opacity. 


TaBLe IV 


Le 


Model Opacity y é 
Mixed Electron scattering 3'0 4°78 
Unmixed Electron scattering 3°0 5°4 
Mixed Kramers’ 7°78 
Unmixed Kramers’ 


As in Paper I we can integrate a law of energy generation through the 
convective core, and, by comparing the energy generation at the core surface 
and just beyond with the mean generation in the core, we can estimate the 
importance of energy production outside the core. In Model VI we find a 
rapid increase of energy production as we enter the zone of increasing hydrogen 
content just outside the core, and in this model shell burning of hydrogen outside 
the core is very important. In both Models IV and V the energy production 


3 


M 
Mz | 


34 R. J. Tayler Vol. 116 


is small at the core surface, and continues to drop in the transition zone ; energy 
production outside the core is possibly important in Model V but not in Model IV. 
In Models II and III there is again an increase of energy production as the 
transition zone is entered, because the change in composition is almost dis- 
continuous, but this is soon balanced by the decrease further out. 

Finally, with the knowledge that but for Model VI most of the energy is 
produced in the core, we can estimate the increase of central temperature which 
we have neglected previously. Assuming that all the energy of Model V is 
produced in the core, we find that the central temperature will have increased 
by about one-fifth. ‘This is of the same order as the increase found in (14). 
It is likely to be an overestimate as some energy is produced outside the core. 
Such a change does not cause much alteration in the evolutionary tracks discussed 
in the following section, especially when one considers the uncertainties discussed 
above. 

6. Conclusions and comparisons.—Using the results shown in ‘lable II the 
positions of the sequence (relative to an arbitrary starting point) can be plotted 
in the H.R. diagram. ‘The results can be compared with those of (14) and at 
first sight it appears that the replacement of a discontinuity of composition by 
a zone of variable composition has made a clear difference to the models. 
However there appears to be an error in the results of (14). Equation (3.10) 
in (14) (equivalent to our equation (5.1)), 


Py, 


u O 
Gp,M, 


0 


N 


“ S 


=InO-u4 | du + | du, (6.1) 
must be handled carefully if the range of integration of the second integral 
contains a point of discontinuity of p. At such a point N (=1— P/p(dp/dP)) 
is momentarily infinite and this leads to a finite contribution to the integral. 

In fact 

where now N can be considered finite everywhere. 

We have taken Roy’s solutions in homology variables and have recalculated 
their physical properties by using equations (5.1) to (5.6) of the last section. 
‘This involves the assumption of constant central temperature so that the results 
are not exactly equivalent to those found by Roy. ‘The results of these calculations 
are shown in ‘lable V._ It can be seen that the introduction of the factor In (j4/p,) 
has virtually removed the apparent differences between the models of (14) and 
the present paper. Fig. 4 shows the tracks of the two sets of models in the 
H.R. diagram, and it can be seen that they only differ for the last two, and most 
uncertain, models of the sequence. ‘Thus it seems that, although the assumptions 
of (14) neglect the zone of variable composition 1d imply the convection of 
energy against the temperature gradient, they have ve. , ‘ittle effect on the properties 
of the star as a whole during the early stages of its evolution. Added to that the 
assumption of a discontinuity leads to much simpler mathematics. 

Also shown in Fig. 4 is the evolutionary track for fully mixed stars. Corre- 
sponding models in the mixed and unmixed sequences have the same mean 
molecular weight. Finally we show the track followed by the models of Sandage 
and Schwarzschild (15); in these evolution to a Schénberg—Chandrasekhar 
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limiting isothermal core (16), is followed by gravitational contraction. ‘These 
models also neglect the possible importance of a transition zone. 


TABLE V 
I Il III IV V VI 

al 1°00 1°42 1°65 2°33 

©7000 0°008 0°029 0°049 O'105 
Pe! Pel 1°00 0°96 0°37 0°73 
Moore/M, 0°146 O°142 0'137 0°136 0'135 
O'170 O'150 O'124 O'110 O'094 
X core ©°9900 0°7777 0°5138 0°3790 0°2300 0°1459 
Mass burned/M, 0:000 0°030 0'067 0'084 0°103 


T 


+0.05 0.00 -0.05 -0.10 015 
$LOG ole 
Fic. 4.—The evolutionary tracks in the H.R. Diagram. The straight line to the left shows 
the evolutionary track of a well mixed star. The upper curve is the H.R. Diagram for the models 
of Sandage and Schwarzschild. The lower curve is the H.R. Diagram for the sequence of models 
of the present paper; the points Vp and VIR show the deviations of models with a discontimuty 
in composition from the solid curve. 


Fig. 5 shows the instantaneous H.R. diagram for stars of different masses 
which started their evolution at the same time. Assuming that the central 
temperature is effectively constant for different masses, we can apply the 
homologies 


LoM®, rocM, Timex<(M/L). 

In Fig. 5 the most massive star has completed the evolution from Model I to 
Model VI, while the lightest has, in the same time, scarcely started its evolution. 

It will be seen that Fig. 5 is very similar to the one obtained for massive stars 
in Paper I and to that obtained by Sandage and Schwarzschild. In this case 
the “‘ knee” in the diagram occurs at the point of 8 per cent burning of hydrogen. 
This is rather lower than the results obtained previously. ‘The curve is in 
qualitative agreement with the lower part of the observational curve for globular 
clusters ((1) etc.); but it seems that almost any series of models that assumes 
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a low enough degree of mixing in stellar interiors will obtain such qualitative 
agreement. At the moment physics rather than observational astronomy must 
be used to differentiate between the theoretical results. 


0.2 0.1 0.0 
SLOG ole 


Fic. 5.—Instantaneous H.R. Diagram for the stars of varying masses. The straight line at 
the left of the diagram is the original main sequence; the curve at the right is the locus for stars of 
varying masses after a time depending on the unit of mass and other parameters of the problem. 
The connecting curves are evolutionary tracks of individual stars whose masses (in arbitrary units) 
are shown at their starting points. 


This problem would have presented far more difficulties had Dr Roy not 
deposited detailed tables of the models of (14) with the Royal Astronomical 
Society. I am grateful to Dr Roy, and to the Royal Astronomical Society, for 
the loan of these tables. 

Note added in proof—Since this paper was submitted, Roy has published a 
further note* in which the error in his equation [14, (3.10)] has been corrected 
‘by the introduction of the factor in y/u,, and his evolutionary tracks then 
essentially agree with those found above. 

A.E.R.E., 
Harwell, Didcot, 
Berkshire : 
1955 November 9. 
* M.N., 118, 396, 1955. 
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Summary 


A series of line-of-sight velocity measurements, made on spectra of the Sun 
obtained at Oxford in 1953 May, is described. Spectra some 15 mm in 
height, taken with a solar image of 85 mm diameter, covered an arc of length 
4 “10°km in the surface extending inwards from the limb. Measurement 
was made of two plates, taken an hour apart at the East equatorial limb, each 
carrying four exposures obtained in rapid succession. Each of the eight 
spectra was measured at 140 points o‘1O mm apart, corresponding to an 
average spatial resolution of 4 » 10* km on the solar surface. ‘The error to be 
attributed to a velocity measurement is carefully considered ; with the fine- 
grain plates used in the present investigation, a value + 0°09(0)+ 0:02(8) 
km sec”! is found for the error of measurement of a single velocity. 

An investigation of the equatorial velocity field is based on the 1120 
measurements of velocity. Measurements from the four exposures on a plate 
were combined to give the best values of the velocity, thereby reducing the 
uncertainty of the velocity at any point on the surface to + 0°04(5) + 0°01(4) 
km sec”'. ‘The residual velocities remaining after the limb effect and the 
component of solar rotation have been removed show a dispersion of + 0°16(1) 
km sec™'. This dispersion exceeds the error of the velocity at a point on the 
surface by some eight times the uncertainty of the latter, thus providing 
conclusive proof of the existence of a velocity field. An irregular distri- 
bution of the velocity residuals over the solar surface is found ; correlo- 
graphic analysis gives 2°6 * 104 km as the average linear extent of a velocity 
fluctuation. ‘The stability of the field is established from the similarity in the 
velocity distribution at the beginning and end of the period of observation. 
Finally, a suggestion is put forward for the possible origin of the field. 


Introduction.Measurements of the rotational velocity of the Sun, using 
Doppler displacements of spectral lines, have always shown much larger variations 
from one determination to another than could be explained by the random error 
of measurement (2, 9, 16, 22, 24). Whether these discrepancies, amounting to 
a few tenths of a kilometre per second, were the result of a variation of the 
photospheric circulation with time or were due to local variations in velocity 
was not, however, clear. 

In a previous paper by the writer (11), measurements of velocity, made at 
intervals along arcs extending inwards up to 2°5 x 10° km from the limb in 
the equatorial region, were described. A dispersion in velocity of + 0-17 kmsec™! 
was found, which significantly exceeded the error of measurement, estimated as 
+ 0°08 + 0-02 km sec! fora single velocity. Since it was also shown that scattering 
of light in the terrestrial atmosphere and by the telescope (14) could not produce 
spurious velocity differences of this size, as had once been suggested (7, 8), clear 
evidence was obtained for the existence of motions in the photosphere. 
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The linear scale of the velocity fluctuations was fairly large, varying from 
about 2°5 to 8-5 x 10'km, and the observations indicated that the motion was 
quite stable, persisting perhaps for 3 hours. As a precise determination of the 
velocity field the investigation was, however, inconclusive for several reasons. 
‘The true velocity variation along an arc of surface was partially masked by the 
random velocities due to errors of measurement. Furthermore, even the value 
+o-o8kmsec"! was only achieved by measuring over such a length of solar 
surface that features extending less than 2 x 10'km could not be resolved. Finally, 
although some velocity sequences covered a length of 2:5 x 10°km, most were 
restricted to about 1-3 x 10°km near the limb, so that the general distribution of 
velocity, even in the equatorial region, was hardly established. 

Attention in the present paper has, therefore, been focused on a detailed 
determination of the velocity variation along an arc of length 4 x 10° km, in the 
equatorial region, over a period of an hour. Satisfactory spatial resolution has 
been achieved by measuring velocities at 140 positions along the arc, corre- 
sponding to an average interval of 4 x 10% km between two successive positions. 
Spectra of moderate height covering the necessary length of arc in the solar 
surface have been obtained by reducing the size of the solar image. Errors in 
the measured velocities have been minimized in two ways. ‘The actual accuracy 
of measurement was increased by using finer-grain plates. Secondly, both at 
the beginning and at the end of the period, repeated observations were made 
of the velocity field. Altogether the velocity at 1120 positions has been obtained. 

In the sections which follow, an account is given of the observations made 
at Oxford in 1953 May, the measurement of the spectra obtained, and their 
reduction to velocities on the Sun. An analysis of the results is carried out in 
the last section and their physical implications are briefly considered. 

Papers in this series ‘“‘ Motions in the Sun at the Photospheric Level ”’ are 
referred to by their serial numbers. 

1. Observations.—In its normal form the Oxford solar equipment consists of 
a vertical telescope of the Cassegrain type (17) with an equivalent focal length 
of 198m, feeding a fixed horizontal Littrow spectroscope (18) whose slit lies 
in an E~W direction. ‘lhe diameter of the image produced by the telescope 
is 180mm, and in the wave-length region 5900 A the dispersion of the spectro- 
scope is 1:5Amm''. A guiding plate (19) placed in the beam in front of the 
spectroscope permits accurate focusing, and a precise determination of the 
position of the solar image on the spectroscope slit. 

For the present set of observations, the size of the solar image was reduced 
by modifying the telescope. ‘The hyperboloid and its carriage were removed 
and were replaced by a lens carried in vertical ways. ‘lhe lens, supplied by 
Sir Howard Grubb, Parsons and Co., is a visually corrected achromatic doublet, 
10 cm in diameter and of focal length g-1 m. A parallel beam of light was fed 
to it in the usual way by the coelostat and secondary mirror. ‘The converging 
beam from the lens was then deflected in a horizontal direction by the coudé 
mirror, the main paraboloidal mirror of the telescope, just below, being kept 
covered. An image of the Sun 85 mm in diameter was formed at the spectroscope 
slit. ‘The telescope was thus used as a simple achromatic refractor. 

‘The centre of the lens cell was moved laterally until a plumb-line from it 
intercepted the centre of the coudé mirror. By means of squaring-on screws 
and a sensitive level the lens cell was then made perpendicular to this vertical 
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axis. By moving the lens up and down in its ways, and observing the definition 
of the image in the eyepiece of the guiding device, the position of the lens for 
which the image lay precisely in the plane of the spectroscope slit could be 
detected. ‘This position of best focus was sharply defined and remained unchanged 
throughout the fortnight when observations were made. 

The spectroscope was adjusted so that the wave-length 5900 A was central 
on the plate. A diaphragm over the spectroscope slit limited its height to 
approximately 19mm; about 15-5 mm of this was covered by the solar image, 
so guided that the slit intersected the limb close to the equator. While, say, 
the East limb was accurately maintained in position, four exposures were made 
on a plate. Each exposure lasted 10-20 seconds and all four were generally 
made within 34 minutes. ‘The procedure was repeated at the same limb position 
at intervals of approximately an hour. Similar plates were obtained at the 
West equatorial limb. Ilford thin film panchromatic half-tone plates were 
used; they were developed in Parkhurst two-solution developer. 

The first fortnight in 1953 May was a period of exceptionally good weather. 
During that time, six West limb plates and three East limb plates were obtained 
with seeing varying between o and 2 on a scale of 5 (19). ‘I'wo East limb plates, 
separated by an interval of one hour, were selected as most suitable for 
measurement. Details of the exposures are given in the first two sections of 


Table I. 


TaBLe 


Spectrum Date and U.T. Exp. time Sp. ht Slit width 
sec mm mm 


1953 May 5° 


1953 May 


Centre 1953 May 21°4558 
standard 


1mm on spectrum normal to the dispersion = 22°°49 on Sun. 
1mm on slit 22°°33 on Sun. 

When all the limb spectra had been obtained, the telescope was restored to 
its usual form. With the solar image focused 50 cm in front of the slit, so that 
local irregularities of velocity are smoothed out on the spectra, several exposures 
were made of the centre of the disk. ‘The last row in Table I gives details of 
the spectrum selected as a standard for the velocity measures. 


2. Measurements 

2.1. General method of measurement.—Reliable values of velocity from Doppler 
displacements of solar lines are obtained when terrestrial atmospheric lines are 
used as standards of zero velocity (16,19). Six unblended solar (Fe) lines, 
bracketed by six water-vapour lines of comparable strength, were measured in 


362 | 15°42 0°04 
| Ev §°3935 15°44 0°04 
| 5°3943 15°41 0°04 
| 5°3649 15°44 0°04 | 
| Ew 54020 10°5 15°43 | 
§°4027 10°5 15°43 0°04 
54033 10°5 15°46 0°04 
5°4044 10°5 15°39 0°04 
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allthe spectra. ‘Table II gives the wave-lengths recommended by the International 
Astronomical Union (3, 10) and the Rowland intensities of the lines (23). 


TABLE II 


Terrestrial lines Solar lines 


Wave-length | Rowland Vel. factor 


Wave-length | Rowland 
Element 
A intensity A intensity km sec”! mm 
5go1'468 | 6 5905 ‘6080 4 Fel 74°69 
5919°054 5 §916°257 3 | Fel | 75°18 
5919°644 7 5927°797 2 Fer | 75°74 
5932092 5 5929682 2 Fer | 75°82 
5946-006 3 5930°191 6 bel 75°55 
5947066 2 | 5934°665 5 | Fer 76°06 


All the measurements were made on a modified Hilger photo-measuring 
micrometer. ‘lhe plate is adjusted on the stage with the dispersion parallel to 
the run of the main measuring screw. A diaphragm in the eyepiece limits the 
height of spectrum visible and a single measuring wire, perpendicular to the 
dispersion, lies in the plane of the diaphragm. ‘The stage carrying the plate 
can be moved by a second screw in a direction perpendicular to the dispersion, 
permitting measurement at various heights in the spectrum. From a series of 
fine scratches parallel to the dispersion, previously made at both ends of the 
spectrum at intervals of exactly 1 mm, the position on the solar image at which 


measurement is made can be determined. 

Velocities corresponding to the measured line positions in a spectrum are 
accurately found, independently of asymmetry in the apparatus function of the 
spectroscope, by reference to a standard spectrum. Measurement of the centre 
spectrum (‘lable 1) was made in the forward and reverse directions at 19 heights 
Imm apart, using a diaphragm height of 0-32 mm. ‘lhe standard spectrum 
was then taken as the mean of the 38 readings for each of the 12 lines. It gives 
a velocity (calculated from the Hartmann dispersion formula) of + 0°17(5) km sec"! 
and corresponds to a point situated o-o5 mm to the East of the centre of the 
disk. Averaging measurements from different heights in the spectrum largely 
eliminates the effects of irregular silver grain distribution, as the low value, 
+ 0°03(3) km sec™', for the root mean square error of this velocity shows. 

2.2. Preliminary measurements.—Spatial resolution on the solar surface 
decreases but the consistency of measurement increases with the height of 
spectrum over which measurement is made. With a solar image 180 mm in 
diameter, a diaphragm height 0-75 mm and coarse-grain plates, an average 
value +0°25kmsec™! for the root mean square error of a single solar line was 
obtained in the previous investigation ((11)——Paper IV). 

In order to find the best spectrum height in the present conditions, diameter 
of the solar image 85 mm and fine-grain plates, a short series of preliminary 
measurements was made. Each spectrum on the selected pair of plates (‘Table 1) 
and eight spectra on a similar pair of plates from the West limb were measured 
at 15 points, with a diaphragm height 0-32 mm (chosen to give about the same 


| 
| 
| 


42 A. B. Hart Vol. 116 


spatial resolution as the earlier measurements). For 60 of these velocities the 
value of the root mean square error of a single solar line (Paper IV, Section 4.1) 
was calculated, the average value being +0-14(0)kmsec!. Remembering that 
the diaphragm height in the present case is less than one-half that used in the 
earlier measurements, it is clear that a factor of at least two in the accuracy has 
been obtained by the use of the finer-grain plates. 

Attention in the final results is centred on the error attributed to the velocity 
at some point on the solar surface. As the velocity at any point is the mean 
from six lines and four independent exposures (‘Table I), the final errors should 
be smaller than the error of a single line-velocity by a factor of nearly \ 24. 
A diaphragm height of approximately o-1 mm was therefore decided upon for 
the final measurements. 

2.3. Measurements of limb spectra. Definitive measurements were made of 
the eight spectra on the two East limb plates listed in Table I. Over the one or 
two millimetres of spectrum height nearest the limb, darkening towards the 
limb caused a noticeable decrease in density. As the consistency of measurement 
depends to some extent on the apparent contrast between the line and its adjacent 
continuum, the first point of measurement was taken a distance 1-34 mm in 
from the limb. (The limb appeared quite sharp on the spectrum owing to 
moderate seeing conditions and the use of contrasty plates; its position was 
estimated under low magnification with a cross-wire.) Each spectrum was 
measured direct and reversed at 140 points over a height of 0-13 mm at intervals 
of o-10 mm, giving in all 1120 velocity measurements. 

lor any position on the solar image, let A, be the difference between corre- 
sponding water-vapour line (wave-length A,) readings in the limb and the 
standard spectrum and let A,, A, be the quantities for the solar lines. ‘Then the 
dispersion of the limb spectrum is converted to that of the standard by the 
A,~ A, curve; in practice a straight line is taken through the centre of gravity 
of the six points A,, A,. ‘The difference, 5,, between A, and the value of the 
adopted A, — \, curve at wave-length A, is the displacement of the limb line relative 
to the same line in the standard. ‘The factor converting linear displacement to 
line-of-sight velocity can be calculated from .the Hartmann dispersion and 
Doppler formulae; values corresponding to the dispersion of the standard are 
listed for the six solar lines in the last column of ‘Table I]. Multiplication of 6, 
by the appropriate factor from this table gives line-velocities, u,, relative to the 
standard spectrum. ‘The spectrum-velocity, L/, at the solar position considered, 
is the mean of the six line-velocities u,. 

‘The variation of line-of-sight velocity with position on the solar image is 
shown in Figs. 1 and 2. ‘The four spectra E,,, Ey, E,3, Ey, in Fig. 1 are those 
obtained at the beginning, E,,, Ex2, E.,, Ey, in Fig. 2, those obtained at the end 
of the interval of one hour (see Columns 1 and 2 in ‘Table 1). Let U,, be the 
velocity of the nth point of measurement occurring a distance h,, on the plate 
from the limb. a runs from 1 to 140 as the velocities were measured at 140 points 
on each spectrum. Abscissae are then the distances A, along the E-W chord 
of the solar image defined by the spectroscope slit. For the purpose of presen- 
tation the velocity data have been smoothed in order to reduce the effects of 
random errors of measurement. Values of 4(U, ,+U,,+ U,,,,) are therefore 
plotted in the figures as ordinates; the velocities are referred to an arbitrary 
origin, common to all the spectra. 
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2.4. Accuracy of measurement.—Errors due to measurement were investigated 
for a sample of 112 from the 1120 measured velocities. ‘The calculations were 
made at corresponding heights in each spectrum, for 14 positions distributed 
at regular intervals over the spectrum height. In this way, any systematic 
variation of consistency with photographic density is taken into account (see 
Section 2.3). It is clear that the accuracy of measurement depends on the 
deviations of both atmospheric and solar lines forming the velocity. ‘The 
discussion in the present section is outlined in Paper IV, Section 4.1, where the 
basic formulae are derived. 

If s, is the root mean square error of measurement of a single solar line- 


velocity u,, 


n 
(1) 
where n=6, since the spectrum-velocity is the mean of six line-velocities 
Let &, be the average of the 112 values of s,. ‘The line-velocities are regarded 
as a random sample from a hypothetical infinite population of all possible 
measures. An estimate, o,, of the root mean square error of a single line-velocity 
in this population is 


(2) 


the standard error of this estimate being approximately o,/[2(m—1)]". 
The deviations of the six water-vapour lines, at wave-lengths A,, from the 
A,—A, curve (Section 2.3) contribute to the error of measurement in a way 


analogous to the solar line deviations. A ‘ velocity-factor” can be attributed 


to the wave-lengths A, converting linear deviations to velocity deviations, v, say, 
comparable with the solar U—u,. ‘The root mean square error of measurement 
of a single atmospheric line is then 


= 6, (3) 


the mean of the u-values being zero. §&,, o, and o,/[2(n—1)]"* can be calculated 
as in the solar case. 

Fig. 3(i) and (ii) shows the distribution of the 112 values of s for the solar 
and atmospheric lines, together with the predicted theoretical distribution 
(Paper IV, equation (15)). Comparison of the two sets of values of 4, o and 
a/|2(n—1)|'* is made in Table III], columns 2, 3 and 4. Column 7 gives the 
displacement in microns on the plate to which the value of o corresponds. 

While dispersion of the solar lines about their mean appears as a random 
error of measurement, dispersion of the water-vapour lines about the A,— A, curve 
introduces a common error into the individual line-velocities and their resultant 
spectrum-velocity. ‘This systematic error, attributed to the centre of gravity of 
the six points A,, A,, is the root mean square deviation, o,/\/6, in A,. Hence 
the root mean square error of measurement of a single line-velocity is 


fa,2/[2a(n—1)] x 2(n—1)]}"?, 
and the root mean square error of measurement of a spectrum-velocity is 


+ + x 2(n—1)]+0,7/[6 x 2(n—1)])"?, 
a ai j 


I n 
n—1 


Vol. 116 


A. B. Hart 


ay] fe Suopy qual) ua Mf IDSSLISG 


‘LAwtt topy £861 yooda uvaw potaojonba ysogq avjos uo uotisog yjtor py ut 4- 


SHA NO GWITD JONVISIC 


| 
L510 


yBO 


LN 


VAN 


tus 


00 
20+ 


00 


cUee 
cv 
g 
VV a 
4 
20+ 
2 


- 
| 
~ 
* 
~ _lo 
: 
oO 
ro) 
| 
+ 


pS Avpy £861 uvawu ‘qust) potaoqonba avjos 24} uo uoytsog ut 1—t 


WY TOS NO 3HL JONTLSIO 
0S! Ovi 02 0 00! 06 


9S -wSE-y6O 


- 


= 


= 

= 

~ 

8 
~y 
§ 
= 


| 11-59! 


e 
No. 1, 1956 ee 45 
~ = 
aad w ws 
= 
~ 
— 
wy 
oo oo 6 Oo aa 


A. B. Hart 


(i) (ii) 
DISPERSION IN ! DISPERSION IN WATER - 
SOLAR LINE ERRORS. ' VAPOUR LINE ERRORS. 
Yo FREQUENCY. 
—-20- 
— 15 — 


10 


00 01 01 0.2 03 
ROOT MEAN SQUARE ERRORS, S 


Fic. 3.--Frequency histograms of the distribution of the root mean square errors of solar and 
water-vapour lines in a sample of 112 velocities. 
Abseissae: Root mean square error of measurement of a single line. 
Ordinates; Actual and percentage frequencies of occurrence cf a given error. 
Observed distribution, irregular curve; theoretical distribution, smooth curve. 


Values of o,,=0,/1/6 and the corresponding root mean square error of the 
solar spectrum-velocity are listed in column 5 of ‘Table II]; standard errors, 
o,/{2(n—1)]"*, of the values o,, are given in column 6. Extracting the relevant 
values from ‘Table II] gives as the error of a single line-velocity 

+ 0°15(7) + 0°05(0) km sec"! 
and for a spectrum-velocity 

+ 0°09(0) + 0°02(8) km sec™!. 


Taare [Il 


km sec km sec ' km sec '! km sec ! km sec ' | microns 


Solar + 0'13(4) +t o'14(1) + 0°04(5) + 0°05(8) + o'o1(8) 
Water-vapour | + 0°16(o0) + 0'16(8) + 0°05(3) + 0'06(9g) + 0°02(2) 


Comparing the value 4,= + 0°13(4)kmsec~! with the value +0-14(0)kmsec™! 
obtained in the preliminary measurements, it is reassuring to find that reduction 
in the height of spectrum measured has not greatly changed the accuracy of 
measurement. ‘The apparent gain in accuracy is illusory, the preliminary 
measurements having been made near the limb, where, as explained in Section 2.3, 
measurement error is highest. 

Calculations of errors of measurement made in Paper IV only took account 
of the dispersion of the solar lines about their mean, so the values given there 
should be somewhat’ increased. Rather similar considerations apply to the 
error quoted in Section 2.1 for the measurement of the standard spectrum. 
Here, however, dispersion of the water-vapour lines contributes to the error 
in a more complicated manner; the value already derived can be regarded as 
a lower limit to the actual error. 
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3. Solar velocities 

3.1. Reduction to velocities on the Sun.—The method of converting positions 
on the solar disk to positions on the Sun and measured velocities to line-of-sight 
velocities on the Sun is fully explained in Paper | (19). 

For the standard centre spectrum the calculated velocity due to Earth’s 
motion and solar rotation can be compared with the observed velocity found 
from the Hartmann formula (Section 2.1). A difference 0-07(1) km sec”! exists. 
As common vertical motions of the photospheric matter amounting to as much 
as 0-15 km sec~ may take place over several degrees of longitude (20), the observed 
velocity was adopted as the standard velocity V, to which all the limb velocities 
are referred. Possible error in the wave-lengths of the solar lines relative to the 
atmospheric ones (1, 20) and asymmetrical distortion of the line profiles by the 
horizontal apparatus function of the spectroscope (5, 20) may falsify this 
velocity, but in this case the amount is at least determinable. 

Heliographic latitudes, B, and longitudes, L, were computed for positions 
on the limb plates at which measurement had been made. Corrections, V,, 
were made to the measured velocities, V, for the effect of the observer’s motion 
at the moment of observation, giving 1120 values of V + V,+ V4, the velocity 
on the Sun. 

3.2. The velocity field—The four spectra on either plate (Table 1) were 

measured at exactly the same 140 disk positions (see Figs. 1 and 2). Furthermore, 
the interval of 3} minutes between the first and fourth exposure is so short that 
the four sets of measurements refer to the same positions on the solar surface. 
At each of the 140 positions a mean velocity was, therefore, calculated. From 
the dispersion of the four measures at any given point, a value can be found for 
the root mean square deviation of a single velocity, which can be compared with 
the root mean square error of measurement given in Section 2.4. ‘The deviations 
were calculated for the same 14 spectrum positions for which the measurement 
error had been calculated. The average of the 28 values for the two plates 
was +0°09(3)kmsec'. This value falls within the estimated uncertainty, 
+ 0°02(8)kmsec', of the measurement error +0-09(0)kmsec™! for a single 
velocity. Hence it may be assumed that no significant changes in velocity occurred 
during 34 minutes, and taking mean velocities for each plate is justified. ‘The 
values so found should represent, more accurately than those from individual 
spectra, the velocities actually occurring on the solar surface. 

Interest in the present work centres mainly on the variation in velocity from 
point to point, so that it is convenient to remove from the velocities the limb 
effect and the component due to solar rotation. Accordingly, the residuals 6V 
of the plate means V + V,+V, were calculated for the 280 solar positions from 
the equation 

bV=V 4 V, Ucosy,—A (4) 
(Paper IV, equation (1)). U, A are the velocity of solar rotation and the limb 
effect at the point considered, y, is the angle between the parallel of latitude at 
the point and the line-of-sight. Unfortunately, the variation of limb effect 
with position on the solar disk is poorly determined, interferometric measurements 
(1) disagreeing with values found from velocity measurements (11). However, 
in the present case no measurements were made very near the limb and the 
values found previously by the writer were used (Paper IV); at most the 
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correction amounted to o-16(0)kmsec”!. Since the zone of latitude covered 
by the observations fell within the range |B|<4°-o, the linear rotational 
velocity U was taken as U,cos B. For U,, the equatorial velocity, 1-91 km sec"! 
was adopted (Paper IV). ‘Table IV(i) contains the 140 velocity residuals so 
calculated from the mean of spectra E,,, E,,, Ey, E,, (Table 1) together with 
the heliographic latitudes and longitudes of the positions on the Sun where 
velocities were measured (Section 3.1). Similar information from spectra E,,, 
Foo, Egg, Egg is given in Table IV (ii). Columns headed £,, contain the values 
of 6V for E,,, ..., Ey, and Ey), ..., Ey, respectively. It should be emphasized 
that these values of the velocity residuals have not been smoothed. 

‘The mean velocity residuals, E, = — 0-00(8)kmsec~!and E, = +0°01(2)kmsec™, 
were calculated for each plate. ‘The small difference 0-02 km sec~! between them is 
very satisfactory in showing that no spurious velocity difference exists between 
measurements on the two plates. That the mean of E£, + E, is so closely zero is 
to some extent fortuitous, since it depends on the velocities adopted for U, and V, 
(cf. Section 3.1). From the point of view of having separated the velocity 
fluctuations from the systematic velocity, the position is entirely satisfactory. 

The variation of residual velocity with heliographic longitude (tabulated in 
‘Table IV (i) and (ii)) is plotted in the lower half of Fig. 4. A given difference 
in longitude covers a varying length of arc on the surface, depending on the 
latitudes of the points on the arc. The linear distance between two points of 
observation can be calculated from the B and L entries in ‘Table 1V. (The 
observations being in the equatorial region, the linear distance corresponding 
to 1° L varies only very slightly; on the average 1° L ~ 1-232 x 104 km.) 
Precise linear intervals corresponding to the angular longitude ones are shown 
in Fig. 4. 

Between the mean time of observation of the velocities E, and FE, a period of 
nearly an hour elapsed. Inspection of Fig. 4 shows that many features of the 
velocity field have persisted for this time. For example, the features between 
L= 236° and 240° are obvious in both FE, and FE, and the much smaller feature 
centred at 248°°8 is quite recognizable in both sequences. As entries in 
Table IV (i) and (ii) show that, at a common longitude, observations F, and EF, 
were made at latitudes differing by not more than 0°-2, a mean velocity curve 
was constructed. At any longitude L where a velocity, E, say, had been measured, 
the value of FE, was found by linear interpolation between two measured velocities ; 
the mean, E, of these gave a point L, FE on the mean curve. A similar procedure 
gave velocities at longitudes of observation of E,. Successive points on the mean 
curve were joined by straight lines; this is the curve E shown in the upper part 
of Fig. 4. It can be regarded as the quasi-steady velocity field for the period of 
observation. 


4. The solar velocity field 

4.1. Analysis of the velocity field.—The root mean square deviation of a 
single velocity, as calculated directly from the residual velocities given in 
Table IV, is +0°16(1) kmsec™'. Since each velocity FE, or E, in the 
table is the mean of four values, the root mean square error of measurement to be 
attributed to a single velocity E, or E, is + 0-09(0) + 0°02(8)/+/4 km sec (Sections 
2.4 and 3.2), i.e. a value +0-04(5)+0°01(4)kmsec~. The dispersion of the 
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IV (11) 
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velocities from point to point thus exceeds the error of measurement by some eight 
times its estimated uncertainty. Any suggestion that the velocity dispersion is 
a spurious effect due to errors of measurement can, therefore, be totally rejected ; 
the variations in velocity shown in Fig. 4 essentially represent the motion 
occurring in the Sun. Moreover, the observational procedure has so reduced 
the absolute value of the uncertainty of a single velocity ( + 0°04(5)kmsec~'), 
that the curves FE, and F, actually depict the motion with a considerable degree 
of precision. 

The linear scale of the velocity fluctuations is now required. ‘his was found 
from a correlogram analysis of the mean curve E (Section 3.2, Fig. 4) since it 
may be taken as the steady field occurring during the period of observation. 
Values of E to o-o1 kmsec™! were tabulated at each o'-1 interval of longitude 
over the 34° of longitude covered by the curve. Serial correlation coefficients r,, 
defined by 
n—k 


ly = (4) 


\J 

(13), were calculated for k= 1, 2, 3, ..., 30, 32, 34, 30, where is the order of 
the correlation coefficient and n= 340 is the total number of tabulated values. 
(‘his simple expression for the coefficient can be used when the mean residual, F, 
is approximately zero, cf. Section 3.2.) Values of r, are plotted in Fig. 5 against 
the longitude interval k/10 ', corresponding to the order of correlation k. ‘The 
correlogram is formed by joining successive points; ry=1 by definition. Apart 
from this trivial initial maximum, the correlation coefficients have a maximum 


at a longitude interval 2-1. ‘This then gives, using the mean conversion factor 
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1 1:232x10'km (Section 3.2), a distance 10'km as the average 
linear extent of a velocity fluctuation. No great significance attaches to the 
value, however, as the curve E shows clearly that the field is extremely irregular. 
The broadness of the maximum displayed by the correlogram and the very small 
value, +0-oo1, of the correlation coefficient at maximum reflect this absence of 
a single well-defined size of the motion element. 

Questions concerning the duration of the field remain to be discussed. ‘The 
general similarity of curves E£, and E, in Fig. 4 is quite marked; one wishes to 
know whether it is a chance resemblance or whether it denotes the persistence 
of the field over the hour’s interval between the observations. ‘The correlation 
coefficient (27a) between the two curves E, and E, (as calculated from the 208 
pairs of points from which the mean curve E is constructed ; cf. Section 3.2) 
is +0°607. Application of Fisher’s z-¢ test (27b) shows that the probability of 
this correlation arising by chance from an uncorrelated population ts less than 
1 in 1072, Clear evidence of the stability of the velocity ficld is thus provided, 

Some comparison can now be made of the present results with those obtained 
in the previous investigation. ‘The picture of an irregular field, with velocity 
fluctuations of amplitudes upto + kmsec extending some tens of thousands of 
kilometres, is very similar in the two cases. A precise comparison between the 
spatial distribution of the motion is difficult to make because of the great 
ditference in spatial resolution; the average value 2-6 x 1o'km found here for 
the extent of a velocity fluctuation lies close to the original limit of resolution, 
2x10'km. Owing to this and the irregular nature of the phenomenon, general 
similarity of the larger-scale motions is all that can be expected. ‘lhe earlier 
value obtained for the dispersion in velocity was +o:17kmsec!. ‘lhe close 
agreement of that value with the present one, +0°16kmsec'', is entirely satis- 
factory. Perhaps the most pleasing feature is, however, the conclusive 
demonstration of the stability of the velocity field from the present work. ‘Uhis 
result, foreshadowed in the earlier investigation, is of particular importance in 
that it is the only information at present available concerning the time scale of 
the motion. 

Refinements in observational technique and the accumulation of a large 
number of accurate velocity measures have thus led to a much clearer conception 
of the dynamic state of the equatorial region of the Sun. ‘This information may 
help in the identification of the physical processes causing the motion. 

4.2. Origin of the velocity field.—A simple calculation shows that the energy 
of the velocity fluctuations is not negligible compared with the flux of radiant 
energy near the surface of the Sun. ‘The energy density of black-body radiation 
is given by 407"/c, where o is Stefan’s constant, 7 the absolute temperature 
and ¢ the velocity of light; the kinetic energy per unit volume, on the other 
hand, is pu*/2, where p is the density of the material and o its velocity. ‘Taking 
108deg. K, p=1'2x10' gem", for r=0°3 at 5485A (21), and 
v=o'16kmsec"! (Section 4.1), values of 7-4 and 15-4ergem™ are obtained for 
the radiation density and kinetic energy density, respectively. While the precis« 
ratio of the two quantities is obviously somewhat uncertain, owing to 
the sensitivity of the radiation density to the assumed value of 7, it is clear 


53 


54 A. B. Hart Vol. 116 


that for any likely values of 7 and v the quantities are comparable in magnitude. 
It is interesting to see whether the magnetic energy density, »H?/87, makes any 
significant contribution. If the magnetic permeability p» is unity (6) and H, the 
magnetic field, is 1 gauss (4), the density is 0-04 ergcm™* 
would require a field of some 16 gauss. It would therefore appear that outside 
the zone associated with a sunspot (4), the energy resides mainly in the radiation 
and velocity fields. 

The question arises of what mechanism is capable of producing such a large 
quantity of kinetic energy. ‘The irregularity and large linear scale of the 
phenomenon (Section 4.1) rule out the obvious answer of convective instability. 
It appears more likely that the internal and gravitational energy of the gas are 
being converted into kinetic energy according to a suggestion made by Margules 
(15). ‘lwo contiguous columns of gas, of height A say, having a temperature 
difference AT between their bases are in unstable equilibrium and they tend to 
change to the stable configuration of the less dense gas on top, with the release 
of kinetic energy. If é is the average velocity developed, 7 the mean temperature 
of the two columns and g the acceleration due to gravity, it can be shown (12) 
that For €=o:16kmsec™', g= 2-74 x 10'cmsec™*, h=300km and 
T=5:6~ 10° deg. K, AT deg. K; for €=0-5kmsec™!, AT=7odeg.K. Thus 
the mechanism appears at least to be a possible one. 

The problem of identifying the origin of temperature differences having the 
necessary linear extent remains. From an analysis of limb-darkening observations, 
it has been shown (21) that the observable region of the photosphere consists of 
two layers. Above optical depth r= 1 (at A 5485 A) is a thermally stable layer of 
small temperature gradient in which the large-scale motions, considered in the 
present investigation, may take place; immediately below, granulation occurs 
in the form of convection cells in a narrow layer of steep temperature gradient (20). 
Beneath this again and as the ionization of hydrogen becomes important, Unséld’s 
hydrogen convection zone is encountered (25, 26). Variations in temperature 
will occur over the upper surface of this hydrogen convection zone owing to 
the differences in temperature between the ascending and descending material 
in the convection cells. ‘These inequalities in temperature will be transmitted 
to the layers vertically above, though radiative transfer in a lateral direction will 
tend to reduce them somewhat. ‘The linear extent of the observed motions may, 
therefore, give some indication of the size of the cells and hence of the thickness 
of the hydrogen convection zone. 
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II. 


Summary 


Data on the brightness and polarization of the solar corona given in a 
previous paper are analysed to separate the # and K components and to give 
the electron density in the outer corona and inner zodiacal light regions. ‘The 
electron density is supposed to be a continuous function of distance from the 
Sun with a sharp cut-off at a variable position corresponding to a sudden 
decrease of ionization. Let pf be the polarization of the F component at 
elongation 6. The following results are demonstrated. (1) If no cut-off 
is assumed and pf,,—o, then pf'<o for 6>35°. But if p¥,. is sufficiently 
large and positive, pk -o for all 6. (2) If a cut-off is imposed just inside 
the orbit of the earth and pf, o, then pf--o for 635°. The relation 
between the zodiacal light and the corona is discussed. It is shown that 
there is probably a sharp change of electron density gradient near a position 
30RO. from the Sun, 


1. [ntroduction.—It is now well established (1, 2, 3) that the light of the 
solar corona is due to the scattering of sunlight by free electrons and by dust 
particles surrounding the Sun. ‘The observed corona can therefore be considered 
as the superposition of two components, an electron (K) corona and a dust (F) 
corona. ‘The brightness of the K component decreases with increasing elongation 
because both the electron density and the flux of sunlight on the electron cloud 
diminish with increasing distance from the Sun. ‘The brightness of the F 
component also diminishes with increasing elongation; it is determined by the 
distribution of particles around the Sun, and by their average phase functions 
(which give the angular distribution of scattered light). 

The chief problem of the solar physicist is to separate the two components 
and calculate the electron density as a function of distance from the Sun. ‘The 
separation is most easily made for the outer parts of the corona by making use of 
the observed polarization, because the K component should be strongly polarized 
and the F component weakly polarized. In a previous paper (4) data are given 
for the brightness of the solar corona in the equatorial regions at sunspot minimum 
to an elongation, 6, of 13°°5 and for the polarization to an elongation of 5”, all 
at a wavelength of 6400 A. In Fig. 4 these polarization data are shown combined 
with the data obtained from zodiacal light studies by Behr and Seidentopf (5) 
which begin at elongation 35°. It is the purpose of this paper to analyse these 
data to give the electron density as a function of distance from the Sun. 

The polarization curve falls with increasing distance from the Sun, passes 
through a minimum at an elongation of about 6° and then rises to high values 
again. Because of severe technical difficulties no values for the polarization 
are available for the region 5’ <6<35', but if we are to satisfy both the corona 
and the zodiacal light measures the rise in polarization must be very steep. ‘There 
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are good theoretical and observational (6) reasons for supposing that in the 
neighbourhood of the Sun (@- 1°) the polarization of the / component is very 
small in comparison with the polarization of the A component. ‘lhe gradually 
diminishing polarization with increasing distance from the Sun in this region 
is due to the increasing predominance of the F corona over the K corona. It may 
be conjectured that the rise in observed polarization beyond 6=6° is due to the 
onset of polarization in the F component. However, it is unlikely that the 
F component is appreciably polarized at such small elongations. For very small 
particles the polarization is given by Rayleigh theory of scattering, independently 
of the refractive index of the particle material, at elongation 4, as 


1 —cos* 
) 
1 + cos*@ 
For 6=5° (the limit of the present observations), p= + 0°46 per cent. For very 


large transparent drops the polarization is readily calculated by combining the 
reflected ray with the refracted ray, taking account of their phase ditference 
This has been done by Bricard (7) for water drops with the result that for 4 © 5 
the polarization is +0-40 per cent. Polarization values for intermediate sizes 
can be calculated by the Mie theory of scattering, but they vary very rapidly 
between wide limits as the drop size is altered * and it is very difficult to obtain 
a realistic average, but the effective value is unlikely to be as high as the Rayleigh 
limit. ‘The two extreme values quoted may be assumed representative of the 
maximum likely polarization for all materials at this elongation but it cannot be 
stated definitely whether the polarization will be positive or negative. ‘These 
values are much smaller than the observed values for the corona (2°7 per cent at 
an elongation of 5). A more satisfactory explanation for the minimum at 
6~6° is therefore that the ratio of the K component to the / component has 
passed through a minimum and that as the elongation increases this ratio is 
increasing. 

2. Electron densities in the equatorial regions of the outer corona.t+-—Vhe 
procedure for the calculation of the electron density over the range of distance 
from the Sun where the polarization has been measured is simple in principle. 
At a particular distance from the Sun let the observed degree of polarization be 
Po and the degree of polarization of the electron corona only p,. At the same 
distance the contribution to the total intensity of the electron component is A 
and of the dust component is F. ‘Then, assuming that the polarization of the F 
component is negligible in comparison with the observed polarization we have 

PolPe= + F) (1) 
Now we proceed by successive approximations. We first assume plausible 
values for p, as a function of distance and so calculate A as a function of distance. 
This brightness distribution corresponds to a unique variation of polarization 
with distance which can be calculated by van de Hulst’s approximate method (3). 
If the initial guess was correct the calculated and assumed polarization values will 
be identical. If they are not, the assumed values must be changed until there is 
agreement. 

Having now decided upon the true polarization of the electron corona we can, 
in principle, pass to the electron density distribution by solving an integral 


* See, for example, the values obtainable from the tables of Gumprecht ef al. (8) 
t The problem of the electron density in the polar regions is treated in a later section 


| 
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equation. If N(r) is the electron density at distance rv from the Sun and HK(x) 
the tntensity of the electron corona at a distance x in the projected image, where 
all distances are measured in terms of the solar radius as unit and H is the mean 
surface brightness of the Sun, then we have, following van de Hulst (3), 


2 2 


K(x)=D | - =) A(r) + B(r)> = (2) 


A and B are complicated functions of r, which, for all distances from the Sun 
beyond 6R©, approximate to 1/2r* and zero respectively. At this stage, however, 
it is necessary to draw attention to a fundamental difficulty. ‘The brightness and 
polarization of the electron corona at distance x on the projected image depend 
upon the electron density at distances far beyond x=r. Supposing the model 
developed later in this paper to be correct, we find that about one quarter of the 


Fic. 1.—Variation of electron density with distance from the Sun in the region of the ecliptic 
and solar equator. 


total brightness observed at a distance of 20R© from the Sun in the image comes 
from parts of the corona which are more than 40R© from the Sun. ‘This limits 
the information obtainable about the corona at any specified distance from observa- 
tions which extend only to this distance in the projected image. _ If no information 
is available beyond x =r, this difficulty can be met (but not surmounted) only by 
representing the electron density by an explicit function of distance and assuming 
that this function does represent the electron density well at distances of up to 
r=3ry. Fortunately, the observations of the inner zodiacal light by Behr and 
Siedentopf (§) provide sufficient information about the electron density in these 
regions to allow the reduction of the outer corona observations to be carried through 
with reasonable accuracy. ‘lhe interpretation of the zodiacal light observations 
will be considered in a later section and we here only assume that at the distance 
of 128R0O* the electron densityt is of the order of 10* electrons cm~ and that the 


*1A.U.—213R 
t This value is justified in Section6. 
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electron density for the outer corona and the zodiacal light join in the way shown 
in Fig. 2. 

This curve shows that there is a rapid change of gradient of electron density 
at about 25Ro0, and if this is neglected erroneous results will be obtained in 
interpreting the polarization measurements in the outer corona. Using the 
extra data provided by the zodiacal light observations the calculations have been 


) 50 
DISTANCE FROM SUN 


2 


Electron densities in corona and zodtacal light 


TOTAL CORONA 


ELECTRON COMPONENT 
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Fic. 3.—The observed brightness of the solar corona and zodiacal light, and the computed 
intensity of the K component for three positions of cut-off of electron density at 500RO, 215RO 
and 165R©. The coronal observations are those of Blackwell(4) and the zodiacal light observations 


are those of Roach et al. (g). 
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completed by van de Hulst’s approximate method with the results shown in 
part of Fig. 1. Results of calculation for the range of observation are given in 
Table 1. ‘This table separates the observed corona into its components K and F 
and also gives the polarization of the K corona alone. ‘The K component and 


TABLI 


K and F « 


smponents of the solar corona and electron densities 


(Unit of brightness is to Bo) 


RRO | fo 
| > | | | 


oO 060, 
0535 
O40, 


043 


0404 54 d 
12 | 0°037e 0°53 | 0'070 159 145 0°66 
| 
13 | | 0066 133 | $878 | 12 


032, 


030», 


029 
Wi 0'028, 45 | 0°060 


1A O28 0°47 0°'060 69°2 


the sum K + F are plotted in Fig. 3, in which the part as far as 5 relates to the 
present discussion. We have assumed in the computations that at 5 elongation 
the polarization of the / component is zero. We have already shown that at 
this elongation |p’ | + 0-40 per cent, and the maximum error in electron density 
introduced by the forced assumption p’ =oisonly - 15 percent. ‘The assumption 
of zero polarization is made because the sign of the actual polarization is unknown. 
‘The probable error in the measurement of polarization at this elongation 1s of the 
order of + 10 per cent, giving a total uncertainty in electron density of 18 per cent. 

3. Electron densities in the zodiacal light region.— Vo continue this analysis 
into the zodiacal light region we need accurate values of the intensity and 
polarization of the zodiacal light at the wavelength already used for the coronal 
observations (6300A). ‘The best measurements of polarization are those of 
Behr and Siedentopt (§) which we use even though they were made at a wavelength 
of 5400A. We also use those measurements of the total brightness of the zodiacal 
light by Roach et al. which have been made at the wavelength 6300A (9g). ‘The 
observations of total intensity have been plotted in Fig. 3 with a plausible inter- 
polation curve linking them to the present measurements of the outer corona. 

‘The observations of Behr and Siedentopf start at an elongation of 35 . Some 
previous discussions of intensity and polarization data at greater elongations 
(for example Siedentopf, Behr and Elsasser (10)) have assumed that the polariza- 
here denoted by p} 


tion of the / component is very small for all elongations, 
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although van de Hulst (11) has considered the possibility of appreciable 


polarization of the F component. In this discussion we make no assumptions 
about except that 


This is not a drastic assumption. From the work of van de Hulst (11) and 
Bricard (7) it seems likely that |p4,.| +5 per cent. ‘This value is much less 
than the 23 per cent polarization of the zodiacal light at this elongation, but 


even if pi, is as much as +15 per cent the calculated electron density is only 


halved. ‘Ihe consequences of assuming p‘,. 40 are discussed later. 


An elongation of 35 corresponds to a distance of closest approach to the 
Sun of Ry. =213s8in35 =123R0. We now make the further assumption that 
if n, is the electron density at this distance from the Sun then 


dlogn, 
const. 


dlogr 


for all distances beyond 123R© that make an appreciable contribution to the 
intensity. ‘Ihis expression is the gradient of the graph in Fig. 2. In the graph, 
the electron density at 20R© has been calculated to be 0-26 x 10cm, and the 
electron density at 128R© is shown later to be of the order of 10%cm™*. ‘The 


gradient of electron density between these two points is small, and the only 


reasonable way of joining them without supposing a minimum in the electron 


density curve is shown. We are assuming that the curve continues beyond the 
point r= 123R© as it approached it. In addition, we shall later consider the case 
in which a sharp cut-off of electron density is imposed at an arbitrary distance from 
the Sun.* 

The brightness of the K corona can now be calculated from equation (2), 


and the polarization of K corona using equations 3, 4, 5, following 


7 A (x)- 

K(x) + K(x) (3) 

1) 

K (x) =D | N(r) ———= (4) 
re —x 


K (x)— K {x) D\N(r) (5) 


—x* 


In these equations, which have been adapted from van de Hulst, A, and K, 


represent the transverse and radial components of the polarized light. 


We first consider the case in which dlogn,/dlogr is a constant between 


r=123R© and 500R©, and has the value corresponding to interpolation between 
4OR© and 123R© (see Fig. 2). The extrapolation involved is no greater than the 


range of interpolation and the electron density tor distances greater than 500R 


does not affect the computations. Indeed, the electron density at 500R~ has 


only a very small weight in the computations. ‘This density distribution gives 


the value 


ps 0°56 


By observations 


This do« J 


udden diminution of 1onization 


not imply a sudden reduction in space density of interplanetary was, but only a 


| 
= 0°23. 
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Hence 
= 
K+F °* 
corresponding to an electron density at r= 123RO0 of 
n, = 1°08 x 10° 


The brightness of the K corona alone corresponding to this distribution is shown 
in Fig. 3 and the computed polarization of the sum of the two components, making 
the initial assumption that at all elongations p} =o (i.e. supposing only dilution 
of the K component by the F component) is shown in Fig. 4. The observations 
of Behr and Siedentopf are alse plotted. Even assuming pj =o the calculated 
polarizations are too high and to obtain agreement with observation using this 
model it is necessary to assume negative values for pf. In Fig. 5 is shown the 
negative polarization that has to be assumed. From elongation 6=35° the 
negative value increases steadily. 


Observed polarization of corona (Blackwell 1954 


/ 
pserved polarization of zodiacal light (Sredentopf ) \ 


mputed polarization, cut off at 165 Ro 
COO omputed polarization, cut off at 215 Ro «0 / 


Computed polorization ut off at 500 Ro p,*9 i/ 


ELONGATION IN DEGREES 


Fic. 4.—Observed and calculated polarization in the solar corona and zodiacal light. In the 
coronal region the points represent the model of Blackwell (4). In the zodiacal light region the 
points represent the data of Behr and Siedentopf (4). The broken lines give the calculated 
polarization of the zodiacal light for positions of cut-off of electron density at 500RQO, 215RO 
and 165RQ©) assumin py for all values of @. 


A negative polarization is not impossible but in the absence of any proof to 
the contrary it is desirable to try to find a model which needs only a positive value 
for pf. If it is supposed that p{,.=o0 this can be achieved only by imposing a 
sudden cut-off on the electron density distribution, 
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Now let us suppose a sudden cut-off at 165RO (0°78 A.U.). This gives 


= 
and 
n,= 1:20 x 10°cm~ (at r= 123RC). 
The brightness of this model K corona is shown in Fig. 3 and the computed 
polarization for pf,.=o is shown in Fig. 4. ‘The polarization of the F corona 
that must be assumed to account for the observed polarization is shown in Fig. 5, 
together with another curve for a K corona with a cut-off at 215RO. Ifa cut-off 
at 165R© is assumed, the F corona has a positive polarization, but a cut-off at 
215R© still requires a negative value for for 
We have so far assumed that p/,.=0, and we should now investigate how 
far these conclusions are changed by assuming another value for pj,.. 


+30 


Fic. 5.—Calculated polarization of F component for various models of K component, 
We now suppose p= 0-12 and no cut-off in electron density giving 


Piss = 0°57 
and 
n,=0°66 x 10% (at r=123RO). 

Even this extreme value for the F corona polarization alters the computed electron 
density very little. If plausible values of the polarization of the sum of the 
components are interpolated between the observations of the corona and the 
zodiacal light over the range 20° —4< 35° (a computed interpolation is given in 
Fig. 4) p} can be calculated for this range. ‘The results are shown in Fig. 5, 
Here maximum and minimum interpolated values for p/*+* have been assumed, 
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and the corresponding maximum and minimum values for p£ have been shown. 
It will be seen that p/ shows a maximum for elongations in the neighbourhood 
of 30°-35 and that even at = 20° the polarization cannot very well be less than 
5 per cent and may be as much as 14 per cent. Both of these results conflict 
with the expectations of theory. In particular the results of Bricard make it 
likely that |p!,.| 3 per cent. The ratio of the electron component to the 
combined components is plotted in Fig. 6 for two models. 


\. ko 


INC ATION IN 


hic, 6 Ratio K (KF) jor models of K components narving cut-off at 500R© and 165R 


and © 


4. Summary of conclusions about electron densities in the zodiacal light region.— 
We emphasize that this analysis depends upon the assumption that the intensity 
and polarization of the zodiacal light are accurately known. If the assumption 
that p‘, =o is made and no cut-off is imposed upon the electron density, then 
pio for 035. But if pi, is taken to be sufficiently high and there is no 
cut-off then pj} —o for all values of 6. However, p/ is then also unreasonably 
high for @~20. If a cut-off is imposed just inside the orbit of the Earth and 
pi, =o, then p} for An analysis of the data for the extreme cases 
of no cut-off and a cut-off at 165R©, both with pi... =o is given in Table II. 

5. The form of the electron density curve..-We should at this stage consider 
whether the derived electron density curve is a reasonable one. 

It has been demonstrated that ander certain conditions the observed data 
may be compatible only with a sharp reduction of electron density near the 
position of the Earth’s orbit. Opik (12), however, has considered theoretically 
the problem of the ionization of an interplanetary gas and concluded that the 
ionization is complete to distances far beyond the orbit of Jupiter. 

It is disturbing that at distances beyond about 30R© the density of the 
interplanetary gas seems to decrease so slowly with increasing distance from the 
Sun. An increased gradient can be obtained only by supposing that p4,,.~0-19 
but this assumption leads to an unlikely form for pj} (see Fig. 5 where derived 
values of p} are plotted for p‘,.=o-12). The derivation of the form of pj 
is the outstanding problem of the zodiacal light. At the present moment we have 
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to choose between a negative value for p! over most of the range of elongation, 
and a cut-off of electron density near the orbit of the Earth. ‘The first assumption 
is the more likely. 


Models of the outer corona and inner zodiacal light regions 


‘TABLI 


(Unit of brightness 10°'* « average intensity of solar disk) 


No cut-off Cut-off at 165 RO .(0°78 A.U.) 
Elongation K+F |— 
| 24'0 0'032 | 
10 0°087 0°039 | 
13 7°05 0°120 0°054 | 
16 3°95 0°617 0'072 
20 2°21 0°457 395 | 12 
25 0°33! 0°29 o'16 0°257 | 0°22 | 
| 30 0°246 0°33 0°175 21 
35 0°49 0°35 | 0°23 | O'24 0°23 
| 40 0°39 0'25 | O83 | 
| | | 


6. Comparison with other data.—'Vhe electron density data are plotted in 
Fig. 1. On this diagram there are also shown the data of van de Hulst’s model 
for the minimum corona (3), the results of Michard’s computations from his 
observations of the eclipse of 1952 February 25 (13) and the coronal model of 
O’Brien and Bell (14) deduced from radio observation over the period July 1951 
to July 1952. ‘The agreement with Michard’s data is within the limits of experi- 
mental error, but the electron densities of van de Hulst are probably too low. 
On the basis of Fig. 1 it is suggested that the electron density values of ‘Table ITI. 
are more representative of sunspot minimum than the model of van de Hulst. 


Tasie Ill 


Suggested electron densities in the equatorial regions at sunspot minimum 


| 3°39 10° 
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} 40 1°35 
0°96 
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‘ 


‘The larger electron densities obtained by O’ Brien may only be due to the difference 
in phase of the solar cycle. As Allen (15) has pointed out the radio results are 
not strictly comparable with the optical results because the former are a measure 
of [n,2]'*. In Fig. 2 the corona data are related to the datum obtained from 
zodiacal light observations at an elongation of 35°. ‘There is naturally good 
agreement with Siedentopf’s electron densities (§, 10) in the region near 123R 
because practically the same data have been used, but in the absence of 
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observational data about the polarization of the F component at great elongations 
it is unjustifiable to draw conclusions about the electron densities at greater 
distances. Storey (16) has derived the electron density of 400 cm~* for regions 


near to the Earth’s orbit. ‘This value is of the same order of magnitude as those 
shown here. 


7. Connection between corona and zodiacal light.—I\t has often been discussed 
whether or not the zodiacal light is continuous with the solar corona. The solar 
physicist is chiefly concerned with the electron component and he wishes to know 
whether or not the electron component of the zodiacal light can be considered 
to be part of the solar corona. ‘The F component is only an unwanted diluent 
of this component and for this reason it is impossible to decide this question 


simply by an examination of total brightness, because except very close to the 
Sun and perhaps at elongation 40 ~50° the total brightness is dominated by the 
F corona. Any arguments based on total intensity are therefore very misleading. 
But there is no a@ priori reason to expect a smooth gradation of electron density 
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Pic. 7.-Gradient of electron density in the region of the ecliptic and solar equator. 


between corona and zodiacal light. Looking at the electron densities of Fig. 2 
it is evident that there is no discontinuity in electron density, although there is 
a rapid change in density gradient near 30R0. In Fig. 7 the gradient has been 
plotted as a function of distance and there is clearly a discontinuity at this point, 
although it should be realized that the part from 20R© onwards has been derived 
from an interpolated curve. Even bearing in mind the fact that the electron 
density curve is not well defined in this region, it is difficult to avoid the conclusion 
that there is a sharp change of gradient here. We do not suggest at this stage, 
however, that this discontinuity represents a real difference in origin between the 
corona and the zodiacal light. 

8. Electron density in the polar region of the corona.—'The polarization data 
here are very uncertain because the polarization islow. Supposing (from paper I) 
that at 6R© the polarization is 0-006, the calculated electron density is of the order 
of 15x 10%cm™*, ‘This value for the electron density is an order of magnitude 
only, but it is not incompatible with the value that would be expected from an 
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extrapolation of van de Hulst’s model for the polar region at sunspot minimum 
(10?cm™). On the other hand it is considerably less than the value found from 
a reasonable extrapolation of Michard’s model (8-5 x 10*cm™*), and we suggest 
that this difference is due to the difference of phase of the solar cycle between 
the two eclipses. ‘There can certainly be no doubt that the electron density in 
the polar regions at sunspot minimum is very low in comparison with densities 
in the equatorial regions. All available electron densities are shown plotted in 
Fig. 8. Van de Hulst’s model for sunspot minimum in the polar regions may 
be taken to be a good guide to the clectron density and we suggest the addition 


of the following data (‘Table IV) 


T 


5 
DISTANCE FROM SUN 


Fic. 8.—Electron densities in the polar regions at sunspot minimum. 
TaBLe IV 


Electron densities in the polar regions at sunspot minimum 


9. Further observational work.—\t is now possible to make a reliable 
interpolation between the intensity of the corona and the intensity of the zodiacal 
light, but there is still a wide gap between the two polarization curves, the closing 
of which will greatly help their interpretation. ‘lhe practical difficulty in bridging 
this gap from the sunward side lies in the fact that the quantity to be measured 
on the photographic plate to deduce the polarization of the corona is made small 
by the presence of the sky background. ‘Thus in the work described in Paper I 
it was necessary to measure a polarization of 0°63 per cent on the plate in order 
to deduce a polarization of the corona of 2-8 per cent at an elongation of 5°. ‘This 
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work was done at an altitude of 30000 ft. Using a likely polarization curve 
interpolated between the outer corona and the zodiacal light we can now calculate 
the height at which observations should be made to give polarization data of the 
same accuracy as before at various distances from the Sun. ‘The distances and 
corresponding heights are given in ‘Table V. It can be seen from this table that, 
using this technique, it will not be possible to push the measures of polarization 
much beyond 40R© from the Sun, 
TABLE V 


Heights at which the difference between the polarization of the sky alone and the combined sky 
and corona is 0°0063 


Distance Height of 
from Sun 


| klongation Height of 


Observer Observer 


30 000 ft. 15 63 000 ft. 
| 7° +s 53 000 ft 73,000 ft. 
| 10 ‘Oo 62 000 ft 25 "0 7600c0 ft. 


It is even more desirable that the zodiacal light data should be extended 
back towards the Sun. In the region 20° <@<30° the polarization is probably 
of the order of 20 per cent and the uncertainty in interpretation introduced by 
the assumption p}, =o is a minimum. It is hoped to present some data for this 
region in a later paper. 

Acknowledgment.—'The author is much indebted to Professor R. O. Redman 
for discussion of these problems. 
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Note added in proof.—This paper inadvertently omits reference to Ohman’s 
discussion (Stockholm Observatory Annals, 18, No. 8, 1955) of his observations 
of negative polarization in the / component of the corona (6). 
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Summary 


The paper describes the occulting disk method of coronal photography 


and photometry used at the solar eclipse of 1954 June 30 Although the 
eclipse was obscured by thin cloud it has been possible to measure the corona! 
brightness, colour and ellipticity in the range 1.2 to 12 radu 


1. Introduction.—kclipses of the Sun provide the only opportunities for 
making those optical observations of the Sun’s outer corona which are required 
for studying the physics of the solar corona and interplanetary particles. Some 
of the more serious difficulties for photometric observations have been 

(a) to allow for the light scattered from the brilliant inner corona, 

(4) to allow for sky light, and 

(c) to cope with the great range of relative intensity from the inner to the outer 

corona. 
‘The apparatus and observations of the University of London Observatory eclipse 
expedition of 1954 June were designed to determine or eliminate these effects 
as well as possible. ‘lhe eclipse was observed from the island of Syd Koster, 
Sweden. 

The observations were made by photographic photometry in two colours. 
The novel feature of the apparatus was an occulting disk, a little larger than the 
camera aperture, driven at the diurnal rate about 15 feet in front of the equatorially 
mounted camera. ‘lhe disk prevents any light from the precise direction of the 
camera axis from passing into the camera. ‘lhe effective camera aperture was 
thus zero for the axis (for the umbral region in the centre of the plate), but increased 
outward (in the penumbra) and became full aperture at about a quarter of the 
full radius of the field. Very little light from the bright inner corona reached the 
camera lens and hence scattering within the optical system was reduced to a 
minimum. At the same time the change in effective aperture counterbalanced 
to some extent the change in surface brightness of the corona, 

The eclipse of June 30 was observed through thin cloud and therefore the 
photometric results cannot be as accurate as was hoped. Nevertheless it has 
been possible to determine from the photographs the distribution of coronal 
brightness from 1-1 to 12 solar radii and to give some estimate of the colour and 
ellipticity of the corona inthis range. Because of cloud absorption it was necessary 
to utilize the inner corona to obtain absolute intensity calibration. ‘The corona 
was of extreme minimum type. 

2. Equipment.—Two cameras A and B were attached to the ends of the 
declination axis of a weight driven polar axis mounting as shown in Fig. 1. Each 
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camera was provided with interchangeable gelatine filters F; Wratten 47 for 
blue, and Wratten 27 for red. ‘The stops 5 used in the lenses were considerably 
smaller than full size with the result that there was no vignetting over the field of 
the plate. Each lens was enclosed with a hood H to prevent general sky light 
from reaching it, but of course not obscuring the outer corona field. 

The trolley mechanism ‘T driving the disks was mounted on the top of a 
14-ft. wooden tower. Its motor was operated by lead storage batteries and the 
speed controlled by means of a rheostat. ‘The occulting disks D were blackened 
on the under surface. 


nd BB and occultin 


top, | filter 


‘The dimensions of the optical systems for cameras A and B are given in ‘Table I, 
where ‘‘umbra”’ and “‘ penumbra” refer to the out-of-focus shadow of the disk 
on the photographic plate. 

The reason for having two cameras was that it seemed important to obtain 
photographs with the umbral radius both greater than and smaller than the solar 
radius. Inthe former case (camera A) the inner corona is lost completely but the 
intensity variations of the outer corona are well compensated by the disk shadow. 


If on the other hand the whole inner corona is to appear, the umbral radius 


must be smaller than the solar radius. In that case the compensation is not as 
good and photographs have a large intensity range. For the camera B the 


umbral shadow was only slightly smaller than the solar disk and the intensity 
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range was satisfactory (not more than 20:1). A more realistic and informative 
photograph of the corona is obtained if the inner corona is just visible and it 
would now appear that the ideal arrangement is to have the umbra about the 
same size as the solar image. 


Optical dimensio 


Camera A Camera B 
Focal length of camera lens 67°4.cm 75°5em 
Diameter of stop cm cm 
\perture (seen from outside) 9°71 em 6°80 cm 
Disk diameter 14°50cm 10°70 cm 
Distance, disk to aperture 447 cm 433 em 
Plate size (square) 
Field diamter (to 15°0 cm) *4 
Radius of Sun (on plate) 308 mm 3°45 mm 
Radius of umbra (on plate) 3°62 mm 3740 mm 
Radius of penumbra (on plate) 18°25 mm 15°2mm 
Distance, disk to dec. axis 55 em 


‘The duplication of the cameras introduced some difficulties that were not 
entirely overcome at this eclipse. Relative to the cameras the disk system 
rotates around the camera axis at sin(solar dec.) times the diurnal rate and in 
consequence it 1s not possible to have each disk moving precisely with its camera. 
It was thought that this effect could be tolerated during the 154 seconds of totality, 
but during the adjustment period before totality an allowance had to be made for 
the misalignment of the Sun, disk andcamera. ‘This could be done approximately 
on rehearsal days by observing the somewhat obscured solar image on a ground 
glass screen, but in the partial eclipse phase this procedure became very inaccurate. 
As a result the disks are not well centred on the Sun in the eclipse photographs 
and this gives rise to asymmetry of the photometry. It has not been difficult 
to allow for this effect, but it has prevented any opportunity of observing the 
real asymmetry of the coronal radiation. 

3. Eclipse photographs. —\ welve plates were exposed during the 154 seconds 
of totality. Half of these were blue-sensitive plates Kodak I1-O exposed through 
the blue Wratten 17 filter giving an effective wave-length ot about 4500 3 The 
others were panchromatic Kodak P300 plates exposed through the red Wratten 
27 filter and giving an effective wave-length of about 6200 A, Plates were exposed 
simultaneously on both cameras and the same exposure time given to the blue 
and red plates. ‘lhe three groups of four plates (blue and red in cameras A and B) 
had exposure times of 8sec., 2sec., and 20sec. in this order. ‘I ransparency 
was improving during totality. 

The plates were developed in pairs (blue and red) in a Dobson tank (1) for 
§ minutes with Kodak D76 fine-grain low-contrast developer. ‘he usual contrast 
of coronal images was reduced in this case by both the occulting disk and the 
cloud absorption. Consequently the coronal pictures had a rather dull appearance 
but at least they gave the intended result of showing (in the longer exposures) 
the complete inner and outer corona on one picture. 


' 
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4. Photometric standardization.—\t was not convenient to expose the intensity 
calibration marks on to the same plates as the coronal photographs. Instead 
these were put on to plates from the same batch and developed in similar conditions. 

‘The characteristic curve for the plates was obtained with the help of a tube 
photometer having 12 tubes and entrance holes ranging in diameter from 0-39 to 
4:99mm. ‘The same filters were used in the photometer as with the eclipse plates. 


Blotting 
paper 
2. lrrangement for absolute intensity standardization 


Absolute intensity standardization was provided by exposing the eclipse 
camera A to sunlit white blotting paper. ‘lhe arrangement is shown in Fig. 2. 
‘The blotting paper was shaded from general sky illumination by a hood and set 
with its normal about 45° from the direction of the Sun and the camera. Rough 
tests with an exposure-meter showed that the intensity of scattered light from 
the paper was almost independent of the viewing angle and hence the surface 
brightness of the illuminated paper could readily be computed in terms of the 
Sun’s mean surface brightness and the albedo of the paper. ‘The albedo could 
have been determined by laboratory photometry if justified by eclipse success, 
but for our purpose we assume the albedo to be 0-7. ‘The surface brightness of 
illuminated white paper is about 10 times greater than the inner corona and about 
10‘ times greater than the eclipse sky. ‘l'o make plate densities comparable it was 


therefore desirable to reduce the lens diameter by nearly 100 when used with 
illuminated paper. ‘The diameter of the stop actually used was 2°4mm as 
compared with 55mm tor the corona. 


‘ests were also made to determine the variation of intensity of light reaching 


the plate with position on the plate. Geometrical considerations indicate a 
variation of intensity proportional to cos*@ (@ from Fig. 2) which is almost 
negligible. ‘lo test this, some plates were exposed on an evenly bright twilight 
sky (during a rehearsal). Photometry revealed that the variation of intensity 
was much greater, being approximately proportional to cos*@ and that the 
eclipse photographs required correction for this effect. ‘lhe effect itself is 
possibly associated with the method of plate development. 

5. Photometric measurements..-Only the four 20-second exposure plates 
obtained towards the end of totality were measured. 

All plates were photometered with a new Hilger microphotometer. Runs 
were made through the centre of the Moon’s shadow in four directions: 
(i) equatorial, 4°N of E to 4°5 of W, (ii) polar, actually in the N-S line, 
(iii) and (iv) at 45 to the NS line. No attempt was made to discriminate 
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between intensities on either side of the Sun, these being averaged in all cases. 
Likewise the two 45 runs were averaged. 

The blue photographic plates with blue filter were about twice as sensitive 
to sunlight as the red with red filter. Photographic densities were converted 
into intensities from a calibration curve in the usual way. ‘The intensity scales 
of the curves were chosen to make the exposures on sunlit white blotting paper 
give the same intensity for both blue and red plates. Any difference of blue and 
red intensity measurement thus represents a departure from solar colour. No 
distinguishable departure from solar colour was found in any part of the corona. 
For intensity measurements the blue and red intensities have been averaged. 

‘The sky was far more intense than the outer corona. ‘lo correct for sky 
light it was assumed that the light in the outer part of the field (over 20 radii 
from the Sun) was entirely due tothe sky. ‘The sky brightness was then evaluated 
using the cos” 4 distribution mentioned in Section 4, but also taking unevenness 
of sky illumination into account as far as possible. ‘he estimated sky intensity 
was subtracted from the measured intensity to give the coronal values. ‘The 
reliability of this procedure may best be judged by the agreement of measurements 
in Fig. 4. 


Fic. 3 The 


shaded area represents effective aperture of lens occulted by disk 


Within the penumbral region of the disk shadow the measured intensities 
required adjustment for the variation of effective aperture; see Fig. 3. ‘To 
compute this we note that the common area of two intersecting circles of radii 
r, and r, and with centres separated by a is 


c=r,*« + 7,78 — 4(2a*r,? + + 2r,*7,? — a* 

where cos =(a® + r,* —1r,”)/(2ar,) 

cos =(a® —1,? + 1,")/(2ar,) 
Identifying r, with the radius of the lens aperture, r, with that of the disk, and 
a with sd where is the angular distance between the relevant point in the sky 
and the line joining the centres of the aperture and the disk (which should bi 
lined on to the centre of the Sun but was not precisely so in this eclipse) and d 
is the disk-to-aperture distance, we may compute the effective aperture (77r,* — c) 
as a function of %. Measured intensities were then increased in the ratio of 
full/etfective aperture. 
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Because of errors in alignment was not precisely the same as the angular 
distance p from the centre of the Sun. In order to estimate the misalignment the 
corrections were made for several assumed differences between 4 and p. ‘The 
value that gave a symmetrical coronal intensity distribution was chosen as the 
correct one. However it was found that the results, averaged on both sides of the 
Sun, were not at all sensitive to any ambiguity in estimating the difference of 


and p. 

6. Distribution of coronal brightness.—Vhe coronal measurements of surface 
intensity are shown in Fig. 4, where the data for the equatorial, 45°, and polar 
directions have been separated to avoid superposition, ‘lhe differences between 
the curves are too small to show in the diagram, but appear in ‘Table II, where 
they are compared with data from van de Hulst’s model corona (2). ‘he sky 
intensity on the same scale as the corona (with absorption) is shown in Fig. 4. 


radii 


Coronal surface intensity 


© Blue sky intensity on same ale. 
X Red sky intensity on same scale 


Since clouds absorbed light from the corona by an amount that could not be 
determined directly, some assumption had to be made to give the results on an 


absolute scale. For this purpose the intensity of the inner corona was fitted 
to van de Hulst’s model for sunspot minimum (sp-min). It is then found that 
there is good agreement between the observations and the medel over the whole 


range of radial distance. 

lhe coronal colour may be estimated by comparing the blue and red intensities 
in Fig. 4. As no systematic difference is found over most of the corona it is 
concluded that the colour is the same as the Sun. ‘This was known already for 
the inner corona (3) but the diffraction theory (4, 5) predicted a small reddening 
for larger radii. No such reddening is detectable in the data. The polar data 
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appear to indicate the outer regions (6 to 12 radii) to be bluish but the blue plate 
microphotograms in the polar direction were rather irregular and this effect is 
almost certainly due to inaccuracy. When interpreting these observations it 


TABLe II 


Logarithm of coronal surface intensity in terms of the mean solar intensity 


Observations Model | 
Radius (van de Hulst, sp-min) 


Equatorial) 45 Polar Equatorial | Polar 


should be remembered that the sky background was considerably more intense 
for the blue than the red light and hence some systematic differences might have 
appeared, All we can conclude from the present observations is that there ts 
no severe change of colour in the outer corona as far as 10 radii. 

The ellipticity of the corona may be determined by comparing the equatorial 
and polar intensity. ‘lhe equatorial polar intensity ratio of 2-9 In the inner 


corona is typical for sunspot minimum. For the outer corona beyond 3 radii 
the intensity ratio remains about 1-<* as compared with the van de Hulst model 


of 1-1. ‘This may mean either that the intensity of the diffraction F corona is 
less than supposed or that it is somewhat elliptical as close to the Sun as 5 or 10 
radii 

7. Absorption and sky brightness.—Vhe sky brightness determined by 


comparison of plate density of the coronal sky with the exposures made with 
sunlit paper are 


Sky (in blue) = 1-g mean solar intensity 
Sky (in red) =0-g x 10°-* mean solar intensity 


‘These agree very closely with the corresponding measurements at the 1940 eclipse 
(1-9 x 10 * and 11 x 10-") (6) despite the fact that one eclipse was clear and the 
other somewhat cloudy. 


| 
| 
r | | | 
loa | 7 7°6 | 6-01 | 7°75 7°06 | 
| die | 7°65 3! 7°15 | 7°65 7°29 8-86 | 
7°30 | 7°49 8-96 8°70 | 
| | 7°16 76 8:70 | 714 8:73 | 
1°5 8-99 55 8-52 | 8°95 | 856 B43 | 
1°6 | 8-82 40 8-38 | 8:75 B42 8-32 | 
1°8 | 8-61 15 B15 8-so | 8-20 8-14 | 
2°0 8-41 9°99 9°99 8-27 8:04 8-00 | 
2°5 9°95 9°68 9°69 | 9°74 9°72 | 
| 3 | 9°66 | O44 9°47 | 960 | 9°52 | 
} 9°29 | | 9°25 | 9°20 | 
5 10°87] 10°83 10°99 | 10°95 | 10°95 
6 10°77 10°62 | 10°55 | 10°75 10°75 10°75 
| s 10°40 | 10°24) 10°20 | 10°47 | 10°44 | 10°44 | 
|} 10 11°96] 11°93 | 10°23 |} 10°20 | 10°20 ; 
| 12 | 
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By comparing the coronal brightness given in Table II with values determined 
by absolute photometry with the sunlit paper, the absorption towards the end 
of totality was found to be about three quarters of a magnitude. It is evident that 
the absorption was greater during the early stages of totality. 

8. Conclusions.—'Vhe various intensity measurements in the outer corona 
cannot be applied to physical arguments until the two components K-corona 
(electron scatter) and F-corona (particle diffraction) can be segregated. ‘This 
can be done by measurements of Fraunhofer intensities or polarization and tt 1s 
hoped that such measurements will be available in due course for the present 
eclipse. In the meantime it is advisable not to attempt any re-analysis but 
merely to look for tendencies that appear in the present data. 

‘The observations agree well enough with van de Hulst’s model (2) to ensure 
that only minor ammendments to that model will be needed. ‘lhe departures 
from the model that appear to be detectable are (a) no reddening of the outer 
corona is found, (6) the outer corona is more elliptical than expected, and (c) the 
outer corona is fainter than expected. ‘hese results suggest an F-corona some- 
what fainter than the values used in the model. 


‘The University of London Observatory eclipse expedition to Syd Koster, 
Sweden was made possible by a grant from the Joint Permanent Eclipse 
Committee. Much assistance was given the author both in setting up the 
equipment and making the exposures by the Observatory’s Chief ‘Technician, 
Mr. C. R. Spratt. ‘The London expedition shared a site selected by Professor 
R. O. Redman of Cambridge from whom much assistance was obtained. Willing 
help trom the residents of Syd Koster and the Swedish astronomers was also much 


appreciated, 


University of London Observatory, 
Vill Hill Park 
London, N.U 
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Addendum on Eclipse of 1955 June 20 

\n attempt to repeat the measurements, but with only a single camera and 
disk, was made at the eclipse of 1955 June 2oin Ceylon. However, the eclipse 
was obscured by heavy cloud. li has now been decided that if the same principle 
is to be applied at future eclipses the disk should be placed inside the camera 
since it is probably easier and safer to allow for the light scattered in the lens than 


to erect and align the travelling disk. 
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AN ATTEMPT TO DETECT DEUTERIUM IN THE 
SOLAR ATMOSPHERE 


T. D. Kinman* 


(Communicated by the Director, University Observatory, Oxford) 


(Received 1955 November 26) 


Summary 


\ search is made for the A 6561-023 deuterium line in eight centre and 
limb spectra of the Sun using an interferometer with a resolving power of 
1:2 10°. Near this waveslength there are two atmospheric water vapour 
lines at A 6560-562 and A 6561097, the former of which is shown to be blended 


with a weak solar line at A 6560°596. An analysis of the line A 6561-097 for 
centre, East and West limbs yields equivalent widths of the deuterium line 
of 0°03 + 0°5, 0'09+ 0°7 mA respectively. ‘There is, therefore, no 


evidence for deuterium on the Sun. The observational errors, however, 
are such that a weak line of equivalent width 1-1 mA (corresponding. to 
*H/'H=4~10°°) might remain undetected. de Jager’s claim to have 
detected a line of 6 mA was based on a single centre spectrum from which 
sufficient photometric accuracy could not be obtained. 


Introduction.—Vhe solar abundances of deuterium, lithium and beryllium 
are of particular interest because these elements are rapidly destroyed by thermo- 
nuclear reactions at the temperature of the solar interior. Greenstein, Richardson 
and ‘l'andberg—Hanssen (8, 9) have observed that the solar abundance ratio Be/Li 
is about 30 times the meteoric value; this is taken as evidence that mixing may 
occur between the photosphere and the depth where the Li but not the Be abun- 
dance is much reduced (about 0-5 solar radii) but no deeper. It is important 
therefore to try to detect deuterium in the Sun since Salpeter (23) calculates that 
it should have a mean life of 5 x 10% years at a depth of only o-14 solar radii 
(100000 km) and so should be absent on the mixing hypothesis. 

The atomic spectrum of terrestrial hydrogen, which contains 0-015 per cent 
*H (10, 25), shows the deuterium lines as faint satellites to the blue of the hydrogen 
lines; the separation is proportional to the wave-length of the hydrogen line and 
equals 1-785 A at Ha. ‘The first suggestion that deuterium could be detected on 
the Sun was made in 1931 by Menzel (13) who showed that an unidentified line 
of intensity -3 is listed in Rowland’s Revision close to the predicted wave-length 
of the deuterium satellite to HB. ‘The region, however, is heavily blended with 
solar lines and he could find no other coincidences of unidentified with deuterium 
lines. ‘The deuterium line which should be most intense in the visible solar 
spectrum is the 2*P°-3?D transition (A6561-023) in the wing of H« and lies close 
to two atmospheric lines. Although Claas (2), from an inspection of the single 
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spectrum of the Utrecht Atlas (4), finds no evidence for a depression greater 
than 2mA at this wave-length, de Jager (3, 4) claims that this spectrum shows a 
depression of 6mA corresponding to his predicted profile of the deuterium line. 
However, to detect so weak a line in the presence of the water vapour lines clearly 
requires more than one spectrum and for this purpose new observations have 
been obtained with high resolving power using an interferometer. It is possible 
to allow for the effect of the water vapour lines (which are formed by exponential 
absorption and have symmetrical profiles) and so derive the true wing of H« from 
the observed profile; the deuterium line, if present, should then appear as a dip 


| 


in this continuum. Although such a dip is likely to be very weak, a comparison 
between spectra of the solar centre and limb (where the solar lines will be shifted 
by solar rotation) and between spectra obtained on days of different atmospheric 
water vapour content should distinguish between a weak solar line at A 6561-023 


and the random errors of observation. 


1. Observations and photometry 


1.1. Method of observation.—Spectra of the violet wing of Ha were obtained 
using a Fabry-Perot etalon placed behind the Oxford 19°83 m. telescope (1g) and 
spectroscope (20) using ‘Treanor’s method of circular channels (28). In the 
arrangement of this method for line profile measurement suggested by Dr M. G. 
Adam, the etalon is adjusted to give a broad central heterochromatic fringe and 
then tilted slightly so that the line falls on to this fringe; in this way the diffi- 
culties of photometry around a narrow circular channel are avoided. 

‘The apparatus is identical with that used by Hindmarsh (11), A 1-5 mm 
interferometer spacer was used in conjunction with a slit width of o-20mm,; in 
this case the order-to-order scattered light (which increases with slit width and 
spacer distance) corresponds to only about 2:5 per cent. ‘lhe general scattered 
light in the primary spectroscope (20) with the small slit height used (3:2 mm) is 
calculated to be only 0-06 per cent. By mounting the interferometer with its 
wedge angle perpendicular to the dispersion, the ghost images were thrown 
vertically clear of the central fringe; thus a line height of about 3 mm was avail- 
able for photometry. Spectra of the Cd line (A6438) which were taken by 
Hindmarsh shortly after the present observations showed a resolving power of 
1:2 x 10° with the 1-5 mm spacer. 

Fifteen spectra were taken of the centre of the disk, and seven of the East 
limb and five of the West limb of the Sun. ‘The limb spectra were taken with 
the slit along the E—W line and the slit centre at the point inside the limb cor- 
responded to cos @=0°233. ‘lo calibrate the plates the etalon was removed and 
a spectrum of the centre of the disk was taken through a calibrated step-wedge 
(1) mounted just in front of the spectroscope slit. A similar spectrum without 
the wedge was taken for the correction for non-uniform illumination of the slit. 
Diaphragms in the focal plane of the camera allowed either three limb or three 
centre interference spectra to be taken on each plate in addition to the two calib- 
ration spectra. ‘lhe exposure times for all spectra on each plate were made equal 
by adjusting the slit width for the calibration spectra and by using neutral tint 
filters. ‘The backed Ilford Rapid Process Panchromatic plates used were stored 
for 24 h after exposure to equalize the fading of the latent image on all exposures 
before brush development in Parkhurst developer. 


‘ 
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Three plates which were judged to have the most suitable photographic 
density were selected for measurement: 
1. Three spectra of the centre of the disk (cos#= 1-0) on 1954 Oct. 28 at 
and 1o"20™ U.T. 
2. ‘Three spectra of the East limb (cos @= 0-234) on 1954 Oct. 30 at 
9"47™, and 1ro"o8™ 
3. ‘I'wo spectra of the West limb (cos @ = 0-233) on 1954 Nov. g at 10"05™ and 
ULT. 


1.2. Microphotometry of the plates.—Microphotometer (18) tracings of the 
central heterochromatic fringe were made parallel to the dispersion with high 
magnification ( x 48°2) using a photometer slit 0-25 mm high and 0-057 mm wide. 
In order to reduce the effects of plate grain, these tracings were made at five heights, 
0°25 mm apart, in each spectrum. A low-magnification tracing ( x 6-98) of the 
whole of each spectrum parallel to the dispersion was also made; this was used to 
find the dispersion curve across the central fringe (see Section 2.1 below). 

‘lo find the characteristic curve for each plate, tracings were made parallel to 
the dispersion through each of the seven steps of the wedge. A tracing perpen- 
dicular to the dispersion of the spectrum without the step wedge gave the cor- 
rection for non-uniform illumination along the slit. Seven points in both the 
lines and continuum in each step were measured so that forty-nine points were 
available to give a well determined characteristic curve for each plate. No wave- 
length dependence of the curves was found over the wave-length range used 
(AA6562°8-74'2), so these curves were also applicable to the measured region of 
the interference spectra (AA6558—62). 


2. Reduction of the microphotometer tracings 


2.1. Determination of wave-lengths.—'Treanor (28) has shown that in the 
central heterochromatic fringe the dispersion along the radius z from the mono- 
chromatic centre is: 


dA/dz = — Az/F? (1) 


where F is the camera focal length. If p is the distance between monochromatic 
and heterochromatic centres and A is the distance from the heterochromatic centre 
along the same radius, then z=(p +h) which on substituting in (1) and integrating 


gives : 
A\ ph 


where Ay is the wave-length at the heterochromatic centre. With sufficient 
accuracy this may be written: 


ph h* 
A=A 
0 ( I I 2 2k :) ( 3) 


In the low-magnification tracings the heterochromatic centre can be found 
from the outer heterochromatic fringes so that the distances (A) of lines from 
this centre in tracing units may be found. Four lines of known wave-length 
(shown with an asterisk in ‘Table 1) may then be used to find the unknowns p 
and Ay in eqn. (3) by least squares; the resulting dispersion curve gives the wave- 
length at A tracing units from the heterochromatic centre in the low-magnification 
tracing. 


= 
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As the position of this centre is not known in the high-magnification tracings, 
the low-magnification dispersion curve must be written in terms of distances (h*) 
measured from the ‘Tim line (A6559°580) for which a corresponding distance (H*) 
can be measured on the high-magnification tracings. ‘The ratio (H*/h*) of high 
to low magnification is known, so that the high-magnification dispersion curve 
may be found. ‘This curve, however, refers to the dispersion along a diameter 


through the heterochromatic centre, while the tracings refer to a set of chordal 


lines. ‘True dispersion curves may be found by plotting a curve of errors for each 
tracing (derived from the four lines of known wave-length) and by using this to 


correct the radial dispersion curve; any error in the assumed value of (H*/h*) is 
also corrected by this procedure. In this way the wave-length may be found for 


each point measured on a high-magnification tracing. 

2.2. Determination of the mean profile I(A). g(A).—'Vable | gives details of the 
identified lines in the region covered by the tracing and also of the deutertum 
line. ‘The first and second columns give the element and the Rowland number 
in the disk from Rowland’s Revision (22), while the last three columns give the 
wave-lengths for the three plates from Rowland’s Revision corrected according to 
1.A.U. recommendations (27) and also corrected for solar rotation and Earth’s 


motion. 


"TABLE | 


Identified lines and the deuterium line 


| 

Corrected wave-length 

klement Rowland number 

Centre East limb | West limb 

Atm. wv.* | IN 6558°154 6558°154 6558°154 

Ti 6559°567 6559°535 6559°602 

} \tm. wy | INd ? 6560°562 6560°562 6560°562 
H 6561-010 6560°978 6561°045 

| Atm. wv.* 2 6561°097 6561°097 6561 °097 | 

'HE (Ha) | 40 6562°799 65602°767 6562°334 

\tm. wy | 2 6553°525 6563°525 6563°525 

\tm. I 6564°067 6564°067 6564°067 
\tm. wv.* | 6564°212 6564°212 6564°212 


(* used to determine wave-lengths) 


‘The high-magnification tracings were measured from AA6558-8-6561°8 and 
converted to intensities by means of the characteristic curve. In the vicinity of 
the deuterium line (AA6560-2-6561-6) the measures were made at each tenth of a 
tracing square (approx. 0°013 A apart) on the unsmoothed tracings. In the region 
to the blue and red of this (AA6558-8—60°2 and AA6561-6—-61-8) a smooth curve was 


drawn through the tracings and measures of this curve were made at every fifth 


of a tracing square. No smoothing was attempted in the central region in order 
that no distortion should be introduced or weak lines smoothed out. — In the other 


regions, however, which are needed to find the continuous background it is 


advantageous to smooth out the plate grain. ‘The scale in tracing units was con- 


verted to wave-lengths as described above so that for each of the five tracings in 


each spectrum the intensity in arbitrary units was known at each measured point. 


i 
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This observed profile O(A) may be expressed as the product /(A). f(A). g(A), 
where J(A) is the profile of the continuous background (the wing of Ha), g(A) is 
the profile of any absorption lines that may be present and /(A) is the profile of the 
central heterochromatic fringe. At the three wave-length regions close to 
A6559°0, A6560°2 and A6561-6 there are no absorption lines, so that here we may 
put g(A)=1. ‘The factor f(A) may therefore be found for these wave-lengths by 
comparing the observed intensities with those of a standard profile of the H« 
wing (3). The values of f(A) at intermediate wave-lengths may be found by 
fitting a power series interpolation formula. ‘The central fringe was purposely 
made broad so that f(A) was a relatively slowly varying function and could always 
be represented with sufficient accuracy by either a quadratic or cubic equation. 
Since the tracings correspond to different cross-sections of the central fringe, 
f(A) had to be determined separately for each; the intensity in the H« wing 
including superimposed absorption lines (i.e. /(A)g(A)) could then be derived at 
the measured points in each tracing. A profile was drawn and the intensity in 
the profile tabulated at o-2A intervals; these intensities for the five tracings were 
then combined to give a mean profile in each spectrum. As the time interval 
between the first and last spectra on each plate did not exceed twenty minutes, 
little variation in the intensities of the atmospheric water vapour lines may be 
expected between different spectra on the same plate. No systematic difference 
between these spectra could be seen and so they were combined to give a mean 
profile. 

‘The intensity at each point in a spectrum is the mean of five tracings; the 
r.m.s. deviation of a single tracing from this mean was calculated for about fifteen 
points in each plate and the average value of this r.m.s. deviation (0,) is given for 
each plate in the second column of ‘lable II. ‘The intensity at each point in the 
mean profile in each plate is the mean of two or three spectra; the r.m.s. deviation 
of a single spectrum from this mean was calculated for about fifteen points in 
each plate and the average value of this r.m.s. deviation (o,) is given for each plate 


II 


Average r.m.s. errors expressed as a percentage of the true continuum 


Plate 


Centre of disk (3 spectra) 0°73 per cent 0°54 per cent 


East limb (3 spectra) 1°00 per cent o'48 per cent 


West limb (2 spectra) 0°'gO per cent 0°42 per cent 


in the third column of Table II. It is seen that o, is significantly smaller than a, 
as would be expected if the systematic differences between the spectra on the 
same plate are appreciably smaller than the random errors in each spectrum, 
‘Taking a mean profile from the spectra on each plate is therefore justified, 


6 
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2.3 Determination of the mean profile g(A).—Our aim is to detect any depression 
in the Ha wing which remains when we have allowed for the water vapour lines 
(A6560°562 and A6s561-097). ‘To do this we must find g(A), the observed profile 
of these lines, then allow for the absorption profile of the water vapour lines and 
finally examine the nature of any residual absorption. 

We obtain g(A) by dividing /(A) g(A) obtained above by the standard Ha wing 
J(A) ; the resulting profile is thus independent of any assumption about the function 
I(A) except that both it and f(A) vary smoothly. The only additional assump- 
tion is that g(A) is unity (no absorption lines) in the three regions A6559°0, 
A6560°2 and A6561°6. It is possible that faint absorption lines may be present 
although a search in these spectra gave no indication of them. ‘The systematic 
error from this cause at either A6560-2 or A6561°6 is unlikely to exceed the average 
inter-spectrum scatter (o,) of o-5 per cent and this would lead to a spurious 
difference in intensity between the red and blue wings of the water vapour line 
A6561-097 of only o-1 per cent, which is well within the random errors of these 
observations. 

Fig. 1 shows the profile of g(A) obtained for the centre of the disk. 


T T 
090) 
OBO | 
6559 6560 656! 6562 


WAVELENGTH IN A 


Fic. 1.—The intensity profile g(A) for the centre of the solar disk. 


2.4. Scattered light.—For the purpose of detecting a weak line, only the 
distortion produced by the variation of the scattered light across the observed 
profile is important. Hindmarsh (11) has shown that the order-to-order scat- 
tered light for a 1-5 mm spacer is 2°25 per cent at the hetochromatic centre and 
varies as the inverse square of the distance from the adjacent order. Since the 
whole measured range of g(A) extends over only one fifth of the interorder distance 
near the heterochromatic centre, the variation of the scattered light across g(A) 
must be negligibly small. 
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3. Discussion of the observed profile g(A) 


£ 

In this profile the region of the deuterium line contains two blending lines 
(A6560°562 and A6561-097) which it is convenient to consider separately. 

3.1. A6560°562.—This line although identified as being due to atmospheric 
water vapour in Rowland’s Revision, is listed as Nd? (‘‘ broad for its intensity ”’ ; 
““a double whose components are difficult to separate’’ (21)). The present 
observations confirm that the line is broad and also that its wave-length is variable ; 
this suggests that the water vapour line is blended with a solar line. ‘The com- 
ponents of this blend will be weak and therefore on the linear part of the curve of 
growth, so it may be assumed that the equivalent width (IW) of the w.v. com- 
ponent is k, times that of (W,,,) of the w.v. line (A6561-097) and that the equivalent 
width (W,) of the solar component is A, times that of the Ti 11 line (A65 59-580). 
The w.v. component is formed by exponential absorption, so that if 7,(A) and r,(A) 
are the respective residual intensities at wave-length A of the two components 
when unblended, the equivalent width (W,) of the blend is given by 


which on expansion gives: 

If we assume that both components have a Gaussian profile of half-width o, then 
eqn (1) bécomes :— 


664k, 
W, =k, 0°664k, W. 


(2) 


Since the third term (e) on the right-hand side of eqn. (2) is small, it may be 
ignored in a first approximation (i.e. the equivalent width of the blend is taken 
to be the sum of the equivalent widths of its components), and with this equation 
of condition a least-squares solution for k, and k, is made. From these pre- 
liminary values of k, and k, and by taking o = 0-13A (the mean value for A6561-097 
and A6559°580), the third term (€) may be found; this term can now be included 
in the equation of condition and a final solution made for the factors k, and k, and 
hence for W, and W,. Values of these equivalent widths for the three plates are 
given in ‘l'able Illa; the last column gives the observed minus calculated residuals 
for the blend, which are seen to be satisfactorily small. ‘The plate-to-plate varia- 
tion of the atmospheric component is due to the variation in atmospheric water 
vapour while the solar component shows the centre-limb variation. 


(a) Equivalent widths in mA 


observed 


calculated 6560°562 W, 
Plate 


obs-calc 
W; f 


Ww, 


49°06 8:04 | 45°31 13°75 
East 28°54 58°45 21°48 51°75 13°25 3°23 61°77 3°32 
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(b) Wave-lengths in A 
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Plate Blend wv comp. | solar comp. blend blend 
Aobs Acale Acale Acale Aobs-Acale 
Centre 0'013 6560°509 | 6560°499 6560°583 6560°514 0°005 
Kast 0°045 "513 499 551 "510 0°003 
West 0'022 538 ‘499 618 "540 0°002 


The shift of the solar lines (A) is known for each plate (Table IIIb, column 
one) so that knowing the intensities of the two components, their wave-lengths 
may be deduced It is assumed that the observed wave-length centre of gravity 
of the blend(A,, ,, ) is the sum of the centres of gravity of the components weighted 
according to their equivalent widths. For this purpose it is necessary to use 
the approximate values for which (W,+ W,) equals W,, and this gives for the 
equation of condition A,,W, — (Ag+ A)W,=A,,., for each plate. The least- 


_ squares solution for the wave-lengths of the water vapour component (A,,), the 


solar component (Ax) and the resulting blend are given in columns four, five and 
six of ‘Table IIIb; the observed minus calculated residuals for the wave-length of 
the blend are given in the last column. 

The smallness of both the equivalent width and wave-length residuals in 
these approximate calculations confirms the hypothesis that A6560-562 is a blend 
of an atmospheric water vapour with a solar line. ‘The only line near to the solar 
component (A6560°596) which is listed in the Revised Multiplet Table (15) is the 
4p*D,-7d4*F,, transition of Sil observed by Kiess (12) in an arc spectrum at 
A6560°68. While the coincidence in wave-length is within the error of the labora- 
tory wave-length of o-1A estimated by Kiess, much stronger lines of Si 1 of similar 
excitation potential are not found in the solar spectrum, so that the identification 
with Sit is highly unlikely. 

The complexity of A6560°562 and in particular the presence of a solar line 
make it unsuitable for further analysis to detect the deuterium line in its red wing. 

3.2. A6561-097.—This line isan overtone line of the P branch (A6524 band) of 
the fundamental vibration spectrum of atmospheric water vapour (7) and lies in 
the region where the red wing of the deuterium line may be expected. Ideally 
the deuterium wing might be detected in the water vapour line by an analysis 
similar to that used for A6560°562. Since, however, the equivalent width of the 
deuterium wing cannot exceed a few milliangstroms, the comparison water vapour 
line would need to be unblended to a much higher accuracy and it cannot be 
certain that any of the neighbouring water vapour lines (obtainable with a simul- 
taneous exposure) have the required purity. ‘The deuterium line will be stronger 
in the blue than the red wing of A6561-097, however, so that it should be possible 
to detect the deuterium line by analysing the symmetry of the observed line, if 
the water vapour line can be taken to be symmetrical. ‘The Doppler half-width 

of the water vapour line will be approximately o-02A, while the intrumental half- 
width found by Hindmarsh is 0-o51A, and consequently the observed half-width 
of o-105A must be due to the combined influence of the interferometric profile 
and collisional broadening. Van de Hulst (29) has shown that asymmetry effects 
predicted by the Lenz formula are not found in the atmospheric oxygen vibration 
lines (e.g. Panofsky (16)), and Elsasser (5) considers that there should be little or 
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no asymmetry in water vapour lines resulting from non-electronic transitions. 
It is, therefore, legitimate for such a weak line (equivalent width 10-20 mA) to 
assume that its profile is symmetrical. 

Let the observed profile g(A) of A6561-097 be considered as a symmetrical 
water vapour line «(A,’ + AA) superposed on a background f(A,’ + AA) so that: 


+ AA) = + AA)B(Ag’ + AA) 


and A,’ is the wave-length centroid of «(Ay’+ AA). Now if Ay is the observed 
wave-length centroid of g(Ay + AA), a function s(A, + AA) can be found which is 


+ AA) = + AA) + — AA)}. 


Since B(Ay’ + AA) is very nearly unity, Ay and A,’ will only differ by about 0.001 A, 
so that for these purposes it is permissible to put Ay =A,’ and obtain 


+ AA) = + AA){B(Ag + AA) + (Ay — AA)}. 


The quotient Q(A,+ AA) obtained by dividing g(Ay+AA) by s(Ay+AA) may be 
calculated from the observed profile for each wave-length; it equals 


2B(A + AA) 
+ AA) + B(Ay ‘i AA) 


and is seen to be a function of the unknown background only. 

‘The deuterium line will be broad and shallow, because the combined Stark 
and Doppler etfects will produce a broadening of its core comparable to that in 
Ha, so that for simplicity the deuterium wing may be considered as a straight line 
of slope m. ‘lhe background profile containing this hypothetical line is shown 
in Fig. 2 and has the form: P(A, + AA) = 1-0 + m(A— A,) in the wave-length region 
(A<A,) which is being considered. For this background, O(A, + AA) is a straight 
line equal to unity at A, and of slope M where 


m 


M= 2 m to sufticient accuracy. 
m + 1-0 


Now M can be found by passing a straight line through the calculated points of 
O(Ay + AA) by least squares and so the slope m of the background may be found. 
These results are given in ‘Table IV for the three plates: column one gives the 
plate, column two gives the slope m as a percentage of the continuous background 
per A, column three gives the r.m.s. error (0,,) of this slope, while column four 
gives the maximum slope (m,,,,) corresponding to a depression which could lie 
within this error (m,,,,=m-+9,,). ‘These background slopes are very small, vary 
at random from plate to plate and are well within their r.m.s. errors; they are 


‘TABLE TV 


W “Hi 
Plate m Cu Mmax | (MA) (mA) 


Centre 
East Limb | + 0°94 1°00 | 0°03 x 107° 
West Limb 


| 
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(b) Wave-lengths in A 


Plate \ Blend wv comp. | solar comp. blend blend 
Aobs Acale Acale Acale Aobs-Acale 
Centre 0'013 6560°509 | 6560°499 6560°583 6560°514 0°005 
Kast 0°045 "513 499 ‘551 ‘510 0°003 
538 "499 615 ‘540 0'002 


‘The shift of the solar lines (A) is known for each plate (Table IIIb, column 
one) so that knowing the intensities of the two components, their wave-lengths 
may be deduced _ It is assumed that the observed wave-length centre of gravity 
of the blend(A,, ,, ) is the sum of the centres of gravity of the components weighted 
according to their equivalent widths. For this purpose it is necessary to use 
the approximate values for which (W,+ W,) equals W,, and this gives for the 
equation of condition A,,W, —(As+ A)W,=A,,., W; for each plate. ‘The least- 
squares solution for the wave-lengths of the water vapour component (A,,,), the 
solar component (Ax) and the resulting blend are given in columns four, five and 
six of Table IIIb; the observed minus calculated residuals for the wave-length of 
the blend are given in the last column. 

‘The smallness of both the equivalent width and wave-length residuals in 
these approximate calculations confirms the hypothesis that A6560-562 is a blend 
of an atmospheric water vapour with a solar line. ‘The only line near to the solar 
component (A6560°596) which is listed in the Revised Multiplet Table (15) is the 
4p*D,-74°F,, transition of Sil observed by Kiess (12) in an arc spectrum at 
A6560°68. While the coincidence in wave-length is within the error of the labora- 
tory wave-length of o-1A estimated by Kiess, much stronger lines of Si 1 of similar 
excitation potential are not found in the solar spectrum, so that the identification 
with Sit is highly unlikely. 

‘The complexity of A6560°562 and in particular the presence of a solar line 
make it unsuitable for further analysis to detect the deuterium line in its red wing. 

3.2. A6561-097.— This line isan overtone line of the P branch (A6524 band) of 
the fundamental vibration spectrum of atmospheric water vapour (7) and lies in 
the region where the red wing of the deuterium line may be expected. Ideally 
the deuterium wing might be detected in the water vapour line by an analysis 
similar to that used for A6560°562. Since, however, the equivalent width of the 
deuterium wing cannot exceed a few milliangstroms, the comparison water vapour 
line would need to be unblended to a much higher accuracy and it cannot be 
certain that any of the neighbouring water vapour lines (obtainable with a simul- 
taneous exposure) have the required purity. ‘The deuterium line will be stronger 
in the blue than the red wing of A6561-097, however, so that it should be possible 
to detect the deuterium line by analysing the symmetry of the observed line, if 
the water vapour line can be taken to be symmetrical. ‘The Doppler half-width 
of the water vapour line will be approximately o-02A, while the intrumental half- 
width found by Hindmarsh is 0-o51A, and consequently the observed half-width 
of o-105A must be due to the combined influence of the interferometric profile 
and collisional broadening. Van de Hulst (29) has shown that asymmetry effects 
predicted by the Lenz formula are not found in the atmospheric oxygen vibration 
lines (e.g. Panofsky (16)), and Elsasser (§) considers that there should be little or 
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no asymmetry in water vapour lines resulting from non-electronic transitions. 
It is, therefore, legitimate for such a weak line (equivalent width 10-20 mA) to 
assume that its profile is symmetrical. 

Let the observed profile g(A) of A6561-097 be considered as a symmetrical 
water vapour line a(Ay’ + AA) superposed on a background f(A,’ + AA) so that: 


+ AA) = + AA)B(Ag’ + AA) 


and A,’ is the wave-length centroid of «(A,’+ AA). Now if Ay is the observed 
wave-length centroid of g(Ay+ AA), a function s(Ay + AA) can be found which is 


+ AA) = + AA) + (Ay — AA)}. 


Since B(Ay’ + AA) is very nearly unity, Ay and A,’ will only differ by about 0.001 A, 
so that for these purposes it is permissible to put Ay =A,’ and obtain 


+ AA) = + AA){B(Ag + AA) + AA)}. 


The quotient Q(Ay+ AA) obtained by dividing g(Ay+ AA) by + AA) may be 
calculated from the observed profile for each wave-length ; it equals 


2P(Ay + AA) 
t AA) t AA) 
and is seen to be a function of the unknown background only. 

‘The deuterium line will be broad and shallow, because the combined Stark 
and Doppler effects will produce a broadening of its core comparable to that in 
Hz, so that for simplicity the deuterium wing may be considered as a straight line 
of slope m. ‘The background profile containing this hypothetical line is shown 
in Fig. 2 and has the form: f(A, + AA) = 1-0 + m(A— A,) in the wave-length region 
(A<A,) which is being considered. For this background, O(A, + AA) is a straight 
line equal to unity at A, and of slope M where 


m 

M= ————. = m to sufficient accuracy. 

m (Ay— Ay) + 1-0 
Now M can be found by passing a straight line through the calculated points of 
O(Ay + AA) by least squares and so the slope m of the background may be found. 
These results are given in ‘Table IV for the three plates: column one gives the 
plate, column two gives the slope m as a percentage of the continuous background 
per A, column three gives the r.m.s. error (¢,,) of this slope, while column four 
gives the maximum slope (m,,,,) corresponding to a depression which could lie 
within this error (m,,,,=m+0,,). ‘Vhese background slopes are very small, vary 
at random from plate to plate and are well within their r.m.s. errors; they are 


TABLE TV 


| =| (7H | 
Plate m | (mA) ( | (mA) | 
Centre 0°06 | + 1°10 1°04 | 0°00 1°07 4°0~ 10 
East Limb 0:06 | + 0-94 1700 | 0°03 0'2 10 | 4'4* 10 
West Limb | + 1°27 1°42 | 0'7 10 | 0°79 10 
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due doubtless to the small photometric errors in the observed profiles. ‘There 
is therefore no evidence for a deuterium line in this wave-length region. 

de Jager (3, 4) has computed the deuterium profiles which would be produced 
by different abundance ratios (7H/!H) and these may be interpolated to calculate 
abundances from the present observations. It is assumed that the straight line 
profiles derived above approximate to the de Jager profile which has the same 
equivalent width over the wave-length range observed ; the equivalent width (W) 
of the whole deuterium line and the abundance ratio (7H/'H) which correspond 
to a given slope (m) of the background can then be found. Columns five and six 
of ‘Table IV give W and *H/!H for the mean values of the slopes derived for 
each plate, and columns seven and eight give the same quantities for the maxi- 
mum depression which could be present and be just within the r.m.s. error. ‘The 
mean abundance from the three plates is only 0-3 x 10~° and lies well within the 
errors of observation. ‘The deuterium line should be strongest in the centre 
spectrum and here the strongest line which could just remain undetected would 
have an equivalent width of 1-1 mA, corresponding to a maximum abundance 
ratio of 4 10°°. We may compare our observed abundance ratio of 0-4 x 107° 
and the maximum observable value of 4 x 10~° with de Jager’s claim to have found 
a line of equivalent width 6mA corresponding to an abundance ratio 14°5 < 107° 
from the inspection of a single spectrum. 


A 


Vic. 2.—Background profile fi (Ag+ 4A) containing the hypothetical deuterium line. 


Severny (24) claims to have found a 100mA asymmetry between the red and 
blue wings of Ha which he attributes to deuterium. He assumes that his appar- 
atus function is symmetrical because he observes that other strong solar lines 
including the D, sodium line are symmetrical ; it is known, however, that D, is 
asymmetrical on account of a strong water vapour component in its violet wing. 
Previous photographic (6, 26) and photoelectric (17) observations of H« give no 
evidence of any asymmetry and Severny’s claim for a high deuterium abundance 
must be treated with reserve. 

The deuterium line, if present, must be so weak that any attempt to detect it 
requires many independently reduced spectra; it is found from a series of such 
spectra of both solar centre and limb taken on days of different atmospheric 
water vapour content that no deuterium line is present. It is possible, however, 
that a line exists which is weaker than the limit set by photographic photometry, 
which corresponds to an abundance ratio (7H/'H) of 4 10°°. Any attempt to 
detect such a so weak line will require a new photometric technique. 
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ASTROMETRIC OBSERVATIONS OF EROS 
DURING THE 1954 OPPOSITION 


v. B. Lourens 


(Communicated by H. M. Astronomer at the Cape) 
(Received 1955 November 8) 


Summary 

‘Table I gives the astrometric positions on the FK3 System for the equinox 
of 1950°0 of Eros on 45 occasions as observed at the Cape during the 1954 
opposition 


At its opposition in 1954, Eros, was too far south to be conveniently observed 
from any northern observatory. An effort was therefore made at the Cape to 
secure accurate observations of its position at approximately four day intervals 
from the end of January to the end of July. ‘The observations were made with 
the 24-inch Victoria Refractor (scale 1 mm=30":64). ‘The telescope was used 
west of the pie, and the plates taken, as far as possible, when Eros was close 
to the meridian and thus at a relatively small zenith distance. Ilford Astronomical 
Zenith plates were used and in general four exposures were made on each plate. 
The exposure times varied from five minutes to one minute depending on the 
brightness of Eros, which even at its brightest (11-7) was nearly two magnitudes 
fainter than the reference stars. Stars were used for guiding so that on the 


longer exposures the images of Eros are slightly trailed. 

The plates were measured with the Hilger Long-Screw Machine specially 
designed for astrometric work. Each plate was measured by two ind pendent 
measurers, each measurer measuring both forward and reversed. ‘Ihe reference 
stars for reducing the plate measures were selected from the following catalogues : 

Yale Zone Catalogues, — 22° to — 27° and —27° to —30,, Epoch 1933. 

Cape Zone Catalogues, — 30° to — 35° and —35° to — 40, Epochs 1932 and 

1936. 

Cape Astrographic Zone Bright and Faint Star Catalogues, Epoch 1goo. 
To reduce the effect of uncertainties in the proper motions which are involved 
in bringing the catalogue places up to the epoch of the Eros plates, as many 
reference stars as could be found in the catalogues were used, though a few had 
to be rejected on account of poor distribution over the plate. Altogether 544 
reference stars were used for the 45 plates, an average of 12 per plate, though the 
actual number on a single plate ranged from 7 to 16. Standard coordinates for 
these stars were computed for each plate, the unit for both the calculations and 
the measures being 0-0001 mm. 

The ordinary linear plate constants were found by the Dyson—Christie method 
and were used both to reduce the measured position of Eros and to calculate the 
residuals (observed coordinate computed coordinate) for the reference stars. 
These residuals are relatively large for those stars for which the proper motion 
is poorly determined. A second solution was therefore made for each plate 
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using the second set of measures and the observed positions of the reference 
stars deduced from the first set of measures and reductions. In a few cases 
where the distribution of the catalogue stars was poor, extra stars whose 
coordinates had been determined from the first set of measures were used for the 
reduction of the second set. ‘lhe residuals from this second solution are small 
and presumably arise principally from the error of measurement. 

The residuals obtained from the first solution were analysed to investigate 
the effect of the second and third order terms on the linear plate constants and 
to compute the necessary correction functions. ‘lhis analysis indicated that the 
plates were very slightly tilted in the right ascension direction and that the cubic 
scale term is appreciable. ‘There were insufficient data either to establish the 
existence or to evaluate the presence of a colour-magnitude effect; nor could 
such an effect be eliminated by the choice of reference stars of the same colour 
and magnitude as Eros. ‘The mean photographic magnitude of the reference 
stars used is 9:7 and their mean colour index +0°47, though the means for the 
individual plates range from 8-7 to 10-6 and from —o-o1 to +092. According 
to the ephemeris, the magnitude of Eros varied from 12-7 at the end of January 
to 11°7 at the time of opposition on May 14 and then to 13°8 by the end of July. 

The positions finally deduced for Eros are given for each plate in the third 
and fourth columns of Table I. ‘They are the mean places resulting from the 
two independent measures of all the exposures on the plate to which have been 
applied the plate corrections deduced from the analysis of the residuals of the 
reference stars and the systematic corrections needed to reduce the catalogue 
places of the reference stars to FK 3 system. ‘The U.'T. of the mid-time of the 
exposures is given in the second column of ‘lable | and the estimated probable 
error of the given position, assumed to be the same in both coordinates, is given 
in the fifth column. ‘This error includes both the errors of measurement and 
of the assumed positions of the reference stars. ‘The sixth column gives the 
actual difference between the results obtained from the first and second set of 
measures. 

These differences result mainly from the two measurements of the Eros 
images since the same reference stars were used for both solutions and the plates 
were centred on Eros. ‘hey indicate + 0"-035 in each coordinate as the probable 
error of the mean of the two sets of measures of the Eros images. ‘lhe corre- 
sponding probable error for a star, as found from the differences of the two 
measures of each plate, averages + 0”-040 for all the plates, but varies with the 
exposure time and clearly shows the decrease in precision of measurement as 
the size of the image increases. ‘The average value for plates with one minute 
exposure is + 0”-034, sufficiently close to that found for the much fainter Eros 
images to indicate that this is approaching the optimum precision that can be 
obtained from this series of plates. 

The probable error of a single reference star as determined from the residuals 
of the first set of measures is +0”-309 for all the plates, but there is a variation 
with the catalogue from which the assumed places were derived. ‘lable II] shows 
the average probable error corresponding to the different catalogues used. 
These probable errors indicate mainly the relative error of the proper motion 
contained in the positions of the reference stars at the epoch of the Fros plates. 
It will be noticed that the C.A.Z. proper motions, which were obtained by the 
differential measurement of two sets of photographic plates, are markedly 
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1954 
28 °09473 
01°09651 
11390 
12328 
‘12817 
12000 
12140 
‘12279 
14594 
°14227 
12631 
11101 
"13822 
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Equinox 1950°0 and FK3 system 


Dec. 


29 
15 
32 
47 


36 
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16°38 
47°14 
35°41 
12°40 
11°30 
28°85 
43°23 


‘09 
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14 
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18 
13 
13 
‘99 
‘10 
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Eros Measure 
(1st minus 2nd) 


‘027 
‘018 
‘006 
‘O12 


10 
03 
‘03 


‘07 


"20 
‘10 
‘02 
04 
03 
04 
05 
15 
‘08 
13 
‘Ol 
07 


0°04 
0°32 


superior as regards relative accidental errors to those in the Zone Catalogues 
which were obtained by comparing a photographic place with an earlier meridian 
place. Some information about the systematic errors of the proper motions will 
be obtained when the observed positions of Eros are compared with a definitive 


ephemeris. 


R.A. 
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: 14 23 03°120 26 26°33 
14 33 26-917 28 12°00 “16 
14 41 09°706 29 22°18 
14 48 44°05! 30 06°00 18 
15 03 25°55! 33 99 34°57 208 19 
i 15 21 §6°320 36 05 47°69 "09 209 
7 15 28 30°719 37 08 09°77 ‘09 205 19 
15 36 53°026 38 28 10°20 09 103 
15 44 46°705 39 45 05°97 205 
1§ 52 0O1°527 40 58 06°87 
ae 15 58 34°925 42 07 33°49 ‘08 204 
Se 19 03 00°683 42 57 26°70 ‘09 104 
16 07 02°165 43 46 11°11 ‘09 204 
16 11 33°380 44 47 26°51 ‘09 03 
eee 16 14 19°372 45 31 31°82 06 208 
ou 16 17 32°234 46 39 41°41 ‘06 202 
~~ 16 18 46°103 47 29 37°90 05 105 
16 18 42°142 48 14 53°95 
ne 16 17 17°262 45 54 13°40 08 006 
16 14 30°195 49 26 46°67 ‘06 13 4 
16 O1 43°495 50 10 41°29 ‘08 10 12 
15 58 18-978 5° 11 42°60 06 21 04 
15 38 13°660 49 3° 39°03 08 02 
eee 15 29 41°412 48 52 10°27 ‘07 103 ye) 
15 11 22°956 46 42 08-00 06 08 14 
15 09 34°563 46 24 51°21 ‘08 605 D2 
15 07 48°166 49 06 54°05 ‘09 202 re) 
15 02 50°174 45 10 31°90 ‘08 
= 14 54 35°345 43 08 50°37 ‘10 208 D6 
14 51 24°579 05 57°97 07 08 23 
14 48 49°215 41 OI 52°00 ‘10 04 D3 
14 46 23405 —39 38 35°95 04 D3 
14 44 59°990 38 17 39°28 "10 122 0 
14 44 34°789 37 00 20°78 902 4 
14 44 37°012 41 39°55 004 D4 
14 45 03°192 35 47 202 27 
14 47 15°028 34 08 10 3 
14 48 58:236 33 22 205 
a 14 §1 49°539 32 27 09 4 
15 00 18-703 30 41 4 
15 O1 24°312 30 31 104 6 
gees 15 03 40°373 30 12 102 6 
15 09 48-982 29 28 01 
15 15 09°147 28 57 07 
ae 15 19 24°478 28 a 12 


of Eros during the 1954 opposition 
II 


No. of Mean A 
Plates P.E. Epoch 
years 


Mean P.E. 
AEpoch 


Catalogue 


Yale 27 + 0°369 21 + 0°0176 
Yale 30 5 0°357 21 170 
Cape —: 0°472 22 215 
Cape — 35 40 9 0°306 18 170 
CAZ to —52 22 + 0'238 55 + 0°0043 


My best thanks are due to Mrs M. C. Coetzee, Mr W. A. Rasmussen and 
Mr 'T. W. Russo for their assistance in the measurements and the reductions 
and to Dr R. H. Stoy and Mr J. B. G. Turner for their help in securing the plates. 
It is hoped that a detailed account of these observations will be published in the 
Annals of the Cape Observatory. 


The Cape Observatory: 
1955 October 18. 
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Summary 
Radio echoes have been observed in two narrow-beam aerials over the 
period October 1949 to September 1951. An analysis of the hourly rate of 
echoes has given a measure of the radiant points of sporadic meteors. 
The radiants are found to be concentrated towards the plane of the ecliptic in 
the directions of the apex, Sun and antihelion points. ‘The majority of meteors 
in space are moving in direct orbits about the Sun, and there is a dense 


complex of orbits in the region traversed by the Earth during the months May 
to August. 


1. Introduction.—Both stream and sporadic meteor particles occur in the 
space surrounding the Earth. ‘These particles evaporate on colliding with 
the atmosphere, generating heat, light and ionization. Stream meteors move 
in paths that are almost parallel, and their radiants lie in a small area on the celestial 
sphere, but sporadic meteors are dispersed so that the radiants are distributed 
over the entire sphere. 


Although there was no observational evidence, for many years it was assumed 


that the motions of sporadic meteors in space were completely random. ‘This 


Fic. 1.—The geometry of the radio echo reflection from a meteor trail 
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led to the erroneous idea that most meteors were interstellar in origin, as discussed 
by Lovell (1954). Subsequent observations have been made photographically 
(Whipple, 1954) but the number of sporadic radiants was insufficient and, of 
course, no daytime meteors could be included. ‘This paper describes an investiga- 
tion of the distribution of sporadic meteor radiants on a statistical basis, using 
the radio echo method of observation. 
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Fic. 2 a.—Kelation between time of transit of en active radiant and the time of appearance of 
echoes in each aerial. 
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Fic. 2 b.—Relation between declination of an active radiant and the time of appearance of echoes 
in each aerial, 


2. Apparatus and method.—The equipment, situated at Jodrell Bank, Cheshire, 
Long. 2° 18’ W, Lat. 53° 14’ N, was originally designed for the measurement of 
stream radiants and has been described by Aspinall, Clegg and Hawkins (1951). 
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Radio echoes were obtained when a meteor formed an ionized column, as shown 
in Fig. 1, in a plane perpendicular to a pencil beam aerial system. When the trail 
was projected on to the celestial sphere ABCD, the radiant point was somewhere 
on the great circle EFG. A similar aerial (not shown in the figure) detected 
echoes from radiants on the circle HFJ. Echoes were produced in each beam in 
turn as the celestial sphere rotated and carried a group of active radiants across 
these circles. In this way meteor activity in the northern hemisphere was 
sampled once in every 24 hour period. If echoes appeared at time 7’; in aerial 1 
and at time 7’, in aerial 2 then the time of transit was a function of 7, and T, 
as shown in Fig. 2a; the declination could be found from the time lag 7, — 7, 
as shown in Fig. 25. 

Because the aerial beam had a finite vertical width, echoes were obtained 
at an angle a above the centre of the beam from radiants on the great circle EF’G. 
‘he sector in which radiants could be observed was enlarged in a similar way 
by the horizontal beam width. If an element of area on the sphere supplied N 
radiants per hour from which N’ echoes were actually detected, the sensitivity, 
S, of the equipment was defined as 

S=N’/N. (1) 

‘The sensitivity depends on the position (a, 6) of an element within the sector, 

where 6 is the angular distance from G, and for a narrow beam aerial we may 
write 

S=f,(a) fo(6). (2) 

If we assume that N remains constant as a group of radiants crosses the sector, 

the function f,(a) is given by the range distribution of observed echoes, and 
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Fic. 3.—Plan view of the celestial sphere showing the sensitive sectors of each aerial. 
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f(b) may be evaluated if we observe a stream radiant at different values of 6 with 
a moveable aerial. This has been done by Hawkins (1952), and normalized 
values of S for aerial 1 are plotted in Fig. 3, which is a plane projection of the 
celestial sphere centred on the zenith. ‘These values are the same for aerial 2 
except that the contours are displaced through a clockwise angle of 49°-2 about 
the zenith. ‘The function /,(6) is approximately equal to sind. It should be 
noted that these radiant sensitivity contours, drawn on the celestial sphere, 
differ from those used by Clegg (1948) in which the aerial polar diagram was 
projected on to a layer g5 km above the Earth’s surface. In the earlier method 
the time of passage of a radiant through the aerial beam was taken from the 
appearance of echoes at maximum range; in the present analysis only the rate 
of occurrence of echoes is taken into account. 
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Fic. 4.—Mean monthly diurnal rate curves for 1949-50. 


A meteor producing a trail of 10'* electrons/metre was just detectable with 
this equipment (Lovell and Clegg, 1948). Greenhow and Hawkins (1952) have 
shown that the visual magnitude of a meteor moving with velocity V m/sec, 
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and generating a trail with « electrons/metre is given by the formula 


m, = — 2°5(logy, V3 + logia + logig (3) 


where 7 is the luminous efficiency and p is the number of electrons produced 
by one evaporating meteor atom. For Perseid meteors, with a velocity of 
6 x 10 metres/sec, the equation has been evaluated observationally as 

mM, = 40°1 — 2°5 logy, (4) 
and the limiting magnitude of the equipment with «=10'* electrons/metre 
is +7°5. It is not known how 7 and p vary with the meteor velocity, and hence 
the limiting magnitudes for other velocities is uncertain. 
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Fic. 5.—-Mean monthly diurnal rate curves for 1950-51. 
‘The receiver, transmitter, and other apparatus were checked at frequent 
intervals, and the overall sensitivity was maintained to within +3db. Any 
records which were affected by interference, solar noise or apparatus defects were 
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rejected. ‘There were no appreciable variations of sensitivity due to galactic 
noise. As a final check, measurements were repeated for a second year. 

3. Observational data.—A continuous survey was maintained from 1949 
October 1 to 1951 September 30. In this period 240000 echoes were observed, 
and 5 per cent of these have been attributed to meteor streams (Hawkins and 
Almond, 1952aandb). The few hours of the survey showing echoes from meteor 
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Fic. 6.—Mean diurnal rate curves taken over a period of a year, 


streams have been rejected. The basis of the present analysis is the hourly rate 
of sporadic echoes detected in the sector of each aerial; the mean rate for a given 
hour of Universal ‘lime has been taken over a month. Fig. 4 shows the diurnal 
variation of the hourly rate, plotted month by month, for the year 1949-50; 
Fig. 5 gives the variation for the year 1950-51. In this way the results have been 
subdivided into two surveys which cover two complete rotations of the Earth 
about the Sun. Fig. 6 shows the average of all the monthly curves throughout 
the year, with 6(a) corresponding to Fig. 4 and 6(4) corresponding to Fig. 5. 
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4. Analysis. (a) The concentration of radiants in the plane of the ecliptic. 
‘The monthly curves in Fig. 4 and 5 show three concentrations of radiant points 
at 0130 U.'l. in the winter night-time, at 1030 U.'T. in the summer daytime and 
at o600U.'l, throughout the year. ‘These concentrations are shown more 
clearly in the yearly curves of Fig. 6, where the angle between the sectors has 
displaced the maxima and two intermediate minima by an average vaiue of 2°35 
hours, which, from Fig. 26, corresponds to a mean radiant elevation of 45°-3 
above the horizon. For radiants confined to the ecliptic, crossing a sector at an 
elevation 6, the hourly rate is given by the sensitivity contours of Fig. 3, and varies 
assin#. In this case a mean elevation # would be observed, given by the equation 


0, 6, 
| sinddé 6 (5) 


6, “0, 
where 6, and @, are the minimum and maximum elevations of the ecliptic. For 
Jodrell Bank, 6,= 13° 19’, 6, = 60° 13’, which gives a value of of 41-1. This is 
in good agreement with the observed value, showing that there is a strong con- 
centration of radiant points in the plane of the ecliptic. Between 14" and 22" 
the curves of the sectors are not displaced, indicating a high mean declination. 
‘This implies that when the meteors are observed in the hemisphere containing 
the antapex, there is a more uniform distribution of radiants. ‘The sensitivity 
contours of sector 2 do not pass as close to the celestial pole as those of sector 1 
(Fig. 3), hence the lower value of the evening minimum in sector 2 is further 
evidence for high declination radiants at this time. 

(b) Apparent distribution of radiants in ecliptic longitude.—Yo obtain the mean 
distribution of radiants in ecliptic longitude, time shifts must be applied to allow 
for the angle between the sectors. With radiants at 45-3 elevation, Fig. 2a 
shows that 1°8 hours must be added to the time scale of sector 1 and 0-55 hours 
subtracted from sector 2, so that both sectors are effectively aligned on the prime 
meridian of the observing station. Displacing the yearly diurnal rate curves 
in this way and taking the mean of both sectors, we can obtain the radiant density 
in a band at longitude A, where A is given by the time of transit. ‘This band is 
approximately 6 wide, and its average position is perpendicular to the ecliptic 
with a variation of + 23:5. If A, is the longitude of the Sun, the angle (A —A,) 
may be called the local longitude because the band is in this position at the same 
time each day. Fig. 7 shows the apparent distribution of radiants in local 
longitude, plotted in polar coordinates with the Earth at the origin. It was 
assumed that with equal sensitivities the same number of meteors would have been 
observed in survey | as insurvey II, andthe observed rates in survey | have therefore 
been multiplied by 1°52. 

(c) The distribution of meteor directions in space.—Most radiants, and hence 
most meteor orbits, lie close to the plane of the ecliptic, and we may assume that 
the apparent distribution of orbital directions in this plane is given by Fig. 7. 
‘To obtain the heliocentric distribution of directions at a point on the Earth’s 
orbit we must correct for the Earth’s orbital motion in the following manner. 
If a meteor is moving with a velocity V relative to the Sun, then its apparent 
direction as seen by an observer moving with velocity Vy relative to the Sun, 
may be obtained by reducing the Earth to rest, and applying a velocity — V, to V. 
!f N, radiants/see are observed within a sector de on the ecliptic at an apparent 
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biG. 7.—Polar diagram drawn in thi plane of the Earth's orbit which shows the apparent number 
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elongation ¢ from the apex, then the true elongation e, is given by the equation 


Ve. 
€)=€+8in sine) (6) 


These radiants are located in a sector de which may be found by differentiating 


(6), inverting the differential and substituting for the apparent velocity V’ by the 
cosine rule. 


de V V5 
de, | (7) 


An increase in the apparent velocity V’ increases the number of intercepted 
meteors. If N, is the number of meteors/sec crossing the Earth’s orbit 
clongations between e, and €, + de,, then 


de 


(8) 
N,=N,—. = 
: "it, V24V242V] p COS « 


at 


0 
sy taking V = 40km/sec, 2 km/sec below the parabolic limit, we can correct the 


apparent distribution for the Earth’s motion by using equations (6), (7) and (8). 
The corrected direction distribution is shown in Fig. 8, where V,, is plotted 
againste,. It represents, for any point in the Earth’s orbit, the number of meteors 
which cross at an elongation €, in unit time. 


‘The space density, or number of 
meteor paths per unit volume per radian is equal to N,/V. Since V 


is assumed 
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to be constant, the form of Fig. 8 also gives the space density distribution. This 
analysis is a revision of previous work by the author (Hawkins, 1952) and quoted 
by Lovell (1954). 


ANTAPE x 


SUN 


Fic. 8.—Polar diagram drawn in the plane of the Earth’s orbit which shows the number of meteors 
per unit angle which cross the Earth’s orbit per second. 


(d) The distribution of meteor debris along the Earth’s orlit.—The area under 
the diurnal curve is proportional to the total daily influx of sporadic meteors. 
Let us suppose that the space density of meteors is constant along the entire 
length of the Earth’s orbit. ‘The distribution of apparent radiants in local 
longitude is shown in Fig. 7 and may be written f(A—A,). Because of the con- 
centration towards the ecliptic, the echo rate, N’, at any instant is given by the 
equation 
N'=N.f(A—A,).sin 6. (9) 
The elevation of the ecliptic, 6, changes with the time of day and season of the 
year. Fig. g shows the predicted hourly rates for a sector along the prime 
meridian of Jodrell Bank for each month. When Fig. 4 and 5 are compared with 
Fig. 9 it can be seen that the observed seasonal variation is greater than that 
expected from the gyrations of the ecliptic. It is evident, therefore, that meteors 
are not uniformly distributed along the Earth’s orbit, and the ratio of the areas 
under the observed curves to those under the predicted diurnal curves gives 
the variation in space density from month to month. ‘This ratio is plotted in 
Fig. 10. It shows that meteor orbits congregate in the space traversed by the 
Earth between April and August when the Sun’s longitude is between 40° and 
150. 

5. Discussion. (a) Distortions produced by the ionosphere and other effects.— 
There is no evidence for a difference in meteoric ionization between the night and 
daytime sky. ‘The distribution of radiants in local longitude (Fig. 7) is 
symmetrical; the same number of sporadic meteors are detected by night as 
by day, and this would be expected astronomically. (See section (b) following.) 
Any enhancement of the ionization by the ambient electron density of the FE region, 
or by sunlight, would have permitted the detection of smaller meteor particles 
and would have caused a noticeable day—night asymmetry. 

There is, however, a systematic deviation in the average diurnal curves 
(Fig.6). During the night, between 2000 U.'T’. and o800 U.T.., the rate in sector 1 
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increases and the rate in sector 2 decreases. For radiants near the ecliptic the 
sectors are almost identical, so that the discrepancy cannot be an astronomical 
effect. It represents a systematic change in sensitivity of 1 db between the two 
aerials, 

(b) Significance of the ecliptic distribution.—The concentration of radiants 
towards the ecliptic indicates that meteor orbits are aligned with the plane of the 
Earth’s orbit, and are therefore intimately connected with the solar system. 


OO O06 nrs 


Fic. 9.—-Mean monthly diurnal curves predicted for a sector aligned with the North-South 
meridian. It has been assumed that the space density of meteors along the length of the Earth's orbit 
is constant, and that the distribution of directions is given by Fig. 8. 


It can be seen that the true and apparent longitude distributions (Fig. 7 and 8) 
exhibit a symmetry about the apex. ‘This is a natural consequence of the fact 
that meteor orbits are conic sections, because, if the directions of perihelia are 
uniformly distributed, every night-time orbit has its counterpart in the daytime 
sky. Night-time meteors are approaching a perihelion point inside the Earth’s 
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orbit, and daytime meteors are receding. ‘The probability of meteors colliding 
destructively with Mercury, Venus, the Earth and particles in the zodiacal cloud 
is negligible for one perihelion passage. Hence we would expect, as shown by 
this survey, the same number of meteors incident from two hemispheres with 
poles directed towards and away from the Sun. 

‘The maximum in the radiant distribution at the apex is caused by meteors 
which are swept up by the Earth. ‘The side lobes, at elongations of 66° from the 
apex, close to the Sun and antihelion points give the most frequently observed 
sporadic meteor orbits if the mean velocity is known. We may assume that the 
heliocentric velocity is between 34 and 40km/sec (Lovell 1954, page 234), and 
we may compute the series of orbits that correspond to an elongation of the radiant 
of 66°. ‘These are given in ‘lable I where it can be seen that the most frequently 
observed orbit is probably elliptical, with high eccenticity and with a semi-major 
axis comparable to the short period meteor streams. 

(c) The distribution of meteor directions in space.—When the correction for the 
Earth’s velocity is applied the results show a remarkable preponderance of direct 
moving meteors; orbits to the left of the line between the Sun and antihelion in 
Fig. 8 are retrograde. Points near the longitude of the antapex have been rejected 
because they do not represent orbits in the plane of the ecliptic. ‘The distribution 
could imply that meteors move in almost circular orbits overtaking the Earth. 
If this is the case, a heliocentric velocity of 30 to 35 km/sec should be assumed in 
section 4(c) instead of 40km/sec, and this would further increase the number of 
directly moving meteors in space. 


However, there are uncertainties in this part of the analysis that must be 
borne in mind. Large corrections, amounting to a factor of 42 for antapex 
meteors, have been applied, and the correction is critically dependent on the 
assumed heliocentric velocity. Furthermore, the distributions are for a limiting 
electron density of 10!*/metre in the meteor trail. The corresponding visual 
magnitude and mass of the meteor are undoubtedly functions of geocentric 
velocity, and hence of €,, so that the distribution when taken to a certain limiting 
mass would be different from that in Fig. 8. 


TABLE 


The orbits that correspond to assumed heliocentric velocities between 
34 and 40 km/sec 


V a e q A 
km/sec Semi-major Eccentricity Perihelion Aphelion 
axis distance distance 


(d) Irregularities along the Earth’s orbit.—-Fig. 10 shows the seasonal variation 
in the total meteor flux incident on the Earth, which increases by a factor of 4 
between February and June. ‘The fact that the maximum corresponds to the 
summer solstice might indicate that the measurements had been influenced by 
variations in the height of the meteor layer, in ionization levels in the E region 
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or by other effects dependent on solar radiation. However the annual variation 
of height (Jacchia) is small, and the ionosphere does not produce diurnal distor- 
tions (Section § a) so that seasonal etfects would not be expected. 
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Variation in the space density of meteors alone the leneth of the Earth's orbit (top fisure) 
compared % ith the cometary index (bottom fiwure) 


The variations in the density correlate significantly with the cometary index 
of Herschel and Hoffmeister (1948) which is plotted for comparison in Fig. 10. 
This index is proportional to the number of comets crossing the Earth’s orbit 
and their closeness of approach. ‘lhe correlation suggests that sporadic meteors 
and comets have a similar origin. 
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A RADIO ECHO METHOD OF METEOR ORBIT DETERMINATION 


Jj. C. Gill and F. G. Davies 
(Communicated by A. C. B. Lovell) 


(Received 1955 October 29) 


Summary 


A radio echo system using three receivers spaced on the ground by 4 km 
has been used to determine the orbits of individual meteors. About 200 
orbits can be calculated from the meteor echoes photographed during a 24 
hour period of observation. An analysis of Geminid meteors observed in 1954 
indicates that the accuracy obtainable is limited to about + 2 km/sec in velocity 
by atmospheric deceleration, and to + 3° in radiant position by turbulent 
winds in the meteor region. 


1. Introduction.—Radio observations of meteors have been made since 1946. 
Velocity measurements based on the observation of diffraction patterns have been 
made both by pulsed (1) and continuous wave (2,6) technique. ‘lhese have 
yielded much information concerning both shower and sporadic meteors. In 
the case of showers, radio techniques are also available for delineation of radiant 
positions (3) and this has led to the discovery of daytime streams and the deter- 
mination of orbits. In the case of sporadic meteors no measurement of individual 
radiants has been possible but, although information concerning their orbits 
has been of a statistical nature, it has been established that a large proportion 
move in orbits of short period, and that less than 1 per cent have an interstellar 
origin. For the more detailed analysis of the orbits of sporadic meteors, and the 
examination of radiant structure of showers, it is desirable to be able to measure 
both velocity and direction of motion for individual meteors. Following 
a suggestion by Kaiser, the pulse technique of velocity measurement (1) has 
recently been extended to accomplish this by the use of three spaced receivers, 

2. Theory.—\f two observing stations (A and B in Fig. 1) each with transmitter 
and receiver are used, there will be a time difference between reception of 
corresponding points on the diffraction patterns at the two stations equal to the 
time taken for the meteor to travel the distance between the specular reflecting 
points P,, P,. ‘The velocity V of the meteor can be determined from either 
or both of the patterns, hence the length P,P, can be found if the time delay ¢ 
is measured. ‘The value of cos@, where @ is the angle between AB and a line 
drawn through B parallel to the trail YX, can therefore be determined, since 


and AB can be measured directly. 

In the general case when XY and AB are not coplanar, cos@ becomes the 
direction cosine between the line of flight XY of the meteor and the base-line AB. 
The addition of a third station enables a second direction cosine to be found, 
thus defining the flight direction and hence the radiant of the meteor. ‘The 


_ Ve 
AB AB 


106 F.C. Gill and F. G. Davies Vol. 116 


velocity is already known and, given the time of occurrence, the orbit can be 
computed. 
‘The error in the measured angle 6 is given by 


2712 
64 = cos 6 | (7 + (*)'| 
V t 


where 4V and 4t are the errors of measurement in V and ¢ respectively. ‘The 
time delay is directly proportional to the base-line AB and thus the angular error 
will be reduced by an increase in base-line until this is such that the fractional 
error in ¢ is small compared with the fractional error in V. For a typical case 
the two error terms are equal with a base-line of the order of 1-5 km. 


ic. 1.—Radio stations at A and B observe the meteor trail at the specular reflection points 
P4,Pp. From the ratio of P4 Pp to the bescline AB, the radiant direction can be determined in 
the ti 


o-dimensional case 


In practice it is not necessary to use complete transmitting and receiving 
equipment at all three stations. Ifa transmitter at A is used with a remote receiver 
at B, the point of specular reflection lies at the mid-point of P, P;,, provided that 
the range of the echo greatly exceeds AB. It can be shown that under these 
conditions the diffraction pattern received at B is for practical purposes the same 
as would be obtained at a complete transmitting and receiving station situated 
midway between A and B. ‘The delay between patterns received at A and B 
is therefore half that for the case of two separate transmitting and receiving 
stations. 

A single transmitting and receiving station, with two remote receivers, is 
thus sufficient to provide the three diffraction patterns necessary for orbit 
determination. 

An equipment using a single transmitter and capable of recording the three 
diffraction patterns necessary for orbit determination has been constructed 
and was first used at Jodrell Bank in 1953 December. 

3. Apparatus.—The basic equipment parameters are given below: 

Radio frequency: 36°3 Mc/s. 

Transmitter: 600 pulses per second are radiated; a pulse length of 
15 microseconds and pulse power of 100 kW are usually employed. 

Aerials: ‘The transmitting aerial most frequently employed has an absolute 
power gain of 28, with a beam width of + 12°. ‘The aerial is directed slightly 
north of east at an elevation of 30°. Reception at all stations is by aerials having 
absolute power gains of 12 in the direction of the main beam of the transmitting 
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aerial; the latter determines the effective beam width of the system as a 
whole. 

‘The transmitter and one receiver are located at Jodrell Bank ; remote receivers 
are situated approximately south and east of this at distances of 3-5 kilometres. 
Signals received at the remote stations are amplified, converted to link frequencies 
in the 75-80 Mc/s band, and after further amplification are re-transmitted to the 
home station, where the signals received at all three stations are displayed on 
cathode-ray tubes, a fourth tube being used to record the range of the echo at 
one station. 

The received pulses are recorded photographically using a moving film 
system triggered by reception of the echo itself. ‘he echo pulses are recorded 
for a maximum time of o°5 seconds, to ensure satisfactory reception of the patterns 
from all three stations for the slowest meteors. 

The time of occurrence of the echo is also recorded. A specimen record 
obtained with this apparatus is shown in Plate 2. 

With this equipment the majority of meteors for which orbits are determined 
lie in the magnitude range +6 to +8; in an average period of 24 hours some 
200 sporadic meteors giving patterns satisfactory for orbit determination can be 
recorded. 

4. Reduction of results and computation of the orbit.—Vhe derivation of the 
orbit from the recorded information falls into two stages: the determination 
from the film records of the meteor velocity and the time delays between stations, 
and the computation of the orbit from these results. 

The first stage involves fitting the best theoretical diffraction pattern to that 
observed; the velocity corresponding to this pattern and its separation from 
the other two are recorded. ‘This is now achieved using an analogue computer, 
with considerable saving in time over alternative methods of computation using 
a desk machine. 

The computation of the meteor orbit from the observed velocity and delays 
follows the method described by Portes (4), including corrections for atmospheri 
deceleration, diurnal aberration and zenith attraction. Since the calculation 
is alengthy process and several hundred orbits may be computed from observations 
over a single 24-hour period, it was found necessary to use electronic computation. 
The Manchester University digital computer is used and is arranged to print 
the following information for each meteor: 


Corrected geocentric radiant position in altazimuth, equatorial and 
ecliptic co-ordinates. 

Elongation of radiant from apex, and heliocentric velocity. 

Heliocentric radiant in ecliptic coordinates. 

Longitude and argument of perihelion. 

Longitude of ascending node. 

Inclination of orbital plane. 

Semi-major axis, a. 

Eccentricity, 

\phe lion and pe rihelion di tance 

‘The parameter of the orbit p=a(1 —e?). 


Ecliptical coordinates of aphelion direction and pole of the orbit. 


his process takes approximat« ly 40 SCC onds tor each meteor 
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5. The geminid shower, 1954.—In order to determine the accuracy of the method 
an investigation of the Geminid shower of 1954 December has been made. 
Geminid meteors have previously been extensively observed both by radio and 
photographic techniques. ‘The results of the most important measurements 
are summarized in ‘lable I. 


TABLE I 


Observed 


Observing station Radiant (Sun’s longitude = 260°6’) Velocity spread in 
R.A, Dec, Diameter km/sec. velocity 
, km/sec 

Jodrell Bank 1947 34°4 +1°2 


(Radio) (5) 1949-50 112°3 


35°9 +4°6 


Ottawa 
(Radio) (6) 


35°2 


Harvard 
(Photographic) (7, 8) 112°5 + 32°8 0°22 36°4 +o'2 


In ‘lable I the radiant positions have been corrected for the effects of diurnal 
aberration and zenith attraction; the velocity quoted is the mean observed 
velocity in the case of the radio determinations, and the mean velocity at the top 
of the atmosphere for the photographic results, after removal of the component 
due to the Earth’s rotation. ‘The figure given for the spread in velocity is the 
standard deviation of an individual observation, assuming the distribution to 
be Gaussian. In the case of the radio data the spreads are considerably larger 
than the random errors of measuring the velocity from the diffraction patterns. 
‘The radio results refer to + 6th magnitude and fainter meteors; the Harvard 
data are for meteors brighter than zero magnitude. 

6. Radiants.—Sixty-one Geminid meteors were recorded on 1954 December 
13. ‘lhe radiants of these are shown in Fig. 2 plotted separately for three 
consecutive periods. [From considerations of the accuracy of measurement of 
velocity and time displacement of diffraction patterns errors of + 0-5 km/sec. 
in velocity and 1-5 in radiant position would be expected. It is seen that, 
although the mean radiant agrees with previous photographic and radio results, 
the spread of individual radiants exceeds that anticipated from the errors in 
reading the diffraction patterns, and also the observed shape of the radiant changes 
with its position relative to the base-lines used, showing that the spread is not 
astronomical in origin. 

It has been shown by Kaiser (9) that such a radiant spread can be caused 
by a uniform wind shear in the region of the trail. ‘The chief effects of such 
a shear are found to be :— 


(a) The diffraction pattern observed shows no distortion unless very high or 
second-order shears are present. 


(6) The value of velocity obtained is slightly too high; the error, however, 
is negligible compared with the random errors of determination, for uniform 
shears. 
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(c) The displacements between diffraction patterns are seriously affected. 
‘The fractional error in delay is found to be 
2k du 
where b= 
2-0 V’ ds 


R=echo range, 


= 
| 


= meteor velocity, 


wind velocity parallel to line of sight from station to trail, 
s = distance along trail. 


= 


e 


Li 
105 


F'1G.2.—Kadiant points of Geminid meteors observed on 1954 December 14. Solid points 
represent undistorted diffraction patterns. Open circles represent diffraction patterns showing 
evidence of turbulent winds. The straight lines E and S re present the errors introduced by 10 per 
cent errors in estimating the time displacements of the eastern and southern remote station re spectively 
from the home station. 


(a) Sidereal time 07" 30™ to ogh yom 
(b) Sidereal time og 40m to 10h 30m 


(c) Sidereal time 10 30™ to 124 oom 
Abscissae : R.A. in degrees. 
Ordinates : declination in degrees, 


The effects on radiant position produced by 10 per cent errors in the time 
displacements between diffraction patterns are indicated for the central time of 
each period in Fig. 2. If the radiant is so placed that there is no time displacement 
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between the diffraction patterns observed at a pair of stations, then wind shear 
cannot affect the magnitude of this displacement, and this accounts for the very 
small S vector in Fig. 2(c). In this figure the radiants lie along a line parallel 
to the E vector with a spread of +1. ‘The observed spread in radiants therefore 
does not conflict with Whipple’s measurement of radiant diameter. “lhe Geminid 
observations show that the radiant spread is equivalent to errors in displacements 
of about 5 per cent. Such errors would be caused by wind shears of 
1omsec'km!. Shears of this order are commonly observed by other methods, 
e.g. by photography of persistent meteor trails. ‘The accuracy of radiant 
determination for individual meteors is therefore limited by this effect to +3 
for an average case. 

7. Observed velocities.—Vhe observed velocity will be affected by errors of 
measurement and by atmospheric deceleration. Photographic velocity deter- 
minations by Whipple (10) indicate that for 3rd magnitude meteors at their 
maximum of ionization the loss in velocity due to atmospheric resistance 1s 
about 0-5 km/sec for meteors having an initial velocity of 35km/sec. As the 
height of the meteor relative to its maximum of ionization is not known in this 
experiment, there 1s an uncertainty, of the same order as the loss in velocity at 
maximum ionization, in the correction to be applied and hence in the geocentric 
velocity of the meteor. ‘The theory of Herlofson (11) predicts that the loss in 
velocity of a spherical meteor travelling through the atmosphere ts 

~ In (1 ~ fe") 
where L = latent heat of evaporation of meteor, 
\ = heat transfer coefficient, 
h height of meteor above maximum of ionisation, measured in units 
of atmospheric scale height, 


veloc ity of meteor on entering atmosphere. 


The loss in velocity, following this theory, is independent of the meteor 
magnitude, and therefore the necessary correction should remain o-5 km/sec 
for the fainter Geminid meteors. 

The velocities of meteors from the Geminid radiant recorded during the 
period 1954 December 10-14 are shown in Fig. 3. It is apparent that these 
radio determinations are systematically lower and exhibit a larger spread than the 
photographic results, the mean observed radio velocity being only 34:1 km/sec, 
After applying the correction for the Earth’s rotation, this becomes 34-3 km/sec 
compared with 36:4km/sec found photographically. ‘(he spread of the radio 
observations is much greater than the error of reading the diffraction patterns, 
and exceeds by a factor of 4 the spread due to atmospheric deceleration predicted 
from Herlofson’s theory. 

The radio observations of velocity quoted in ‘Table | have slightly higher 
values than those given above, although they are significantly lower than the 
photographic velocity. A large spread in velocity occurs for all the radio 
determinations. ‘lhe difference between the 34°3km/sec found here and the 
previous mean radio velocities probably arises from the inclusion of sporadic 
meteors in the earlier data. ‘These would be likely to have higher velocities than 
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Geminids, and are excluded from the present results for which radiants have 
been determined. 


34 36 


Histogram of velocities of Geminid meteors recorded het 


en 1954 December to and 14 
lhscissac elocits in km see 


Ordinates : number of meteors 


The difference between the radio and photographic observations can only 
be due to the difference in magnitude range covered. Either the true geocentric 
velocities of the fainter meteors are lower and exhibit a large scatter or they suffer 
much greater retardation in the Earth’s atmosphere than do the brighter meteors, 
Since no systematic shift in radiant position with measured velocity is observed, 
it is unlikely that the cause of the lower velocities can be astronomical in nature. 
In the case of atmospheric retardation, the lowest observed velocities will be for 
meteors seen at the ends of their paths in the atmosphere, and consequently at 
the lowest observed heights. A marked correlation between height and velocity 
should be observed. If the velocity spread is astronomical, the only correlation 
between meteor velocity and height will be that given by the relation 


COS xX 


Ve 


Pmax 


where Pax = atmospheric pressure at maximum of ionization, 
y =zenith distance of radiant, 


V = meteor velocity. 


In Fig. 4 the measured velocities are plotted as a function of height. ‘The latter 
has been obtained by an approximate method based on a knowledge of the radiant 
and range of the meteors. ‘lhe errors in the plotted heights are of the order of 
+2km. It is seen that the correlation between velocity and height is marked, 


and that the slope of the height-velocity curve is much greater than that expected 
if the observed velocity distribution were astronomical in origin. It must therefore 
be concluded that the spread in observed velocity is due to deceleration in the 
Earth’s atmosphere, and that there is no reason for believing that the velocity 
at the top of the atmosphere of the faint Geminid meteors differs significantly 
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from that observed photographically for the brighter members of the stream. 
Assuming Whipple’s value of the velocity to apply to the faint meteors, the median 
value of the correction for atmospheric deceleration is found to be 2:0 km/sec. 
There is, therefore, an uncertainty of this order in the geocentric velocity of 
meteors recorded with velocities around 35 km/sec. 


30 


hic. 4.--Observed velocity of Geminid meteors as a function of height of point of observation. 
The dashed line indicates the relation between height of maximum ionization and velocity for 


sporadic meteors. 


Abscissae : velocity in km/sec. 
Ordinates : height in km. 


8. Conclusions. In the light of these results, it appears that the accuracy 
obtainable using this technique is limited, for 35 km/sec meteors, to + 2 km/sec 
in velocity and +3 in radiant position by atmospheric effects. ‘These effects 


are expected to decrease for meteors of greater velocity. 

This accuracy is amply sufficient for surveys of sporadic activity, but for more 
detailed work, such as investigations of the radiant structure of showers, greater 
accuracy is desirable. Corrections for the effects of wind shear on radiant 


determination could be derived from measurements of the wind components 


along the lines from the observing stations to the points of specular reflection 


on the trail. ‘The error in geocentric velocity due to atmospheric retardation 
could be reduced if the velocity at a number of points spaced along the whole 
trail were measured, ‘his would require a series of receiving stations spaced 


along each baseline, and would allow the velocity of the meteor before entering 
the atmosphere to be found by extrapolation. A simpler, though less accurate, 
method of correction would be to determine the height of the trail, and use the 


relation between height and retardation demonstrated in Fig. 4 to derive the 


necessary correction. 

The apparatus has been used since 1954 May for a systematic survey of 
sporadic and summer daytime activity, the results of which will be published 
later. It is intended to continue and extend these observations in the 


future. 
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DISSIPATION OF MAGNETIC 
IN AN IONIZED GAS 


THE ENERGY 


T. G. Cowling 
(Received 1955 December 2) 
Summary 


A general formula is derived for the rate of dissipation of magnetic energy 
in an ionized gas. ‘The formula is used to confirm Piddington’s value for the 
conductivity effective in producing dissipation. It is shown that, contrary to 
Piddington’s assertions, an effective electric force must be taken into account 
when considering currents in the presence of a pressure gradient; however, 
this effective electric force has little effect on the dissipation. Piddington’s 
results are completed by a general discussion of conductivity in a partially 
ionized gas. The results are applied to dissipation in interstellar gas; the 
effect of the magnetic field on the conductivity is shown to increase the rate of 
dissipation up to a point where the dissipation is not insignificant. 


1. In recent papers J. H. Piddington (1) has studied certain problems of 
the dissipation of magnetic energy when an ionized gas moves in a magnetic 
field H. He starts from the equation 


x (1) 


where j is the current-density, E,’ and E,’ are respectively the components 
parallel and perpendicular to the magnetic field of the electric force on the 
moving material, and o9, 0, and o, are appropriate conductivities. He finds 
that the rate of dissipation of magnetic energy per unit volume when currents 
flow perpendicular to the magnetic field is j?/o,, where 


as 2 
0, > 90; T Oe 


‘That is, o, is the conductivity effective in dissipation in this case. 

I do not query Piddington’s result; I had, in fact, obtained an equivalent 
result a month or two before reading his paper. But his derivation appears 
unnecessarily complicated ; | give below a simpler and more general derivation. 
Moreover, he accompanies his discussion by remarks suggesting that the 
‘equivalent electric fields’’ introduced by myself (2) and Schliiter (3) are 
concepts which correspond to no essential physical reality, the true currents 
being always given by (1). Such a suggestion is incorrect; the relative diffusion 
of ions and electrons due to forces other than electric forces produces electric 


’ 


currents, and the equivalent electric force was introduced simply in order to 
take account of these currents. Schliiter, by skilful manipulation of such 
equivalent electric forces, was able to show that in certain cases the conductivity 
affecting the dissipation of electric energy was o9, not o,. Further investigation 
of the equivalent electric forces accordingly appears desirable, in order to be sure 
that they do not invalidate Piddington’s results. 


(2) 
2 
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Piddington found that in a lightly ionized gas, when positive ions, as well 
as electrons, are free to spiral round the lines of magnetic force between collisions, 
a; is much smaller than og (which is the conductivity in the absence of a magnetic 
field). ‘This means that the magnetic field may increase the dissipation rate very 
appreciably in, for example, H1 clouds in interstellar space and in the chromo- 
sphere above a sunspot. For the exact calculation of this important effect a 
further analysis of current-flow in a partially ionized gas is needed; this is also 
given below. 

2. Simple theory of dissipation.—Assume first the validity of equation (1). 
The electric field E’ of which Ej, and E/ are components is given in terms of 
the field E on matter at rest and the velocity v by 


E’=E+vxH (3) 
(e.m.u. are used throughout). ‘The velocity v is taken to be the mass-velocity, 
defined as such that the total momentum per unit volume is the same as if every 
particle present had velocity v. 

The total magnetic energy is given by 


, _ 
Wy =| go dr (4) 


integrated over the whole of space. Since 0H/dt= —curl E, the rate of change 
of the magnetic energy is 


| \div(E x H)+ E.curlH} dr. 
7 


In this, by Green’s theorem, the divergence term gives no contribution on 
integration over the whole of space, if one assumes a field which tends to zero 
sufficiently rapidly at infinity. ‘Thus, since curlH = 47j, 


dW, 


i | E .j dr. 


The changes in Wy are not wholly due to dissipation; the mechanical force 
j x H of electromagnetic origin does work on the material at a rate (j xH).v 
per unit volume, and this is done at the expense of the magnetic energy. Hence 
dW ,,/dt includes a part —(j x H).v per unit volume due to the work done; the 
part of dW,,/dt due to dissipation is 


— | E.jdr+ | (jx H).vdr-= | (E+vxH).jdr 


- | E’ .jdr. 


That is, the dissipation of energy due to electrical resistance is E’.j per unit 
volume and time. 


(5) 
4a 
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Kor a current parallel to the magnetic field j=0,€’, so that the rate of 
dissipation is j*/o, per unit volume—the usual result. For a current perpendicular 
to the field E’ reduces to E', and 


E’.j=0,F;?. 


But, according to equation (1), j is now the sum of components 0,E), o,E' at 
right angles. ‘Thus 


1 


=(0," + 


and the rate of dissipation per unit volume is 


This is Piddington’s result. 

3. Equivalent electric force.—In his reference to equivalent electric forces 
Piddington has confused the influence of a pressure gradient in causing diffusion 
and its effect in producing a state of motion, which reacts on the current equation 
through the term v x H in E’. Clearly only the diffusion effect can be treated as 
an equivalent electric force; the effect on the state of motion is already adequately 
allowed for by the term v x H in E’. Equally clearly equation (1) is incomplete 
until the diffusion effect is taken into account. 

For a wholly ionized gas the diffusion effect is equivalent to adding to E’ in 
equation (1) the electric force E” given by 


E” = (grad p,)/ne (6) 
where p, is the electron-pressure, n the number-density of electrons, and e the 
electronic charge (2, p. 468). A simple derivation of (6) will be given for a binary 
mixture of positive ions of one type and electrons. From this it will be clear 
that the introduction of E” has no direct relation to Piddington’s discussion of 
steady motion. 

‘The gas is supposed to consist of ions of mass m, and partial pressure p,, and 
electrons of mass m, and partial pressure p,. ‘The number-density of each is n; 
thus if the ion and electron gases are supposed to have the same temperature, 
pi=p,.. Relative diffusion of the two gases occurs if the forces acting on them 
would, in the absence of diffusion and of collisions between ions and electrons, 
give them different accelerations f, and f,; the velocity of diffusion is then 
proportional in magnitude to f,—f,. If the only applied forces are gravitational 
and electromagnetic, 


nm f,=nmg +neE’ — grad p, | 
nm. f,=nm,.g—neE’ — grad p, (7) 


It is to be noted that the partial pressures p,, p, alone enter into these equations ; 
p,, being due to the transport of momentum by electrons, affects ions only through 
collisions, which are neglected in deriving the equations. 

Equations (7) give 


e I I 
f —f + —)€’- rad p, + —— grad p,. 
m, mM, nm,” nm,” Pe 
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Since m,/m, is negligibly small, this equation approximates to 


I 


= e(E’ + E”)/m, 


by (6). The velocity of diffusion, and the resulting electric currents, are 
proportional in magnitude to f,—f,, and so to E’+ E”. ‘Thus, in calculating the 
electric currents, E’ in equation (1) has to be increased by the equivalent electric 
force E”. 

4. Dissipation and €".—When E” is added in, equation (1) is replaced by 


+ Ef) + 0,(E; + + +E). 


An argument similar to that at the end of Section 2 now shows that if E’ + E” 
is perpendicular to the magnetic field 


(E' + E").j 
‘Thus the rate of dissipation per unit volume is given by 
E’ .j =j?/o,— E” .j. (8) 


Similarly if E’+ E” is parallel to the magnetic field 


E’ .j =j?/0,— .j. (8’) 


In either case there is a new term — E” .j representing the effect of the equivalent 
electric force. 


‘The question is whether this term increases or diminishes the total dissipation. 
‘There is no universally valid reason why it should do either; it does so only if 
E” is related in some way to j. And it is easy to see in the cases of most interest 
that E”.j is, on an average, equal to zero. 

Consider, say, a steady discharge. ‘lhe current flowing in the discharge 
produces a magnetic field constricting the discharge, as a consequence of which 
the pressure inside the discharge column exceeds that outside. ‘This gives rise 
to a pressure gradient, and an E”, which are perpendicular to the current. ‘Thus 
E”.j=o and (the magnetic field being perpendicular to j) the dissipation is 
j*/oy, just as if E” were not acting. ‘his agrees with Schliiter’s results (3), since 
in a wholly ionized gas o, is (roughly) equal to a», the conductivity in the absence 
of a magnetic field. 

A case of greater astronomical interest is that in which the lines of force are 
nearly frozen into the material. ‘lhe material motions bend the lines of force 
and the magnetic field reacts back on the motion; the currents flowing are not 
essentially determined by the conductivity, but by the changes in magnetic field, 
in accordance with the equation curlH = 47j. If p is the density, the equation 
of motion of the gas as a whole is 


dv 
Pa = ~gradp +jxH. (9) 


The field reacts on the motion through the force j x H; in part this force helps 
to accelerate the gas, but in part it tends at any time to be balanced by the pressure 
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gradient. ‘hus the part of j which is related to the pressure gradient can be 
expected to be perpendicular to it, and so to the equivalent electric force E”. 
‘his means that, although the new term — E” .j in the dissipation may fluctuate, 
its average value is likely to be zero. 

Of course, as Schliiter and Biermann (4) have pointed out, the equivalent 
electric force E” may be able to build up a small magnetic field from nothing 
(so that — E”.j corresponds to a negative dissipation). When H is very small, 
the connection of j with E” through the term j x H in (g) becomes unimportant. 
On the other hand, the equation of conduction now becomes j=o,(E’ + E”), 
and a connection between j and E” is possible through this equation if E” is not 
balanced by an electrostatic field. However, this effect is important only when 
the magnetic field is very small. 

5. Partially ionized gas.—'Yo study electric currents in a partially ionized 
yas we must consider a mixture of at least three constituents—electrons, positive 
ions, neutral atoms. ‘The general formulae governing diffusion in a multi- 
component gas in a magnetic field are excessively involved. However, relatively 
simple equations can be obtained in a first approximation. 

Let the constituent particles of the rth gas have mass m,, charge e, and number- 
density n,; let p, and p, be the partial pressure and density of the rth gas, and 
V, be its velocity of diffusion relative to the mass. ‘Then the velocities of diffusion 
are to be calculated from the equations * 


dv 


Pr dt =p,g—gradp, +n,e(E’ +V, x H) — (10) 


together with 


up,V, = 0, (11) 
Since in this case j= Ln,e,V,, and &n,e,=o, when the equation (10) is summed 
over all values of r the equation of motion (g) is obtained. 

Equation (7) can be interpreted as follows (3). ‘The velocities of diffusion 
are assumed to be small compared with the velocity v of the mass; thus p,dv/dt 
gives the mass-acceleration of the rth gas per unit volume. Equation (7) equates 
this to the sum of the forces on the gas per unit volume — gravity p,g, the gradient 
of its partial pressure, the electric force on particles moving with average 
velocity v + V,, and finally drag forces due to the diffusion relative to the other 
gases of the mixture. ‘The drag due to diffusion relative to the sth gas is 
proportional to the relative velocity V,—V,, and to the number of collisions 
between particles of the rth and sth gases per unit volume and time. 

Here it is sufficient to consider a gas composed of electrons, of charge —e 
and mass m,, positive ions, of charge e and mass m, and neutral atoms, of mass m,. 
‘The electrons and ions spiral round the lines of magnetic force with spiral- 
frequencies w,, w,, where 

= eH/m.,. (12) 
‘The number-densities of electrons and positive ions are each taken to equal n; 
the neutral atoms are supposed to contribute a fraction f to the total pressure p, 
and a fraction F to the total density. 

* Such equations can be obtained, e.g., from (2), equation (22), on multiplying by m,C, and 


integrating, using equation (19) of (2) and evaluating the collision-integrals to a first approximation 
similar to that used in that paper. 


4 
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The diffusion velocities of ions and electrons relative to the mass are taken 

to be V,, V,+V. Since the mass-density of the electron-gas is negligible, 

equation (11) shows that the velocity of neutral atoms relative to the mass 1s 

~(1—FYV,/F; the velocities of ions and electrons relative to the neutral atoms 
are respectively V,/F and V+V,/F. ‘The current-density is given by 

j= —neV. (13) 


By analogy with this we write 
j, = neV,. (14) 


Thus j, denotes the current-density due to the diffusion velocity of ions (cancelled 
by an equal current-density due to the component V, in the diffusion velocity of 
electrons). 

If the mass-density of electrons is neglected, the electron equation of the 
set (10) takes the form 

o= — grad (15) 

In this, —6,,V denotes the drag on the electron-gas due to collisions with ions. 
Conventionally a collision can be regarded as destroying, on an average, the 
relative motion of the colliding particles. Because of the smallness of m,, this 
means effectively that each collision of electrons with ions results in a loss m,V 
of momentum to the electrons. ‘hus if 7 is the mean time between successive 
collisions of an electron with ions, 


(16) 
Similarly if 7, is the mean time between successive collisions of an electron with 
neutral atoms, 
(17) 
On use of these expressions for 4,, and @,,, equation (15) becomes 
o= —grad 


If j and j, are introduced from equations (13) and (14) in place of V and V,, this 
becomes 


o= —gradp, —neE’ +(j—j,) (18) 
where 
K (w,7) A K, (w,7,) - (19) 
The quantities «, «, are pure numbers; the condition for electrons to be able 
to spiral between collisions is that « +«, shall be small. 
In the equation for ions in the set (10) one may similarly put 


nm 7, (20) 


However, 7, is not exactly the mean time between collisions of an ion with neutral 
atoms. Each collision destroys only a fraction m,/(m,+m,) of its momentum 
mN,/F relative to the neutral atoms; thus 7, is m,/(m, + m,) times the collision- 
interval. On use of (17) and (20), the equation for ions becomes 


(1— (1 — — grad p, + ne(E’ + V, x AH) + —nmz,' FV, 
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Once again V and V, are replaced by j and j,. Then if 
K, (w,7;) (21) 


the equation becomes 
(1—F )p di — F)pg — grad p, + neE’ + j, x H—«Hj—«, (22) 


‘The equation for neutral atoms will not be used explicitly. Its place will 
be taken by (g), the equation of motion for the gas as a whole. 

6. The current equation.—'Yo obtain an equation giving j, we must eliminate j, 
between equations (18) and (22). Add the equations; then 


(1—F)p (1 — "Hj, (23) 


Substitute for dv/dt from equation (9); then, since (p,+p,)=(1 —f)p, one finds 
that 


(x, +«,)FHj, = Fj x —(F —f) grad p +«, Hj. (24) 


If this equation is used to eliminate j, from equation (18), the desired equation 
for j is obtained; it is 


KK, + KK, + KK 


nek’ ‘“Hj+(1—2BF)j gradp 


{F(j 


(F—f)gradp x H} (25) 


where 


(26) 


‘The case of a wholly ionized gas is obtained by making «,, «,, f and F tend 
to zero, their ratios remaining finite; this gives 


neE’ + gradp=xHj+jxH. (25’) 


Since in this case p,=}p, the term }gradp is responsible for the effective 
electric force E”. ‘The term j x H is responsible for the effect of the magnetic 
field on the conductivity (in spite of its apparent form, «/ is actually independent 
of the magnetic field). If the terms j x H and gradp in the equation of motion 
are comparable (which is frequently the case in astronomical applications) the 
term | gradp must not be neglected in any problem in which the effect of the 
magnetic field on the conductivity has to be taken into account. ‘This is the 
reason for insisting earlier on the importance of the equivalent electric force, 
in spite of its normally unimportant effect on the dissipation. 

For a partially ionized gas, any equivalent electric force would have to 
include two terms proportional respectively to gradp and H ~ grad p. 
However, because of the complication of the expression for such an equivalent 
electric force, there is little point in introducing it. ‘The terms involving grad p 
and H = grad p in equation (25) are in general comparable with those of the 
other terms which involve the magnetic field unless the degree of ionization is 
small (f and F both nearly equal to unity). 


{ 
k YK 
i , 
| 
a 
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7. Dissipation in a partially ionized gas.—Equation (25) can, if desired, be 
solved to give an equation like (1), with extra terms involving the pressure 
gradient. However, the dissipation E’.j per unit volume and time can readily 
be calculated direct from equation (25). 

‘The term involving j x H in (25) makes no contribution to E’.j. ‘The term 
proportional to grad fp, like that involving E” in Section 4, can be expected to 
make no contribution on an average, since grad /p, in so far as it is affected by 
the electric currents, depends on them only through a component proportional 
to jx H. Attention will therefore be confined to the other terms; the 


give 
a rate of dissipation per unit volume equal to 


F(F—f)H"' grad p .(jxH)}. (27) 
lor currents flowing parallel to the magnetic field, (27) reduces to j*/a,, where 


ne(K, +K,) 


H(kx, +KK, + (28) 


As can be seen from equation (25), a) (which is independent of 7) denotes the 
conductivity parallel to the magnetic field, as in equation (1). 

lor currents flowing perpendicular to the magnetic field, the term involving 
grad p .(j x A) in (27) can be expected in general to cancel an appreciable part 
of the term involving |j x H|* when Ff, since grad p can be expected to have 
an appreciable part proportional to j x H. (The condition / ~/ implies that the 
ion-mass m, is less than twice the neutral-atom mass m,.) However, for a gas 
only slightly ionized F and f are both nearly equal to unity, and the term involving 
grad p.(j x A) is relatively unimportant. 

If this term is neglected, the dissipation (27) becomes /*/o, for currents 
perpendicular to the magnetic field, where 


Oy fF? 


2¢ 
Oz KK + KK (29) 


‘he expression on the right of this equation gives the ratio in which the dissipation 
is increased by the magnetic field. If the neutral atoms are an appreciable fraction 
of the whole (F not small) and both electrons and ions spiral freely between 
collisions («, «, and « all small) this ratio may be very large. ‘The increased 
dissipation is due to collisions between neutral and charged particles, which 
tend to destroy the common velocity of diffusion V, of the charged particles 
relative to the neutrals (the velocity of ambipolar diffusion, as Schliiter (3) calls 
it). Ascan be seen from equation (24), the “ current’”’ j, due to V, is large compared 
with j when electrons and ions are free to spiral (roughly in the ratio F?/(«, +«,)) ; 
thus the corresponding dissipation can be expected to be large. In cases con- 
sidered by Schliiter, the ‘‘current”’ j, due to ambipolar diffusion was supposed 
not to flow; thus he did not find the increased dissipation given by equation (29). 
In the cases of interest in astronomy, however, there is nothing to prevent j, from 
flowing, and the increased dissipation is certain. 

Piddington obtained a special case of the result (29), when in a slightly ionized 
gas collisions between electrons and ions can be neglected. In this case F=1, 
«=o, and in the increase in dissipation is roughly in the ratio («,«,)"! = w,7,, 
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for free spiralling; which is Piddington’s result. However, in astronomical 
problems, even when the gas is only lightly ionized, it is normally not permissible 
to neglect the collisions between electrons and ions, because of their large 


clectrostatic collision-radius. 


8. Application to interstellar clouds.—'Vhe most immediate application of 
these results is to magnetic fields in interstellar clouds. An H1 cloud can be 
upposed to consist of neutral hydrogen and ionized metal atoms; for illustrative 
purposes we assume the latter to be 10~* of the whole number. ‘The pressure 
in an Hf cloud can perhaps be taken to be much the same as that in an H 11 cloud. 
Because of the lower temperature, this means that the density in an H1 cloud is 
a hundred or so times that in an H11 cloud; we assume a density of 100 atoms 
per cm*, and a temperature of 100 deg. K. 

In the absence of a magnetic field the conductivity would be about 
7x10 '*e.m.u.; in spite of the fewness of the charged particles it is mainly 
determined by the collisions of electrons with ions. ‘The time of decay of a 
magnetic field in a sphere of radius R is of order oR? if the conductivity o of the 
material is unaffected by the field. ‘hus the time of decay would be 6 x 10!” years 
for an H1 cloud of radius 5 pe if the effect of the field on the dissipation could 
be neglected. 

On the other hand, assuming a magnetic field of order 10 ° gauss, one finds 
values of x, «, and «, comparable with 3-5 x 10 and 2:3 x 10°®; the last 
figure does not greatly depend on the mass m, of the ions.* ‘hus the dissipation 
due to currents transverse to the magnetic field is, by equation (2g), increased 
roughly in the ratio 1-2 « 10'! asa result of the effect of the field on the conductivity. 
Since currents must flow transverse to the field for an appreciable part of their 
course, the time of decay of a field in the cloud considered above is reduced in 
about the same ratio, i.e. to 5 x 10" years. ‘These results should be compared 
with those of Gurevich (5), based on a less accurate formula for the reduced 
conductivity. 

A time of decay of order 5 x 10° years is still very large, but it is not, like a 
time of order 6 10!" years, so large as to make dissipation effects altogether 
negligible. ‘The calculated time of decay is admittedly very approximate. For 
example, «, «, and «, are inversely proportional to H; thus if the field were only 
10 ® gauss the time of decay would be increased 100-fold, making dissipation 
correspondingly less important. Likewise the assumed densities and collision- 
radii are somewhat uncertain. Perhaps the biggest ground for uncertainty 1s 
in a possible large increase in the dissipation by reason of turbulence, as suggested 
for different reasons by Sweet (6) and Elsasser (7). For these reasons, the possible 
effects of dissipation cannot be dismissed out of hand as certainly negligible in 
studies of galactic magnetic fields. 

‘The assumption that the whole of the material, neutral as well as ionized, 
moves about with the lines of force is still reasonably obeyed, even in the conditions 
assumed above. ‘lhe charged particles are closely attached to the lines of force 
by their spiral motions. ‘The motion of the neutral atoms relative to these 1s 
comparable with V,, where 

+k,). 


* Details of the calculation are given in the Appendix. 
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But if R is the radius of the cloud, ;~H/47R, since curlH = 47j; thus 
V, ~ Hi 4nneR(«, ~ 10 cm/sec. 


That is, the motion of neutral atoms across the lines of force is fairly slow 
compared with the mass motions actually likely to be present. However, the 
lines of force are less securely frozen into the material than is usual in astronomical 
problems. 

The figures quoted above refer to H1 clouds. For Ht clouds the situation 
is not altogether different, since even in these an appreciable fraction of the 
atoms can be expected to be neutral. ‘The significant term in the dissipation (27) 
is the term As compared with H1 clouds, «, +x, 
is reduced roughly in the ratio F/100 (being proportional to n,, the number 
density of neutral atoms); 7 is increased about 1oo-fold. ‘This shows that the 
dissipation is decreased roughly in the ratio F:1. ‘This even in Hu clouds the 
rate of dissipation need not be altogether insignificant. 


Note added 1956 February 
A referee has invited me to clear up a difficulty implicit in the expression (27 
for the dissipation of magnetic energy. ‘The expression indicates a dissipation 
which tends to infinity like p-* as the density p tends to zero, Such a result is 
clearly absurd; dissipation of electric currents occurs only because of collisions 
between ions, electrons and neutral atoms, and these become negligibly few 
as p->O. 

The dissipation indicated by the expression (27) is proportional to p * only 
if the current j is independent of p. In fact, however, in the case of importance 
In magneto-hydrodynamics jxH is at most comparable with gradp (cf. 
equation (g)). Thus jx H must be regarded as proportional to p, and the 
dissipation given by (27) remains finite at low densities. 

Even so, (27) overestimates the total dissipation at very low densities. 
Equation (10), on which it is based, is a diffusion equation, valid if the gas moves 
as a whole, the different constituents (electrons, ions, neutral atoms) having 
only relatively small velocities relative to the mass. In Section 8 this assumption 
was shown to be justified even in the rarefied interstellar medium. But if the 
density were further reduced, a point would be reached at which the assumption 
is no longer justified. According to equation (24), at low densities j, is proportional 
to j x H/p, and so is independent of the density; by equation (14), this implies 
that the velocity of diffusion V, of ions relative to the mass increases like 1/p. 
‘The increase cannot go on indefinitely. It goes on until V, becomes comparable 
with the mass-motions present; when this is so the process ceases to be one 
ot diffusion in a single gas, and instead there are a neutral gas and an ion gas moving 
independently. ‘The mutual velocity of the two gases remains finite; electro 
magnetic dissipation, which is due to collisions between particles of the two 
gases, now drops off to zero as p*, like the number of collisions. 

It follows, therefore, that the expression (27) does not apply at vanishingly 
small densities. However, even the low densities of interstellar clouds are 
sufficient to make it applicable. 
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Appendix 
‘The equation (2g) for o, can be written 
I I F?H 


(30) 


Under conditions of free spiralling the second term on the right is the important 
one; thus we need to determine «,and«,. On using equations (12), (17) and (19), 


Kk, =6,,/nell. 


ae/ 
Also by reference (2), equation (36), 6,, can be expressed in terms of the coethcient 
of diffusion D,, of electrons relative to neutral atoms in the absence of ions by 
the equation 


“  (n+n,)D,, 
‘The coefficient of diffusion D,, is given in terms of the effective colliston-distance 


r,, by 


(kT(m,+m,)\'* 3 


dD 
27m S(n+n,)r,,,” 


Combining these results 
3eH m, +m. (31) 
Similarly 


a’ at 


* {2mkTm,m, 12 
3eH | m,+m, 


(32) 


‘The values of «, and «, quoted in the text are obtained by taking 


x 10° cm, 10" cm, 
values suggested by mobility measurements. Since «, is large compared with «,, 
for purposes of calculating o, we do not need to know «, very accurately. The 
mass m, is taken to be that of a proton; «, is not sensitive to the exact value of m,, 
and is calculated assuming m,/m, large. ‘Vhe value of « can be calculated from 
a similar formula; the values of « and o, quoted in the text correspond to a 
collision-distance of about 4 x 10-° em for the (electrostatic) collisions between 
ions and electrons. 
The University 
Leeds 2: 
1955 December 1 
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